Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,051,958 C→T Q2341* (CAA→TAA)  yeeJ → inverse autotransporter adhesin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,051,9580CT98.3% 397.9 / ‑2.8 115Q2341* (CAA→TAA) yeeJinverse autotransporter adhesin
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (59/54);  total (60/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 115 total aligned reads displayed.

CATCGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACAGAACTAGCCCTTC  >  NC_000913/2051822‑2052091
                                                                                                                                        |                                                                                                                                     
cATCGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTaaa                                                                                                                                    >  1:80115/1‑140 (MQ=255)
 aTCGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAg                                                                                                                                   <  2:677737/140‑1 (MQ=255)
  tCGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGc                                                                                                                                  >  1:730594/1‑140 (MQ=255)
   cGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCg                                                                                                                                 >  2:526652/1‑140 (MQ=255)
   cGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCg                                                                                                                                 <  2:80115/140‑1 (MQ=255)
         aCTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTg                                                                                                                           >  1:1005821/1‑140 (MQ=255)
             tACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTaa                                                                                                                                     <  2:1151583/126‑1 (MQ=255)
                aGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCt                                                                                                                    <  1:1135016/140‑1 (MQ=255)
                aGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCt                                                                                                                    >  1:850348/1‑140 (MQ=255)
                aGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCt                                                                                                                    <  1:598298/140‑1 (MQ=255)
                 gTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCtt                                                                                                                   <  1:610291/140‑1 (MQ=255)
                  tGAATGGGGAGCTATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTg                                                                                                                  >  2:1408118/1‑140 (MQ=255)
                    aaTGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTggg                                                                                                                >  2:869512/1‑140 (MQ=255)
                    aaTGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTggg                                                                                                                <  2:1430641/140‑1 (MQ=255)
                    aaTGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTggg                                                                                                                >  2:618224/1‑140 (MQ=255)
                       ggggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGaa                                                                                                                         <  1:334127/128‑1 (MQ=255)
                       ggggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGaa                                                                                                                         >  2:334127/1‑128 (MQ=255)
                       ggggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGc                                                                                                                     <  2:599784/132‑1 (MQ=255)
                       ggggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGc                                                                                                                     >  1:599784/1‑132 (MQ=255)
                       ggggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGttt                                                                                                             >  1:382332/1‑140 (MQ=255)
                       ggggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGtt                                                                                                              >  1:427761/1‑139 (MQ=255)
                         ggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGc                                                                                                           <  2:9816/140‑1 (MQ=255)
                         ggAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGc                                                                                                           <  2:1151707/140‑1 (MQ=255)
                            gaTATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTa                                                                                                        >  1:844282/1‑140 (MQ=255)
                              tatGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATg                                                                                                      >  2:1094503/1‑140 (MQ=255)
                                tGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCg                                                                                                    <  2:382332/140‑1 (MQ=255)
                                tGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCg                                                                                                    <  1:380389/140‑1 (MQ=255)
                                 ggggCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGa                                                                                                   <  2:234499/140‑1 (MQ=255)
                                   ggCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACa                                                                                                 <  2:1360759/140‑1 (MQ=255)
                                   ggCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACa                                                                                                 <  1:1148872/140‑1 (MQ=255)
                                   ggCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACa                                                                                                 <  2:1080886/140‑1 (MQ=255)
                                   ggCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACa                                                                                                 >  1:1080886/1‑140 (MQ=255)
                                    gCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACAt                                                                                                >  2:26282/1‑140 (MQ=255)
                                      aTTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGa                                                                                                   <  1:720622/135‑1 (MQ=255)
                                       ttACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAGCAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGtt                                                                                             >  1:954406/1‑140 (MQ=255)
                                        tACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTta                                                                                            >  2:72027/1‑140 (MQ=255)
                                               cTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAAc                                                                                     >  1:1422222/1‑140 (MQ=255)
                                                  aCGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCtg                                                                                  >  1:1343608/1‑140 (MQ=255)
                                                  aCGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCtg                                                                                  >  2:741593/1‑140 (MQ=255)
                                                       ggTTTTCAATCAAATATGTATTGGTTATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGAtt                                                                             >  2:1290805/1‑140 (MQ=255)
                                                         ttttCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGAtttt                                                                           >  2:228871/1‑140 (MQ=255)
                                                                 caaATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTtata                                                                   <  1:1094503/140‑1 (MQ=255)
                                                                   aaTATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAAt                                                                 <  2:321955/140‑1 (MQ=255)
                                                                   aaTATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAAt                                                                 <  1:26282/140‑1 (MQ=255)
                                                                        gTATTGGTTATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCtt                                                            <  1:1290805/140‑1 (MQ=255)
                                                                        gTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGACTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCtt                                                            <  1:1291297/140‑1 (MQ=255)
                                                                          aTTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGGTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTaa                                                          >  2:1282280/1‑140 (MQ=255)
                                                                           ttGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTaaa                                                         <  2:985748/140‑1 (MQ=255)
                                                                           ttGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTaaa                                                         <  1:945097/140‑1 (MQ=255)
                                                                              gTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAagag                                                      >  1:1282759/1‑140 (MQ=255)
                                                                              gTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAagag                                                      >  1:606088/1‑140 (MQ=255)
                                                                              gTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATATTATAAAAACCTGTGATTTTCCGTTATAATTTCtt                                                            <  1:413316/134‑1 (MQ=255)
                                                                              gTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATATTATAAAAACCTGTGATTTTCCGTTATAATTTCtt                                                            >  2:413316/1‑134 (MQ=255)
                                                                              gTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACAAGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAagag                                                      >  2:1375757/1‑140 (MQ=255)
                                                                              gTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGCGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAagag                                                      >  2:1117392/1‑140 (MQ=255)
                                                                               tCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAaga                                                       >  2:1120103/1‑138 (MQ=255)
                                                                               tCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAaga                                                       <  1:1120103/138‑1 (MQ=255)
                                                                                   cTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCtt                                                            >  1:779451/1‑129 (MQ=255)
                                                                                   cTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCtt                                                            <  2:779451/129‑1 (MQ=255)
                                                                                   cTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAaga                                                       <  2:1374932/134‑1 (MQ=255)
                                                                                   cTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAaga                                                       >  1:1374932/1‑134 (MQ=255)
                                                                                   cTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTc                                                 >  2:516996/1‑140 (MQ=255)
                                                                                   cTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTc                                                 >  1:247755/1‑140 (MQ=255)
                                                                                      gTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCAt                                              >  2:1166198/1‑140 (MQ=255)
                                                                                       tCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCAtt                                             <  1:526652/140‑1 (MQ=255)
                                                                                        cccGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTc                                            >  2:957510/1‑140 (MQ=255)
                                                                                            cAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCgggg                                        >  2:1377740/1‑140 (MQ=255)
                                                                                            cAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCgggg                                        <  1:72027/140‑1 (MQ=255)
                                                                                            cAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCgggg                                        <  1:869512/140‑1 (MQ=255)
                                                                                            cAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCgggg                                        >  1:365314/1‑140 (MQ=255)
                                                                                             aaaCTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGc                                       >  2:153151/1‑140 (MQ=255)
                                                                                              aaCTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCa                                      <  2:155756/140‑1 (MQ=255)
                                                                                              aaCTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCa                                      <  2:730594/140‑1 (MQ=255)
                                                                                               aCTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCAc                                     <  2:300045/140‑1 (MQ=255)
                                                                                               aCTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCAc                                     <  1:385215/140‑1 (MQ=255)
                                                                                                cTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCAct                                    <  1:1408118/140‑1 (MQ=255)
                                                                                                   aaGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCtt                                 >  2:119281/1‑140 (MQ=255)
                                                                                                    aGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCttt                                <  1:957510/140‑1 (MQ=255)
                                                                                                         aCAATACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTAc                           >  1:602567/1‑140 (MQ=255)
                                                                                                           aaTACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACAt                         <  2:1422222/140‑1 (MQ=255)
                                                                                                           aaTACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACAt                         <  1:153151/140‑1 (MQ=255)
                                                                                                             tACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCtt                                 >  2:792986/1‑130 (MQ=255)
                                                                                                             tACGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTc                       >  1:655644/1‑140 (MQ=255)
                                                                                                              aCGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTcc                      <  2:1343608/140‑1 (MQ=255)
                                                                                                              aCGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTcc                      <  1:741593/140‑1 (MQ=255)
                                                                                                              aCGTAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTcc                      <  1:937209/140‑1 (MQ=255)
                                                                                                                 tAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACAtt                        <  1:1053720/135‑1 (MQ=255)
                                                                                                                 tAGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACAtt                        >  2:1053720/1‑135 (MQ=255)
                                                                                                                  aGTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCctct                  >  1:564801/1‑140 (MQ=255)
                                                                                                                   gTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTa                 <  1:1166198/140‑1 (MQ=255)
                                                                                                                   gTTTCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTa                 <  2:1393937/140‑1 (MQ=255)
                                                                                                                     ttCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACa               >  2:492424/1‑140 (MQ=255)
                                                                                                                     ttCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACa               >  2:948860/1‑140 (MQ=255)
                                                                                                                      tCCCTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACAg              >  2:851058/1‑140 (MQ=255)
                                                                                                                         cTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTcc                      <  1:1464328/129‑1 (MQ=255)
                                                                                                                         cTGGCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTcc                      >  2:1464328/1‑129 (MQ=255)
                                                                                                                           ggCAACAGGTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACAGAACTa         <  1:250887/140‑1 (MQ=255)
                                                                                                                                  ggTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACAGAACTAGCCCTtc  >  2:801776/1‑140 (MQ=255)
                                                                                                                                   gTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTa                 >  2:820542/1‑124 (MQ=255)
                                                                                                                                   gTGATTAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTa                 <  1:820542/124‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CATCGGGTCACTGTACAGTGAATGGGGAGATATGGGGCATTACACGACTGACGCTGGTTTTCAATCAAATATGTATTGGTCATCTAGTCCCGCAAACTCAAGCGAACAATACGTAGTTTCCCTGGCAACAGGTGATCAAAGCGTATTTGAAAAGCTTGGGTTTGCTTATGCGACATGTTATAAAAACCTGTGATTTTCCGTTATAATTTCTTAAAGAGTGTTCCATTCGGGGCACTCTTTTTTACATTCCCTCTACAGAACTAGCCCTTC  >  NC_000913/2051822‑2052091

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