Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,972,674 G→T 6.1% E51D (GAG→GAT) ‡ rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,972,6740GT 6.1% 296.5 / 7.6 99E51D (GAG→GAT) ‡rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base G (45/48);  new base T (4/2);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

Only 100 of 108 total aligned reads displayed.

AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAACC  >  NC_000913/3972525‑3972821
                                                                                                                                                     |                                                                                                                                                   
agaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGAAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgag                                                                                                                                                     >  2:887783/1‑150 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    <  2:1355978/150‑1 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    <  2:1675883/150‑1 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    >  2:1284606/1‑150 (MQ=255)
 gaAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGATAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAagcgagc                                                                                                                                                    <  1:1403862/150‑1 (MQ=255)
    aaaTTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCt                                                                                                                                                 >  2:214881/1‑150 (MQ=255)
     aaTTCTGATAACAGGTGGTGCCGGGTTTATTTGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCtt                                                                                                                                                >  2:662370/1‑150 (MQ=255)
        tCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGc                                                                                                                                             <  1:466922/150‑1 (MQ=255)
        tCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCCGGCACCGGTCGCGCAAAGCGAGCGCTTTGc                                                                                                                                             <  1:1235516/150‑1 (MQ=255)
         cTGATAACAGGTGGTGGCGGGTTTATTGGCGAGGCGCGGGAGCGTTATATCATCATCTAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGTCACCGGTCGCGCAAAGCGAGCGCTTTGcc                                                                                                                                            <  1:1362901/150‑1 (MQ=255)
         cTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcg                                                                                                                                            <  1:828054/150‑2 (MQ=255)
         cTGATAACAGGGGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCTCGCTTTGcg                                                                                                                                            <  1:507790/150‑2 (MQ=255)
          tGATAAAAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCt                                                                                                                                           >  1:1393525/1‑150 (MQ=255)
           gATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCtt                                                                                                                                          >  2:654215/1‑150 (MQ=255)
            aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                                                                                                         <  2:1351502/150‑1 (MQ=255)
            aTAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCttt                                                                                                                                         <  2:1426026/150‑1 (MQ=255)
              aagaGGTGGTGCCGGGTTTATTGGCTCGGCGCGGGTGCGTTATAACATCCACGAAACGAGCGACGGGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                       <  1:380664/147‑1 (MQ=255)
              aaCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATTTCATCAACGACACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTga                                                                                                                                       <  1:497089/150‑1 (MQ=255)
               aCAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCAGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTgag                                                                                                                                      >  2:1241678/1‑150 (MQ=255)
                cAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTTCCGTTgaga                                                                                                                                     >  2:1665764/1‑150 (MQ=255)
                cAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTgaga                                                                                                                                     <  1:214881/150‑1 (MQ=255)
                    gggtgCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGt                                                                                                                                 <  1:373781/149‑1 (MQ=255)
                      gtgACGGGTTTATTGGCTTGGCGCTGGTGCGTTATATCATCTAGTAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCTAAGCGTGCGCTTTGCCTTTGAGAAAGTTg                                                                                                                               <  1:589417/150‑1 (MQ=255)
                         ccGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGata                                                                                                                            <  1:810362/150‑1 (MQ=255)
                          cGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGCCGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGatat                                                                                                                           >  2:1180056/1‑150 (MQ=255)
                             gTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTg                                                                                                                        >  2:732857/1‑150 (MQ=255)
                             gTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTg                                                                                                                        <  2:826140/150‑1 (MQ=255)
                              tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGc                                                                                                                       <  1:1627892/150‑1 (MQ=255)
                              tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGc                                                                                                                       >  1:506964/1‑150 (MQ=255)
                              tttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGc                                                                                                                       <  1:409287/150‑1 (MQ=255)
                               ttATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCg                                                                                                                      <  2:940782/150‑1 (MQ=255)
                                 aTTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGAt                                                                                                                    <  2:1210472/150‑1 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTCCGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCGCGCTTTTCCTTTGAGAAAGTTGATATATGCGATCgg                                                                                                                 >  2:1166084/1‑150 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCgg                                                                                                                 >  2:1142423/1‑150 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCgg                                                                                                                 >  2:385584/1‑150 (MQ=255)
                                    ggCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGAGCGCGGGTTGCCTTTGAGAAAGTTGATATCTGCGATAgg                                                                                                                 >  2:596911/1‑150 (MQ=255)
                                      cTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGc                                                                                                               <  1:494080/150‑1 (MQ=255)
                                           cgcTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAAc                                                                                                          <  1:902114/150‑1 (MQ=255)
                                             gtGGTGCGGTATATCAGCAACGAAACGCGCGACGCGGTGCTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGGCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTg                                                                                                        <  1:332502/149‑1 (MQ=255)
                                             cTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTg                                                                                                        <  2:248244/150‑1 (MQ=255)
                                              tGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTgg                                                                                                       <  2:1617525/150‑1 (MQ=255)
                                               ggTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGc                                                                                                      <  2:833661/150‑1 (MQ=255)
                                                    gTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCAcgcg                                                                                                 <  2:797645/150‑1 (MQ=255)
                                                         atcatcAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTc                                                                                            >  2:1030925/1‑150 (MQ=255)
                                                          tcatcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCa                                                                                           <  1:443937/150‑1 (MQ=255)
                                                            atcaACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACt                                                                                         <  1:887783/150‑1 (MQ=255)
                                                             tcaACGCAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTg                                                                                        <  2:70046/150‑1 (MQ=255)
                                                               gacgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAg                                                                                      <  2:1432246/149‑1 (MQ=255)
                                                                acgaaacgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGAGCAATGCCATTGACAAAGTTGATATCAGCGATCGGGCAGAAATGGCACTCCTATTCACTGAGc                                                                                     >  2:947578/1‑150 (MQ=255)
                                                                      cgaGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCATTGCCGTTGAGAAAGTTTATATATGCGATCGGGGCGAACTGGCCCGCGTTTTAACTGAGCATCAGc                                                                               >  1:1661718/1‑150 (MQ=255)
                                                                        aGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCa                                                                             <  2:850292/150‑1 (MQ=255)
                                                                          cGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGa                                                                           <  2:163668/150‑1 (MQ=255)
                                                                            aCGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCCAAGCGCGCCCGTTGCCGTTTAGAAAGTTGAAATGAGCGAACGGGCAGAACAGGCACGCGGTATAACTGAGCATAAGCCCGAcg                                                                         >  1:1084542/1‑149 (MQ=255)
                                                                               cgGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACtgtg                                                                      <  1:368719/150‑1 (MQ=255)
                                                                                  tggtAGTGGTCGATAAGCAGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCTGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCa                                                                   <  1:131957/150‑1 (MQ=255)
                                                                                    gtAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATg                                                                 >  2:267847/1‑150 (MQ=255)
                                                                                     tAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGc                                                                <  1:1108421/150‑1 (MQ=255)
                                                                                          gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTg                                                           <  2:356006/150‑1 (MQ=255)
                                                                                          gTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCATTGCCATTGAGAAAGTTGATATCTGCGATCGGGCAGAACAGGCACCCGTATTCACTTACCATAAGCCAGACTGTGTGATTCATCTg                                                           >  2:786060/1‑150 (MQ=255)
                                                                                            cGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTCATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGc                                                         >  2:894708/1‑150 (MQ=255)
                                                                                                 aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGc                                                                >  1:1602219/1‑138 (MQ=255)
                                                                                                 aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGc                                                                <  2:1602219/138‑1 (MQ=255)
                                                                                                 aGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCg                                                    >  2:411360/1‑150 (MQ=255)
                                                                                                       ccTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGCGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACGGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGcc                                              >  2:818362/1‑150 (MQ=255)
                                                                                                             ccGGAAACCTGATGTCGCTGGCTCCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGTTCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCATCCGAAAGCCATGTTg                                        <  2:1273023/150‑1 (MQ=255)
                                                                                                             ccGGAAACCTGATGTCGCTGGCACCGGTAGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTg                                        >  1:230058/1‑150 (MQ=255)
                                                                                                               ggAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGTGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGAc                                      <  2:1611640/150‑1 (MQ=255)
                                                                                                               ggAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCGTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGAc                                      >  1:741664/1‑150 (MQ=255)
                                                                                                               ggAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCGTTGCCGTTGAGAAAGTTGATATATGCGCTCGGGCAGAACAGGCACGCGGATTAACTGAGCATCAGCCAGACTGTGGCATGCATCAGGCCGCCGCAAGCCATGTTGAc                                      >  2:1445287/1‑150 (MQ=255)
                                                                                                                 aaaCCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAATTTGATATATGCGATCGGGCCGAACTGGCACGCGGATTCACAGAGCAGCAGCCCGCCAGTGTGCTGCCTATGGCAGCCGAAAGCCATGTTTACCg                                    >  2:1253802/1‑150 (MQ=255)
                                                                                                                    ccTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTc                                 >  2:1373221/1‑150 (MQ=255)
                                                                                                                      tGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTa                               <  1:1284606/150‑1 (MQ=255)
                                                                                                                       gATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTAt                              >  1:901752/1‑150 (MQ=255)
                                                                                                                          gTCGCTGGCACCGGTCGCGCAAAACGAGCGCTTTGCCGTTTAGAAAGTTGATATCTGCGATCGGGGAGAACTTGCCCGCGTATTAACTGACCAACAGCCAGACAGTGTCAAGCATCTCCAGCCGAAACCCATGTTGACCCTTCTATTGAc                          >  2:387546/1‑150 (MQ=255)
                                                                                                                            cGCTGGCCCCGGTCGCGCAAAGCGAGCGGTTTGCCTGTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCCCGGAGCATCCGCCAGGCGGTGTCATGCATCTGGCAGCCGAAAGGCCGGTTGGCCGTTCGATTGACg                         >  2:127283/1‑150 (MQ=255)
                                                                                                                            cGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCCCGCGTATTCACTGAGCATCCGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACg                         >  1:641073/1‑150 (MQ=255)
                                                                                                                              cTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCCCGCGTATTCACTGAGCATAAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGCCCGTTCTATTGACGGc                       >  2:435241/1‑150 (MQ=255)
                                                                                                                              cTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGc                       <  1:1117792/150‑1 (MQ=255)
                                                                                                                               tGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGTCGGcc                      <  1:1241678/150‑1 (MQ=255)
                                                                                                                                ggCACCGGTCGCGCAAAGCGAGCGCGTTTCCGTTGAGAAATTTTATATCATCGATCGGGCAGAACAGGCACCCGGATTAACAGAGCATCAGCCCAACAGTGTCAAGCCTATGGCAGCCGAAAGCCAAGTTGACCGTTCTAATGACGGccc                     >  2:657593/1‑150 (MQ=255)
                                                                                                                                  cACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                   >  1:24430/1‑150 (MQ=255)
                                                                                                                                  cACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                   <  1:1142423/150‑1 (MQ=255)
                                                                                                                                  cACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCgg                   <  2:1260089/150‑1 (MQ=255)
                                                                                                                                      tcaaaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATTTTGACCGTGATATTGACGGCCCGGCAGc               >  2:591024/7‑150 (MQ=255)
                                                                                                                                       gTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCg              <  1:1373221/150‑1 (MQ=255)
                                                                                                                                         cgcgCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGtt            >  1:206657/1‑150 (MQ=255)
                                                                                                                                         aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGACATGTTGACCGTTCTATTGACGGCCCGGCTGCGtt            >  1:1488178/4‑150 (MQ=255)
                                                                                                                                         aaggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGAACGCGTATTCACAGAGCATCAGCCAGACTGTGTCATGCATCTCGCAGCCGAAAGCCATGTTGACCGTTCTATAGACGGCCCGGCAGCGtt            >  1:862276/4‑150 (MQ=255)
                                                                                                                                          gcgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGttt           >  2:103458/1‑150 (MQ=255)
                                                                                                                                           ggcAAAGCGCTCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGTGTTTa          >  1:1653723/2‑150 (MQ=255)
                                                                                                                                           cgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTa          >  2:651657/1‑150 (MQ=255)
                                                                                                                                           cgcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTa          >  1:772383/1‑150 (MQ=255)
                                                                                                                                            gcAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAt         >  1:999252/1‑150 (MQ=255)
                                                                                                                                             cAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAtt        <  2:769080/150‑1 (MQ=255)
                                                                                                                                             cAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTAtt        >  1:1320552/1‑150 (MQ=255)
                                                                                                                                               aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGa      <  1:1166084/150‑1 (MQ=255)
                                                                                                                                               aaGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGAAGCCGAAAGCCATGTTGACCGTTCTATTGACTGCCCGGCAGCGTTTATTGa      >  2:339850/1‑150 (MQ=255)
                                                                                                                                                 gcgagcgCTTTGCCTTTGTGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGaaa    <  1:878470/150‑1 (MQ=255)
                                                                                                                                                 gcgagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGaaa    <  1:1665764/150‑1 (MQ=255)
                                                                                                                                                   gagcgCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAAcc  <  1:290056/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AGAAAAATTCTGATAACAGGTGGTGCCGGGTTTATTGGCTCGGCGCTGGTGCGTTATATCATCAACGAAACGAGCGACGCGGTGGTAGTGGTCGATAAGCTGACCTACGCCGGAAACCTGATGTCGCTGGCACCGGTCGCGCAAAGCGAGCGCTTTGCCTTTGAGAAAGTTGATATCTGCGATCGGGCAGAACTGGCACGCGTATTCACTGAGCATCAGCCAGACTGTGTCATGCATCTGGCAGCCGAAAGCCATGTTGACCGTTCTATTGACGGCCCGGCAGCGTTTATTGAAACC  >  NC_000913/3972525‑3972821

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