Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,196,803 A→C 5.9% intergenic (+430/+64) icd → / ← ymfD isocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,196,8030AC 5.9% 309.3 / 7.0 116intergenic (+430/+64)icd/ymfDisocitrate dehydrogenase; e14 prophage attachment site; tellurite reductase/e14 prophage; putative SAM‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (56/53);  new base C (5/2);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 4.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 117 total aligned reads displayed.

TACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAAC  >  NC_000913/1196658‑1196952
                                                                                                                                                 |                                                                                                                                                     
tACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAgcgc                                                                                                                                                   >  1:1035071/1‑150 (MQ=255)
 aCCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGcgcg                                                                                                                                                  >  1:845958/1‑150 (MQ=255)
  ccACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCgg                                                                                                                                                 >  2:1572874/1‑150 (MQ=255)
  ccACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCgg                                                                                                                                                 >  2:343320/1‑150 (MQ=255)
    aCTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGcc                                                                                                                                               >  1:368699/1‑150 (MQ=255)
      tCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGcccc                                                                                                                                             <  1:532286/150‑1 (MQ=255)
       acTTACACTCGCTCAAATATTGTTCAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCt                                                                                                                                            <  2:93918/149‑1 (MQ=255)
         ttACACCCGCTCAAATATTGTTAACTTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                                                                          <  1:1431576/150‑1 (MQ=255)
         atACCCCCGCTCAAAAATTGTTAAATTGCCGGGAGTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAACCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTc                                                                                                                                          <  2:1309068/149‑1 (MQ=255)
          tACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGCTCGCGGcctttttt                                                                                                                                         <  2:802545/150‑7 (MQ=255)
          aacacCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTca                                                                                                                                         <  1:1178267/149‑1 (MQ=255)
           acacCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTcac                                                                                                                                        <  2:1590321/150‑1 (MQ=255)
                cGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCacataca                                                                                                                                   >  2:192529/1‑150 (MQ=255)
                  cTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATc                                                                                                                                 <  1:854395/150‑1 (MQ=255)
                    cAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCtt                                                                                                                               >  2:1645229/1‑150 (MQ=255)
                     aaaTATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                                                                              <  2:190134/150‑1 (MQ=255)
                     aaaTATTGCTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCttt                                                                                                                              <  2:1184089/150‑1 (MQ=255)
                           tGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGTAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACt                                                                                                                        >  2:754809/1‑150 (MQ=255)
                            gTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGTGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTATAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                                                                       <  2:1168951/150‑1 (MQ=255)
                            gTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAGGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                                                                       >  2:1033960/1‑150 (MQ=255)
                            gTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTg                                                                                                                       <  2:1434631/150‑1 (MQ=255)
                             ttAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTga                                                                                                                      >  1:743688/1‑150 (MQ=255)
                               aaaTTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATTCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACAAACATCTTTAGTACTgaga                                                                                                                    >  1:189115/1‑150 (MQ=255)
                                      cGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTa                                                                                                             >  1:1038151/1‑150 (MQ=255)
                                       ggTTTTGTATCAACTACTCACCCGGGACTCGGCAGGGGACAGCCAACATGCATTGGGTGCAATCACCTTAGCGTTCATGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTaa                                                                                                            <  1:888612/150‑1 (MQ=255)
                                          tttGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCt                                                                                                         <  2:743688/150‑1 (MQ=255)
                                           ttGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTa                                                                                                        >  2:196981/1‑150 (MQ=255)
                                            tGTATCAACTAATCACCCGGGACTCTCCAGGGGAAAGCCAACAGGAATTGGGTGCAATCACCTTTCCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAg                                                                                                       <  1:1147998/150‑1 (MQ=255)
                                             gTATCAACTACTCACCCGGGACTCGCCAGGGGACAGTCTACAGGCATTGGGTGCAATCACCTTAGCGTTCATGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                                                                      <  1:954443/150‑1 (MQ=255)
                                             gTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                                                                      <  2:26845/150‑1 (MQ=255)
                                             gTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAgg                                                                                                      >  2:545207/1‑150 (MQ=255)
                                              tATCAACTACTCCCCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTACCCTAggg                                                                                                     >  2:445738/1‑150 (MQ=255)
                                                 cAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAat                                                                                                  >  2:11464/1‑150 (MQ=255)
                                                 cAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAat                                                                                                  <  1:1479638/150‑1 (MQ=255)
                                                 cAACTAATCACCCGGGACTCGCCATGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCATGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAat                                                                                                  <  2:1272484/150‑1 (MQ=255)
                                                    ctactCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAattat                                                                                               <  1:1464525/150‑1 (MQ=255)
                                                        tCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATc                                                                                           <  1:1051701/150‑1 (MQ=255)
                                                         cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAACGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                                                                          >  1:751854/1‑150 (MQ=255)
                                                         cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                                                                          >  1:387581/1‑150 (MQ=255)
                                                         cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                                                                          <  1:391734/150‑1 (MQ=255)
                                                         cACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCg                                                                                          >  2:1086490/1‑150 (MQ=255)
                                                           cccGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGta                                                                                        <  1:321469/150‑1 (MQ=255)
                                                              gggACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGtatta                                                                                     <  1:1696597/150‑1 (MQ=255)
                                                                  cTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATAtt                                                                                 >  2:851654/1‑150 (MQ=255)
                                                                     gCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCa                                                                              >  1:561165/1‑150 (MQ=255)
                                                                     gCCAGGGGACAGCCAACAGGCATTGGGGGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAATAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCa                                                                              <  1:339914/150‑1 (MQ=255)
                                                                       aaGGTGACAGCGAACAGGCATTTTGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAATAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTATACTAGGGAATTATTATCGTATTATATTGCAtt                                                                            <  1:461911/149‑1 (MQ=255)
                                                                         ggggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAACGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTAc                                                                          >  1:260155/1‑150 (MQ=255)
                                                                         ggggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTAc                                                                          <  1:192529/150‑1 (MQ=255)
                                                                          gggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACa                                                                         <  1:565656/150‑1 (MQ=255)
                                                                           ggACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACaa                                                                        <  1:1342383/150‑1 (MQ=255)
                                                                            gACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAAc                                                                       <  1:1193093/150‑1 (MQ=255)
                                                                              cAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTc                                                                     >  2:501622/1‑150 (MQ=255)
                                                                               aGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGAGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTcc                                                                    <  1:774048/150‑1 (MQ=255)
                                                                                gCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCAAGGGAATTATTATCGTATTATATTCGATTACAACTccc                                                                   >  2:314740/1‑150 (MQ=255)
                                                                                 ccAACAGGCATTGGGTGCAATCACCTTAGCGTTCATGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCa                                                                  <  1:186414/150‑1 (MQ=255)
                                                                                  cAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCaa                                                                 >  2:177641/1‑150 (MQ=255)
                                                                                   aaCAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAAc                                                                >  2:1683527/1‑150 (MQ=255)
                                                                                       ggCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg                                                            <  2:933698/150‑1 (MQ=255)
                                                                                       ggCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTg                                                            <  1:1392247/150‑1 (MQ=255)
                                                                                        gCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGTAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTAACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGa                                                           <  1:80185/150‑1 (MQ=255)
                                                                                        gCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGa                                                           <  1:970252/150‑1 (MQ=255)
                                                                                          tttGGGTGCAATCACCTTAGCGTTCAGGTACATGCTGAATGTAAATAAGGCCGCGAGCGCGGCCTCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACa                                                         <  2:845958/149‑1 (MQ=255)
                                                                                          aTTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACa                                                         >  1:902417/1‑150 (MQ=255)
                                                                                            tGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACata                                                       <  2:500264/150‑1 (MQ=255)
                                                                                             gggTGCAATCACCTTAGCTTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCCCAGACATCTTTAGGACTGCGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACattt                                                      >  2:89831/1‑148 (MQ=255)
                                                                                             gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACatat                                                      >  1:1160140/1‑150 (MQ=255)
                                                                                             gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACatat                                                      >  1:1066833/1‑150 (MQ=255)
                                                                                             gggTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACatat                                                      <  2:459240/150‑1 (MQ=255)
                                                                                                  cAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACa                                                               <  1:788319/136‑1 (MQ=255)
                                                                                                  cAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACa                                                               >  2:788319/1‑136 (MQ=255)
                                                                                                  cAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACt                                                 >  2:1074438/1‑150 (MQ=255)
                                                                                                  cAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACt                                                 <  1:1539726/150‑1 (MQ=255)
                                                                                                      cACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCt                                             <  1:343320/150‑1 (MQ=255)
                                                                                                        ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGc                                           >  1:1407108/1‑150 (MQ=255)
                                                                                                        ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGc                                           >  1:92251/1‑150 (MQ=255)
                                                                                                        ccTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATAAATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGc                                           >  1:296312/1‑150 (MQ=255)
                                                                                                               gTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTc                                    >  2:287200/1‑150 (MQ=255)
                                                                                                                       aCATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACtt                            <  2:359177/150‑1 (MQ=255)
                                                                                                                        cATGCGGAATGTAAAAATGGCCGCGAGCGCGGACCCTTCACATACATATTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACttt                           <  2:189115/150‑1 (MQ=255)
                                                                                                                         aTGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTa                          >  2:1135695/1‑150 (MQ=255)
                                                                                                                           gCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAat                        >  1:106288/1‑150 (MQ=255)
                                                                                                                           gCGGAATGTAAAAAAGGCCGCGAGCGCGGCCACTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGAATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACACCACGACTTTAat                        >  2:738493/1‑150 (MQ=255)
                                                                                                                            cGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAata                       >  2:1370761/1‑150 (MQ=255)
                                                                                                                                  gTAAAAAAGGCCGCGAGCGCAGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCaa                 >  2:1526641/1‑150 (MQ=255)
                                                                                                                                   tAAAAAAGGCCGCGAGCGCTGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATTTGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAg                <  1:851654/150‑1 (MQ=255)
                                                                                                                                   tAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAg                <  2:1343206/150‑1 (MQ=255)
                                                                                                                                       aggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTAAATATTTCAAGGCCa            >  1:1108325/4‑150 (MQ=255)
                                                                                                                                       aaaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCa            >  2:1375367/1‑150 (MQ=255)
                                                                                                                                        ggggCCGCGCTCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCAt           >  2:267034/3‑150 (MQ=255)
                                                                                                                                        aaGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTATGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCAt           <  1:842940/150‑1 (MQ=255)
                                                                                                                                         aGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCAtg          <  1:1645229/150‑1 (MQ=255)
                                                                                                                                           gCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCAtgtg        >  2:720332/1‑150 (MQ=255)
                                                                                                                                           gCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCAtgtg        >  1:293279/1‑150 (MQ=255)
                                                                                                                                            ccGCGAGCGCGGCCCCTTCATATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTga       >  2:1498772/1‑150 (MQ=255)
                                                                                                                                            ccGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTAATCCACGACTTTAATATTTCAAGGCCATGTga       >  1:507854/1‑150 (MQ=255)
                                                                                                                                             cgcgAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGCCTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTgaa      >  1:1219382/1‑150 (MQ=255)
                                                                                                                                               ggAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAAACTAGGGAATTATTATCGTTTTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACTACTTTAATATTTCAAGGCCATGTgaaga    <  2:1485508/149‑1 (MQ=255)
                                                                                                                                               cgAGCGCGGCCCCTTCACATACATCTTAAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTgaaga    >  1:1425356/1‑150 (MQ=255)
                                                                                                                                                 aGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAAc  >  2:504451/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TACCACTCCTTACACCCGCTCAAATATTGTTAAATTGCCGGTTTTGTATCAACTACTCACCCGGGACTCGCCAGGGGACAGCCAACAGGCATTGGGTGCAATCACCTTAGCGTTCAGGTACATGCGGAATGTAAAAAAGGCCGCGAGCGCGGCCCCTTCACATACATCTTTAGTACTGAGACTGTTTAACCTAGGGAATTATTATCGTATTATATTGCATTACAACTCCCAACAGTGACATATGAACTTCCTGCTTTACTCCACGACTTTAATATTTCAAGGCCATGTGAAGAAC  >  NC_000913/1196658‑1196952

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