Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,103,105 C→A 6.7% pseudogene (285/446 nt) wbbL ← interrupted rhamnosyltransferase WbbL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,103,1050CA 6.7% 492.0 / 17.1 148pseudogene (285/446 nt)wbbLinterrupted rhamnosyltransferase WbbL
Reads supporting (aligned to +/- strand):  ref base C (69/68);  new base A (5/5);  total (74/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

Only 100 of 149 total aligned reads displayed.

ATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGC  >  NC_000913/2102956‑2103253
                                                                                                                                                     |                                                                                                                                                    
 tAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCt                                                                                                                                                     >  2:1264341/1‑150 (MQ=255)
 tAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCt                                                                                                                                                     >  1:1201316/1‑150 (MQ=255)
     tGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCat                                                                                                                                                 <  1:477133/150‑1 (MQ=255)
         acaATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgatg                                                                                                                                             >  1:1049579/1‑150 (MQ=255)
          caATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCAtgatga                                                                                                                                            >  2:1051193/1‑150 (MQ=255)
          caATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCAtgatga                                                                                                                                            >  2:759271/1‑150 (MQ=255)
            aTAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatt                                                                                                                                          <  1:73964/150‑2 (MQ=255)
              aaaaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATc                                                                                                                                        >  1:413647/1‑150 (MQ=255)
              aaaaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATc                                                                                                                                        <  2:1216721/150‑1 (MQ=255)
                aaTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATgatgatttgct                                                                                                                                      <  2:1270029/150‑6 (MQ=255)
                 aTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAATGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCTCCATATTTTATATATTTCAGCAATTCATCATGCTTCATGATGATATCggg                                                                                                                                     >  2:1095079/1‑150 (MQ=255)
                 aTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTTCTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCggg                                                                                                                                     >  2:1100339/1‑150 (MQ=255)
                 aTCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCggg                                                                                                                                     >  1:521095/1‑150 (MQ=255)
                  tCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGa                                                                                                                                    >  2:967174/1‑150 (MQ=255)
                   cAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGAt                                                                                                                                   <  2:61505/150‑1 (MQ=255)
                    aGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGAtt                                                                                                                                  <  2:1244662/150‑1 (MQ=255)
                     gAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTc                                                                                                                                 >  1:388572/1‑150 (MQ=255)
                      aaaGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCa                                                                                                                                <  1:1064600/150‑1 (MQ=255)
                          cacaGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCaaaaa                                                                                                                            <  2:1184931/150‑1 (MQ=255)
                           acacGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCTGGATTCAAAAAc                                                                                                                           <  2:716456/150‑1 (MQ=255)
                           acaGGAAATTTTCTTACTGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATAATCATGCTTCATGATGATATCGGCATTCaaaaaa                                                                                                                           >  2:264756/1‑149 (MQ=255)
                            caGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACa                                                                                                                          >  1:1107578/1‑150 (MQ=255)
                               gAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACaaaa                                                                                                                       <  1:635148/150‑1 (MQ=255)
                                aaaTTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAAt                                                                                                                      <  2:175062/150‑1 (MQ=255)
                                aaaTTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAAt                                                                                                                      <  2:801592/150‑1 (MQ=255)
                                aaaTTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAAt                                                                                                                      <  2:966362/150‑1 (MQ=255)
                                  aTTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGt                                                                                                                    >  1:1181563/1‑150 (MQ=255)
                                     ttCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAat                                                                                                                 <  2:159227/150‑1 (MQ=255)
                                       cTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatca                                                                                                               <  2:760110/150‑1 (MQ=255)
                                         tACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatc                                                                                                             <  2:981859/150‑1 (MQ=255)
                                         tACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAatcatc                                                                                                             <  1:169282/150‑1 (MQ=255)
                                                taatCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCg                                                                                                      <  1:967174/150‑1 (MQ=255)
                                                  atCATGTAAATATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCggg                                                                                                    <  1:916794/150‑1 (MQ=255)
                                                  atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGGATATCATAATGGTTGATGATGATATGGGTATTCAAAAACAAAATGTAATCATCATCTGgggg                                                                                                    <  2:430251/150‑1 (MQ=255)
                                                  atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCggg                                                                                                    <  2:1233130/150‑1 (MQ=255)
                                                  atCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCggg                                                                                                    >  2:145415/1‑150 (MQ=255)
                                                       gTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtat                                                                                               >  1:309102/1‑150 (MQ=255)
                                                        tAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCtata                                                                                              >  1:151728/1‑150 (MQ=255)
                                                         aaaGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCatgatg                                                                                                                                             >  2:906966/1‑102 (MQ=255)
                                                           aGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAAt                                                                                                                      <  2:221153/123‑1 (MQ=255)
                                                               ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTcc                                                                                       <  2:931201/150‑1 (MQ=255)
                                                               ttCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTcc                                                                                       <  2:384694/150‑1 (MQ=255)
                                                                tCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCt                                                                                      <  1:764593/150‑1 (MQ=255)
                                                                 cGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCtt                                                                                     <  1:1278619/150‑1 (MQ=255)
                                                                  gCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCttt                                                                                    <  1:1316530/150‑1 (MQ=255)
                                                                          tcGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGc                                                                            >  2:431284/1‑150 (MQ=255)
                                                                           cGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTATTCATCATCTGCGGGTCTATATTTTTCCt                                                                                      >  1:1116482/1‑139 (MQ=255)
                                                                              aaCAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATAATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCAcc                                                                        <  2:453017/150‑1 (MQ=255)
                                                                                cAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGc                                                                      <  2:521095/150‑1 (MQ=255)
                                                                                 agtcATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCa                                                                     <  1:1384645/147‑1 (MQ=255)
                                                                                  ggCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCaa                                                                    >  1:29296/1‑150 (MQ=255)
                                                                                   gCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATCTTTCCTTTACATACGCCACCGCAat                                                                   >  1:461736/1‑150 (MQ=255)
                                                                                    cATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAata                                                                  <  1:570039/150‑1 (MQ=255)
                                                                                     aTAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAatat                                                                 >  1:184136/1‑150 (MQ=255)
                                                                                        aTGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattat                                                              >  1:189914/1‑150 (MQ=255)
                                                                                         tGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGTATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                             <  2:337222/150‑1 (MQ=255)
                                                                                         tGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCCGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattatt                                                             <  1:759271/150‑1 (MQ=255)
                                                                                            aCTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGc                                                                      <  1:618537/138‑1 (MQ=255)
                                                                                            aCTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGc                                                                      >  2:618537/1‑138 (MQ=255)
                                                                                             cTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattatta                                                         <  2:1107578/150‑1 (MQ=255)
                                                                                              tAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAAtattattattat                                                        >  2:263731/1‑150 (MQ=255)
                                                                                               aaaaGCATAACTCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATg                                                       <  2:388572/150‑1 (MQ=255)
                                                                                               aaaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATCTTTCCTTTACATACGCCACCGCAATATTATTATTATg                                                       <  2:461736/150‑1 (MQ=255)
                                                                                               aaaaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGAATCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATg                                                       >  2:499298/1‑150 (MQ=255)
                                                                                                 aaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc                                                     <  1:1264341/150‑1 (MQ=255)
                                                                                                 aaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc                                                     <  1:1115370/150‑1 (MQ=255)
                                                                                                 aaGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc                                                     <  1:128601/150‑1 (MQ=255)
                                                                                                    cATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaa                                                  >  2:977789/1‑150 (MQ=255)
                                                                                                     aTAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCaaa                                                 >  2:358335/1‑150 (MQ=255)
                                                                                                      tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGTGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAg                                                >  1:1275525/1‑150 (MQ=255)
                                                                                                      tAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCAACCGCAATATTATTATTATGACCAAAg                                                >  1:1351227/1‑150 (MQ=255)
                                                                                                       aaCGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTTTTATTATGACCAAAGc                                               >  2:941318/1‑150 (MQ=255)
                                                                                                       aaCGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGc                                               <  1:596226/150‑1 (MQ=255)
                                                                                                         cGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg                                             >  1:1040683/1‑150 (MQ=255)
                                                                                                         cGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCg                                             >  1:1045720/1‑150 (MQ=255)
                                                                                                             tACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGtata                                         <  1:145415/150‑1 (MQ=255)
                                                                                                             tACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGtata                                         <  2:1303950/150‑1 (MQ=255)
                                                                                                              aCTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATac                                        <  2:1049579/150‑1 (MQ=255)
                                                                                                                   tgACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTc                                   <  1:264756/149‑1 (MQ=255)
                                                                                                                   cGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTc                                   <  2:413647/150‑1 (MQ=255)
                                                                                                                    gACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGc                                                                      <  1:1305351/114‑1 (MQ=255)
                                                                                                                    gACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGc                                                                      >  2:1305351/1‑114 (MQ=255)
                                                                                                                         atTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTaa                             >  1:1287266/1‑150 (MQ=255)
                                                                                                                          tttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAat                            <  2:1271341/150‑1 (MQ=255)
                                                                                                                           ttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAata                           <  2:1343921/150‑1 (MQ=255)
                                                                                                                           ttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAata                           >  2:839601/1‑150 (MQ=255)
                                                                                                                           ttAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAata                           >  2:878930/1‑150 (MQ=255)
                                                                                                                            tAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAatat                          >  2:490241/1‑150 (MQ=255)
                                                                                                                             aaTATATGTCATCAAATCATCATTCTTCATTATTATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAAtata                         <  2:599494/150‑1 (MQ=255)
                                                                                                                               tatatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGt                       >  1:826569/1‑150 (MQ=255)
                                                                                                                                 tatGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTcc                     >  1:179422/1‑150 (MQ=255)
                                                                                                                                    gTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAgg                  <  2:1181563/150‑1 (MQ=255)
                                                                                                                                      cAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGcc                <  2:689497/150‑1 (MQ=255)
                                                                                                                                      cAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGcc                >  1:1144017/1‑150 (MQ=255)
                                                                                                                                       ccgatatcatcATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc                                                     <  2:222360/107‑1 (MQ=255)
                                                                                                                                       ccgatatcatcATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGAc                                                     >  1:222360/6‑112 (MQ=255)
                                                                                                                                       ccgatatcatcATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCt               >  1:846233/6‑150 (MQ=255)
                                                                                                                                        gCAAATCATCATGCTTCATGATGATATCGGGAATCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTg              <  1:499298/150‑1 (MQ=255)
                                                                                                                                           tatcatcATGAAGCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCAt           >  1:780083/2‑150 (MQ=255)
                                                                                                                                              catcatGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAAt        >  2:43627/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAATATATGTCAGCAAATCATCATGCTTCATGATGATATCGGGATTCAAAAACAAAATGTAATCATCATCTGCGGGTCTATATTTTTCCTTTACATACGCCACCGCAATATTATTATTATGACCAAAGCCGTATACACCTCCACTAATATAGTCCAGGCCTGCATAATGCTGGC  >  NC_000913/2102956‑2103253

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