Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,033,044 G→A N12N (AAC→AAT yedJ ← putative HD superfamily phosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,033,0440GA99.0% 370.5 / ‑4.2 97N12N (AAC→AATyedJputative HD superfamily phosphohydrolase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (53/43);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

TCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACC  >  NC_000913/2032896‑2033191
                                                                                                                                                    |                                                                                                                                                   
tCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGGTTTTCAACCAAtt                                                                                                                                                   >  1:1173745/1‑151 (MQ=255)
 cacaATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAAttt                                                                                                                                                  <  2:475770/151‑1 (MQ=255)
 cacaATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAAttt                                                                                                                                                  <  2:11118/151‑1 (MQ=255)
 cacaATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAAttt                                                                                                                                                  <  1:1276996/151‑1 (MQ=255)
 cacaATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAAtt                                                                                                                                                   <  2:933215/150‑1 (MQ=255)
  acaATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCCTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGGTTTTCAACCAAtttt                                                                                                                                                 >  1:1221613/1‑151 (MQ=255)
   caATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTc                                                                                                                                                <  2:1404348/151‑1 (MQ=255)
    aaTATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCg                                                                                                                                               >  1:1127059/1‑151 (MQ=255)
     atatCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAAtttt                                                                                                                                                 >  1:1644561/1‑148 (MQ=255)
     atatCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGGGCTGATGATGGTTTTTCAAACAATTTTCg                                                                                                                                               >  1:1580755/1‑150 (MQ=255)
         cGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCCCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACtg                                                                                                                                          >  1:883175/1‑151 (MQ=255)
          gTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACtgt                                                                                                                                         <  2:733928/151‑1 (MQ=255)
             aaaaTAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgc                                                                                                                                       <  2:597722/150‑1 (MQ=255)
              aaaTAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgcct                                                                                                                                     <  1:1468695/151‑1 (MQ=255)
              aaaTAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgcct                                                                                                                                     <  1:1228767/151‑1 (MQ=255)
              aaaTAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgcc                                                                                                                                      >  2:475240/1‑150 (MQ=255)
               aaTAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgcctg                                                                                                                                    <  2:827301/151‑1 (MQ=255)
                 tAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCCTCTTCCGCCCGCTTCTTTGCCGTAGCCCAGAAGCGGCGAAAATGGCACAAATCATGGGCCGCGTCCGGGTGCTTATGGTGGATTTTCAACCAAATTTCGAACTTtgcctgcc                                                                                                                                  >  1:1206097/1‑151 (MQ=255)
                 tAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgcctgcc                                                                                                                                  >  1:228527/1‑151 (MQ=255)
                  aaCAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGtgcctgcc                                                                                                                                  >  1:94523/1‑150 (MQ=255)
                      ggCGGTTAAAATCCCCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATTATGGTTTTTCACCCAATTTTCGAACTGTGCCTGCCAgt                                                                                                                               <  1:1241601/149‑1 (MQ=255)
                      ggCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAgtgt                                                                                                                             <  1:658227/151‑1 (MQ=255)
                       gCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGtgtt                                                                                                                            <  2:335713/151‑1 (MQ=255)
                        cGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGtgttg                                                                                                                           <  2:1385893/151‑1 (MQ=255)
                                aTCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTcc                                                                                                                    >  1:457366/1‑150 (MQ=255)
                                  cACCAGCATATCGACATCGTCATCTGCCGCCTGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGACATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGACCTGTGCCTGCCAGTGTTGTAAGTCCAt                                                                                                                  >  1:484167/1‑150 (MQ=255)
                                    ccAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATAcg                                                                                                               >  1:889587/1‑151 (MQ=255)
                                         atatCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCt                                                                                                           <  1:176015/150‑1 (MQ=255)
                                         atatCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTc                                                                                                          >  2:1236485/1‑151 (MQ=255)
                                         atatCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTc                                                                                                          <  2:1656875/151‑1 (MQ=255)
                                           atCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTcc                                                                                                         <  1:1334129/150‑1 (MQ=255)
                                           atCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCg                                                                                                        >  1:903275/1‑151 (MQ=255)
                                           atCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCg                                                                                                        >  1:104544/1‑151 (MQ=255)
                                            tCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGt                                                                                                       >  2:628270/1‑151 (MQ=255)
                                             cGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACTCGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGt                                                                                                       <  2:909568/150‑1 (MQ=255)
                                             cGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGtt                                                                                                      <  2:872918/151‑1 (MQ=255)
                                              gACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTg                                                                                                     >  2:1614676/1‑151 (MQ=255)
                                                 aTCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGcc                                                                                                   >  1:355247/1‑150 (MQ=255)
                                                 aTCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGcc                                                                                                   >  1:569241/1‑150 (MQ=255)
                                                        tCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTc                                                                                           >  1:1492133/1‑151 (MQ=255)
                                                         cTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAagag                                                                                             <  1:1066793/148‑1 (MQ=255)
                                                         cTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCg                                                                                          >  2:422424/1‑151 (MQ=255)
                                                          tGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGc                                                                                         >  1:460193/1‑151 (MQ=255)
                                                          tGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGc                                                                                         >  2:1098160/1‑151 (MQ=255)
                                                          tGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGc                                                                                         >  2:174666/1‑151 (MQ=255)
                                                           gccgccAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCg                                                                                          <  1:1698593/149‑1 (MQ=255)
                                                           gccgccAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGcc                                                                                        >  1:1229582/1‑151 (MQ=255)
                                                            ccgccAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCa                                                                                       <  2:1173251/151‑1 (MQ=255)
                                                              gccAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGt                                                                                            <  1:870311/144‑1 (MQ=255)
                                                                cAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCtt                                                                                   <  2:1234811/151‑1 (MQ=255)
                                                                 aGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCtt                                                                                   >  2:448093/1‑150 (MQ=255)
                                                                  gCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCtt                                                                                   >  2:1501750/1‑149 (MQ=255)
                                                                         tgCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGtttt                                                                          >  2:86691/1‑151 (MQ=255)
                                                                                 cccAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATa                                                                  <  2:945187/151‑1 (MQ=255)
                                                                                 cccAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATa                                                                  <  2:491679/151‑1 (MQ=255)
                                                                                   cAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGc                                                                <  2:412228/151‑1 (MQ=255)
                                                                                       cgcgGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCt                                                            <  2:457404/151‑1 (MQ=255)
                                                                                        gcggcgAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTg                                                           <  2:1708951/151‑1 (MQ=255)
                                                                                         cggcgAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGc                                                          >  1:390082/1‑151 (MQ=255)
                                                                                         cggcgAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGc                                                          >  1:855605/1‑151 (MQ=255)
                                                                                             gAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACg                                                      >  1:463393/1‑151 (MQ=255)
                                                                                              aaaaTGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGt                                                     <  1:968930/151‑1 (MQ=255)
                                                                                                  tGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCgg                                                 >  1:1699171/1‑151 (MQ=255)
                                                                                                    acacacGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGa                                               <  1:1670707/151‑1 (MQ=255)
                                                                                                     cacacGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGAc                                              >  1:197159/1‑151 (MQ=255)
                                                                                                      acacGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACg                                             <  1:1173129/151‑1 (MQ=255)
                                                                                                       cacGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGa                                            <  1:1610029/151‑1 (MQ=255)
                                                                                                          gTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAgg                                          >  2:1400405/1‑150 (MQ=255)
                                                                                                             aTGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCgg                                      >  2:476650/1‑151 (MQ=255)
                                                                                                               gggCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCggg                                     <  2:1524561/150‑1 (MQ=255)
                                                                                                               gggCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGgt                                    <  2:1653718/151‑1 (MQ=255)
                                                                                                                  ccGCGTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTa                                 <  2:818455/151‑1 (MQ=255)
                                                                                                                    gcgTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTAcc                               >  2:619242/1‑151 (MQ=255)
                                                                                                                     cgTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTAcc                               >  1:747261/1‑150 (MQ=255)
                                                                                                                     cgTCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTAcc                               >  1:117326/1‑150 (MQ=255)
                                                                                                                       tCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCa                              <  1:1525178/149‑1 (MQ=255)
                                                                                                                       tCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAgg                            >  1:273442/1‑151 (MQ=255)
                                                                                                                       tCCTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAgg                            >  2:1538775/1‑151 (MQ=255)
                                                                                                                        ccTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGc                           >  1:1028201/1‑151 (MQ=255)
                                                                                                                         cTGGTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGc                           >  1:70220/1‑150 (MQ=255)
                                                                                                                            gTGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAAtg                       <  2:584950/151‑1 (MQ=255)
                                                                                                                             tGCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAAtgc                      >  1:231188/1‑151 (MQ=255)
                                                                                                                              gCTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAAtgct                     <  1:1271809/151‑1 (MQ=255)
                                                                                                                               cTGATGATGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAAtgctg                    >  2:427489/1‑151 (MQ=255)
                                                                                                                                tgatgaTGGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAAtgctgc                   >  2:619921/1‑151 (MQ=255)
                                                                                                                                     atgGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTaa               >  2:1350858/1‑150 (MQ=255)
                                                                                                                                     atgGTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAAc              >  2:409241/1‑151 (MQ=255)
                                                                                                                                        gTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAAt           <  2:563621/151‑1 (MQ=255)
                                                                                                                                        gTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAAt           <  1:199714/151‑1 (MQ=255)
                                                                                                                                        gTTTTTCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAAt           >  2:904035/1‑151 (MQ=255)
                                                                                                                                          ttttCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATg          <  2:569292/150‑1 (MQ=255)
                                                                                                                                           tttCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGcc        >  1:1205581/1‑151 (MQ=255)
                                                                                                                                            ttCAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGccc       >  2:327541/1‑151 (MQ=255)
                                                                                                                                              cAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGcccc      <  1:422390/150‑1 (MQ=255)
                                                                                                                                              cAACCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCg     >  2:1519933/1‑151 (MQ=255)
                                                                                                                                                aCCAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGAc   >  2:1040986/1‑151 (MQ=255)
                                                                                                                                                 ccAATTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGAcc  <  1:483057/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCACAATATCGTGAAAATAACAGGCGGTTAAAATCACCAGCATATCGACATCGTCATCTGCCGCCAGCTTCTGTGCCGTAGCCCAGACGCGGCGAAAATGACACACGTCATGGGCCGCGTCCTGGTGCTGATGATGGTTTTTCAACCAGTTTTCGAACTGTGCCTGCCAGTGTTGTAAGTCCATACGCGCTCCGTTGCCAAAGAGTCGCCACCTTAGCAGTTTTTATTCATAGCGTCTGCTACGTGCGGGACGAGGTCGGGTGTACCAGGCAATGCTGCCTAACAATGCCCCGACC  >  NC_000913/2032896‑2033191

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