Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,165,154 T→A 7.0% M274L (ATG→TTG)  parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,165,1540TA 7.0% 223.4 / 5.8 71M274L (ATG→TTG) parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base A (3/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCG  >  NC_000913/3165005‑3165303
                                                                                                                                                     |                                                                                                                                                     
ctTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCAt                                                                                                                                                       >  2:1874289/1‑150 (MQ=255)
    tCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgc                                                                                                                                                   <  1:1468003/150‑1 (MQ=255)
     ccAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgcg                                                                                                                                                  >  2:985252/1‑150 (MQ=255)
      cAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTgcgc                                                                                                                                                 <  2:1285806/150‑1 (MQ=255)
       aGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGCTGCGCAAATGCgca                                                                                                                                                <  1:202716/150‑1 (MQ=255)
        gATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCgcag                                                                                                                                               <  1:1183025/150‑1 (MQ=255)
        gATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCgcag                                                                                                                                               <  1:1190558/150‑1 (MQ=255)
              tggtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAAttt                                                                                                                                         <  2:301934/150‑1 (MQ=255)
               ggtAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGCTGCGCAAATGCGcaacaaaaag                                                                                                                                        <  2:235810/150‑5 (MQ=255)
                  aGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCgcagc                                                                                                                                              >  2:320886/1‑141 (MQ=255)
                  aGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCgcagc                                                                                                                                              <  1:320886/141‑1 (MQ=255)
                   gCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTcc                                                                                                                                    >  2:507715/1‑150 (MQ=255)
                   gCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTcc                                                                                                                                    <  2:1426930/150‑1 (MQ=255)
                     gAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAg                                                                                                                                  <  1:924654/150‑1 (MQ=255)
                     gAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAg                                                                                                                                  <  1:674807/150‑1 (MQ=255)
                     gAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAg                                                                                                                                  <  2:1168050/150‑1 (MQ=255)
                            tggtTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGGAGATCGTCAACCATCTGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTAcgcgc                                                                                                                           <  2:60999/150‑1 (MQ=255)
                             ggtTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTAcgcgcg                                                                                                                          >  2:1524141/1‑150 (MQ=255)
                                tcatcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAc                                                                                                                       <  2:1511556/150‑1 (MQ=255)
                                tcatcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAc                                                                                                                       >  2:783143/1‑150 (MQ=255)
                                 catcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAcc                                                                                                                      >  2:430734/1‑150 (MQ=255)
                                 catcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAcc                                                                                                                      <  1:1524141/150‑1 (MQ=255)
                                 catcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAcc                                                                                                                      >  2:1487512/1‑150 (MQ=255)
                                 catcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAcc                                                                                                                      <  2:642229/150‑1 (MQ=255)
                                 catcaCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCAcc                                                                                                                      >  1:714427/1‑150 (MQ=255)
                                     aCCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGaa                                                                                                                  >  1:1074806/1‑150 (MQ=255)
                                           tCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGt                                                                                                                                 >  2:39493/1‑129 (MQ=255)
                                           tCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGt                                                                                                                                 <  1:39493/129‑1 (MQ=255)
                                           tCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGa                                                                                                            >  1:751567/1‑150 (MQ=255)
                                            ccATATCAACGCGGTTGGAACGCGGAACAAACACCAGGCTCGACGTGTTCTCCTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGAt                                                                                                           >  2:1137340/1‑150 (MQ=255)
                                             cATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCGGTACGCGCGCACCTGAAACCTGATg                                                                                                          <  2:20762/150‑1 (MQ=255)
                                              atatCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGc                                                                                                         >  2:155442/1‑150 (MQ=255)
                                              atatCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGc                                                                                                         <  2:714427/150‑1 (MQ=255)
                                                    aCGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCTGGTTCTCGTGGTCATATTCATCTCGCAGATCGTCAACCATCGGCAGCTTTTTGTTTCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAAt                                                                                                   <  2:733313/150‑1 (MQ=255)
                                                       cgGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTACAGTACGCGCGCACCTGAAACCTGATGCGGCAATgcg                                                                                                >  1:1797397/1‑150 (MQ=255)
                                                             gAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATa                                                                                          <  2:16058/150‑1 (MQ=255)
                                                             gAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATa                                                                                          >  2:1285979/1‑150 (MQ=255)
                                                              aaCGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATaa                                                                                         <  1:1874289/150‑1 (MQ=255)
                                                                    gCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACcg                                                                                   >  2:1566000/1‑150 (MQ=255)
                                                                      acaATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACcgcg                                                                                 <  1:1696575/150‑1 (MQ=255)
                                                                        aaTCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGcc                                                                               <  2:688212/150‑1 (MQ=255)
                                                                           cACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCAtc                                                                            >  2:257952/1‑150 (MQ=255)
                                                                               aGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAAttt                                                                                                                                         <  2:1146250/85‑1 (MQ=255)
                                                                               aGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAAttt                                                                                                                                         >  1:1146250/1‑85 (MQ=255)
                                                                                    ggtcgggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTtctttctt                                                                   >  2:635193/1‑150 (MQ=255)
                                                                                      gcgggtTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTcc                                                                 <  1:307309/149‑1 (MQ=255)
                                                                                              tcgtcgTCAGAATCATCGCGCAGATCCTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCCCAGAAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACcgcg                                                         >  2:325175/1‑150 (MQ=255)
                                                                                                 tGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTtctttctt                                                                   >  2:376633/1‑137 (MQ=255)
                                                                                                 tGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTtctttctt                                                                   <  1:376633/137‑1 (MQ=255)
                                                                                                  ggTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCa                                                     <  1:430734/150‑1 (MQ=255)
                                                                                                          tCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACt                                             >  2:1216631/1‑150 (MQ=255)
                                                                                                          tCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGAAGCAATTTGCTCCAGTACGCTCGCAGGTGAAACATGATGCGGCAATTGGCTGATAACCACCTCGCCATCTTATTTCTTCCACACCGCTCGCATACGCACt                                             <  2:1963459/150‑1 (MQ=255)
                                                                                                              cgcgCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCCCCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAAc                                         >  2:1803030/1‑150 (MQ=255)
                                                                                                              cgcgCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAAc                                         >  1:1519139/1‑150 (MQ=255)
                                                                                                                 gcAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCAc                                      <  1:1137340/150‑1 (MQ=255)
                                                                                                                  cAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACg                                     <  1:985252/150‑1 (MQ=255)
                                                                                                                       cGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGt                                >  1:1265054/1‑150 (MQ=255)
                                                                                                                          cAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTtct                             <  1:635193/150‑1 (MQ=255)
                                                                                                                             ccATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGt                          >  2:288819/1‑150 (MQ=255)
                                                                                                                                   gCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGAttt                    <  2:751567/150‑1 (MQ=255)
                                                                                                                                      gCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTAc                 >  1:1289891/1‑150 (MQ=255)
                                                                                                                                       cTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACg                >  2:1955247/1‑150 (MQ=255)
                                                                                                                                        tttttGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACgg               <  1:1518051/150‑1 (MQ=255)
                                                                                                                                        aaattgctgCGCAAATGCGCAGCAATTTGCTCCAGTACGCGCGCAGCTGAAACCTGATGCGGCACTGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACgg               >  2:1408483/4‑150 (MQ=255)
                                                                                                                                        aaattgctgCGCAAATGCGCAACAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACgg               >  2:1944529/4‑150 (MQ=255)
                                                                                                                                         ttttGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGa              <  1:1549509/150‑1 (MQ=255)
                                                                                                                                         ttttGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGa              <  1:1487512/150‑1 (MQ=255)
                                                                                                                                         ttttGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGa              <  1:783143/150‑1 (MQ=255)
                                                                                                                                          attgctgCGCAAATGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAt             >  1:1580213/2‑150 (MQ=255)
                                                                                                                                           ttgttgCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGAtc            <  1:149548/150‑1 (MQ=255)
                                                                                                                                                gcgcATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgc       >  1:1290291/1‑150 (MQ=255)
                                                                                                                                                gcgcATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGgc       >  2:1764300/1‑150 (MQ=255)
                                                                                                                                                     tttGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGcgcgcg  <  1:155442/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTTTTCCAGATCGGTGGTAGCGAAGAGGTGGTTCATCACCTGATCCATATCCACGCGGTTGGAACGCGGCACAATCACCAGGCGGGTCGGGTTCTCGTGGTCAGATTCATCGCGCAGATCGTCAACCATCGGCAGCTTTTTGTTGCGCATTTGCGCAGCAATTTGCTCCAGTACGCGCGCACCTGAAACCTGATGCGGCAATGCGCTGATAACCACCGCGCCATCTTCTTTCTTCCACACCGCGCGCATACGCACTGAACCACGTCCGTTCTCGTAGATTTTACGGATCTCGGCGCGCG  >  NC_000913/3165005‑3165303

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