Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,171,417 G→A L850L (CTG→TTG)  mfd ← transcription‑repair coupling factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,171,4170GA97.6% 133.0 / NA 42L850L (CTG→TTG) mfdtranscription‑repair coupling factor
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (12/29);  minor base C (1/0);  total (13/29)
Fisher's exact test for biased strand distribution p-value = 3.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCTTCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCGC  >  NC_000913/1171170‑1171662
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
gCTTCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATtc                                                                                                                                                                                                                                                       <  1:28711/251‑1 (MQ=255)
                  gCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGccc                                                                                                                                                                                                                                     >  1:337260/1‑251 (MQ=255)
                  gCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTATACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGccc                                                                                                                                                                                                                                     <  1:253077/251‑1 (MQ=255)
                   cagcGGTAGTCTTCGCTTTTGGAGGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCg                                                                                                                                                                                                                                    <  2:138119/248‑1 (MQ=255)
                   cATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCg                                                                                                                                                                                                                                    >  1:464755/1‑251 (MQ=255)
                              cccGCTTTTGGATTGGGTGGCAGCAACGATGCATATGCCTGATGATGCGAACGTCCGACGCTACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGc                                                                                                                                                                                                                         <  2:333331/249‑1 (MQ=255)
                                        gATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCt                                                                                                                                                                                                              <  1:14197/251‑1 (MQ=255)
                                                         cATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCTGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGc                                                                                                                                                                                             <  1:468280/251‑1 (MQ=255)
                                                            gCATATGCCTGATGATGCGAACGTCCGTCGCGCCCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGAGCCGCGCGTTCAATGATAATAGTGTTGGCTGACGGGATGTCTATCCCGGTTTCGATAATGGTTGTACAATCCAGCACGTTGAAACGTTGATGATGGAAAGCATTCATCACCCGTGCCAATTCGCGCCCGCGCATCGGCCCGTGACCGATGGCGATCTGCGATTCTGGCACCAGTTCTGCCAg                                                                                                                                                                                          <  2:33476/251‑1 (MQ=255)
                                                              atatGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGGGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGcc                                                                                                                                                                                        <  2:464755/251‑1 (MQ=255)
                                                                           tgCGAACGTCCGACGTGACCCCGGAACTGGTGCATCTGCGCCTGACCTAAGTGAACAGCGGGTTCAAGGAGAATATGGTTGGCATTCGGGATGTCGATCCCGGTTTCGATAAGGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCCATTCGCGCTCGCGCATCGGCCCGTGACCGATGGCGACCCGCGCTGCTTGCACCAGTTCTGCCAGCCGTTCGGCGGCtt                                                                                                                                                                            <  2:195764/250‑1 (MQ=255)
                                                                              gAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGCTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTg                                                                                                                                                                        >  1:84151/1‑251 (MQ=255)
                                                                                aCGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCAATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGaa                                                                                                                                                                      <  1:3679/251‑1 (MQ=255)
                                                                                  gTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGTGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGc                                                                                                                                                                                             <  1:10802/226‑1 (MQ=255)
                                                                                   tCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGt                                                                                                                                                                   >  2:283355/1‑251 (MQ=255)
                                                                                        cgcgACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTcc                                                                                                                                                              >  2:76503/1‑251 (MQ=255)
                                                                                              cgcgTAACTTGAGCCGCGGCGCCAGACCGAAGTGATCCGCGCGTTTAAGTATAAGAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACTTTGAAACGTTGATGATGGAAATCATTCACCACCCGTTCCAATTCGCGCTCGCGCATCTGCCTGGGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCAcatca                                                                                                                                                        <  2:458554/251‑1 (MQ=255)
                                                                                               gcgTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCAcatcat                                                                                                                                                      <  2:91983/251‑1 (MQ=255)
                                                                                               gcgTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCAcatcat                                                                                                                                                      >  1:129441/1‑251 (MQ=255)
                                                                                               gcgTAACTGGGGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTTAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCAcatcat                                                                                                                                                      <  2:66935/251‑1 (MQ=255)
                                                                                                   aaCTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTa                                                                                                                                                 <  2:135437/251‑1 (MQ=255)
                                                                                                    cctggagcaggagcgcCAGACCGAAGTGATCCGCGCGTTCAATGCTAATAGTGTTGGATGCAGGGTTGTCGATCCCGGTTTCCAATATGGTTGTCCAAACACGCCCGTTGAAACGTTCATGATGTAAACCATTCATCACCCGTTCCACTTCGCGCTCACGGATCTGCCCGTGCCCGGTGGCGTTCCGCGCTTCTGGCCCCAGTTCTACCAGCCGTTCGGCGGCTTTTGGGATGTTTTCGCCATCAtttagt                                                                                                                                                 >  2:14211/13‑247 (MQ=255)
                                                                                                         tGCAGCTGCGCCAGACCGAAGTGATCCGCGCATTCACTGATAATAGTGGTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCCCGTTGAAACGTTGATTATGGAAATCATTAATCACCCGTTCCAATTCGCGCGCGCGCATCTGCCAGTGACCGATGGCGATCCGCGCTTCTGGCCCCAGGTCTTCCAGCAGTTCGGCGGCTTTCTGCATGTTTTCCACACACTTGTAGAGat                                                                                                                                            >  2:156609/1‑250 (MQ=255)
                                                                                                                gcgcCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATaaa                                                                                                                                      <  1:451356/249‑1 (MQ=255)
                                                                                                                     aGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGTAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACACCATTGTCGAGGTAATAACCCTGtcc                                                                                                                               >  2:5468/1‑251 (MQ=255)
                                                                                                                     aGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGGCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGtcc                                                                                                                               <  2:108202/251‑1 (MQ=255)
                                                                                                                       accgaatggATCCGCGTGTTAAATGAGAATAGTGTTGGCTGTCGGGATGTCGAGCCCGGTTTCGATAATGGTTGTACAAACCAGCACTTTGAAACGTTGAGGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGAGAATAAACCTGtcctc                                                                                                                             <  2:443943/243‑1 (MQ=255)
                                                                                                                              tGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGa                                                                                                                                                                                                                      >  1:71014/1‑157 (MQ=255)
                                                                                                                              tGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTGCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGa                                                                                                                                                                                                                      <  2:71014/157‑1 (MQ=255)
                                                                                                                                                                          gATCCCGGTTTCGATAATGGTTGTACAAACCAGCTCGGTGAAACGTTGATGTTGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACt                                                                          <  2:383608/251‑1 (MQ=255)
                                                                                                                                                                                         aatggtattACAAACCAGCACGTTGAAACGTTTATGATGGAAATCATCCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCTCGGACAAAGGttt                                                           <  2:212726/243‑1 (MQ=255)
                                                                                                                                                                                         aaTGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGttt                                                           <  1:230936/251‑1 (MQ=255)
                                                                                                                                                                                             gTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGTATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCGCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAAc                                                       <  1:23799/251‑1 (MQ=255)
                                                                                                                                                                                                      aCCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAgacg                                              <  2:223973/251‑1 (MQ=255)
                                                                                                                                                                                                            aCGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGAcgggc                                        <  1:283355/251‑1 (MQ=255)
                                                                                                                                                                                                                            gagggAAGTCATTCATCACCCGTTCCTATTCGCGCCCGCGCATCTGCCCGTGACCGATGGCTATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGTCGACGGGCGGTCGGCGTGGCGAt                         <  2:330399/247‑1 (MQ=255)
                                                                                                                                                                                                                                gAAATCATTCATCACCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCCGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCg                    >  1:277049/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaTCATTCATCACCCGTGCCAATTCGCGCTCGCTCAGCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGACAGCCGTTCGGCGGCTTTCGGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCAGGCGGTCATACGCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGTCGATAATCGat                  <  2:18395/251‑2 (MQ=255)
                                                                                                                                                                                                                                   aacattcctCCCCCGGTCCAATTCGCGCTCGCGCTGTTGCCCGGGACCGATGGCGGATCGCGGTTTTGGCACCCATTCTGGCCGCCGTTTGGCTGCCTTCTGAATGGTTTCCCCATTATTTTTAGGATAATAAACCCGTCCTCCGCGCAAAATTTCCCGCAGGAACGCCTCCCGGACCCCCAAGCTGGCAATCTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACa                 <  2:389334/249‑1 (MQ=255)
                                                                                                                                                                                                                                         catcaCCCGTTCCAATTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCa            <  1:394997/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 gTTCCAATTCGCGCACGCGCATATGCACGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCg   >  1:389047/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ttCCAATTCGCGCTGGCGCATCTGCCCCGGACCGATGGCGATCCGCGCTTCTGGCACCAGCGCTGCCAGCCGTTCGGCGGCTCTCTGAATGTTTTCCACATCATTGCAGAGCTAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGTCTCCCGGAACACCCTGCTGGCATACTCCCGGACAAAGGATTTAAATGCCCGCCGACGGGCGGGCGGCGTGGCGTTAATCGACAGGGCCCTCATTCCGc  >  2:157693/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                        
GCTTCAAGACGTTTTTGTGCATCGGTAGTCATCGCTTTTGGATGCGGTGTCAGCAACCATGCATATGCCTGATGATGCGAACGTCCGACGCGACCGCGTAACTGGTGCAGCTGCGCCAGACCGAAGTGATCCGCGCGTTCAATGATAATAGTGTTGGCTGTCGGGATGTCGATCCCGGTTTCGATAATGGTTGTACAAACCAGCACGTTGAAACGTTGATGATGGAAATCATTCATCACCCGTTCCAGTTCGCGCTCGCGCATCTGCCCGTGACCGATGGCGATCCGCGCTTCTGGCACCAGTTCTGCCAGCCGTTCGGCGGCTTTCTGAATGTTTTCCACATCATTGTAGAGATAATAAACCTGTCCTCCGCGCAAAATTTCACGCAGGATCGCCTCCCGGACCACCATGCTGTCATACTCACGGACAAAGGTTTTAACTGCCAGACGACGGGCGGGCGGCGTGGCGATAATCGACAGGTCACGCATTCCGC  >  NC_000913/1171170‑1171662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: