Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,993,582 G→A 5.9% T375T (ACC→ACT uvrC ← excinuclease UvrABC, endonuclease subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,993,5820GA 5.9% 191.8 / 2.4 70T375T (ACC→ACTuvrCexcinuclease UvrABC, endonuclease subunit
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base A (2/2);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.99e-01

ACTCCGCACGCAGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGAGCAGCGTT  >  NC_000913/1993434‑1993729
                                                                                                                                                    |                                                                                                                                                   
aCTCCGCACGCAGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACg                                                                                                                                                     <  2:761940/149‑1 (MQ=255)
 cTCCGCACGCAGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGCTGTCAAAGCACTCCCTCCGCTTCACTGCCGGCAATTTCAGCACGCTGGCACGCGCGGTCAGTCGCTGGTGACCgg                                                                                                                                                    >  1:1547642/1‑149 (MQ=255)
  tCCGCACGCAGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGt                                                                                                                                                   <  1:1101646/149‑1 (MQ=255)
   ccGCACGCAGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCTCTGGTGAACGGTa                                                                                                                                                  <  2:182106/149‑1 (MQ=255)
    cgcacgcaGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCCCTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAg                                                                                                                                                 >  2:1406084/1‑149 (MQ=255)
         gcaGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATtgc                                                                                                                                            >  2:642412/1‑149 (MQ=255)
           aGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCCCTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATtgctg                                                                                                                                          >  1:217427/1‑149 (MQ=255)
                gCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGaaa                                                                                                                                     <  2:1501991/149‑1 (MQ=255)
                 ccGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTACGGCAATTTCAACACGCTGGCAAGCGCGGTCAGGCGCTGGTGAACGGTAGATTGCTGCGAAAg                                                                                                                                    >  2:646266/1‑149 (MQ=255)
                    ttAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCCACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACAGTAGATTGCTGCGAAAGttt                                                                                                                                 >  1:227849/1‑149 (MQ=255)
                    ttAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGttt                                                                                                                                 >  2:1186392/1‑149 (MQ=255)
                    ttAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACAGTAGATTGCTGCGAAAGttt                                                                                                                                 >  1:221885/1‑149 (MQ=255)
                      aGCATCAAACCCCACACAGGAAGCGACGGTTTGGTCGCCCATGGTATGGCTGATGTCAAAGCACTACATCCGCTTCACTTCCGGCAATTTCAACCCGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGc                                                                                                                               <  1:1020323/149‑1 (MQ=255)
                       gCATCAAACACCACACAGGACGCGACTGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGAGAAAGTTTGCt                                                                                                                              <  2:817195/149‑1 (MQ=255)
                       gCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCt                                                                                                                              <  2:259959/149‑1 (MQ=255)
                        cATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTg                                                                                                                             >  2:548814/1‑149 (MQ=255)
                         aTCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTgg                                                                                                                            >  2:355136/1‑149 (MQ=255)
                            aaaCACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTa                                                                                                                         >  1:154593/1‑149 (MQ=255)
                             aaCACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTaa                                                                                                                        >  1:145616/1‑149 (MQ=255)
                               caccacACAGGAAGCGACGGTTTGTTCGCCCATGGTATGTCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAgg                                                                                                                      <  2:984947/149‑1 (MQ=255)
                                accacACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGc                                                                                                                     >  1:437478/1‑149 (MQ=255)
                                  cacacaGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCg                                                                                                                   >  2:1370178/1‑149 (MQ=255)
                                   acacaGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGt                                                                                                                  >  1:1158752/1‑149 (MQ=255)
                                                          gCCCATGGTATGGCTTATGCCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGtt                                                                                           <  1:213815/149‑1 (MQ=255)
                                                          gCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGtt                                                                                           >  1:1241083/1‑149 (MQ=255)
                                                           ctcATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCCCTTCCGGCAATTTCAACCCGCTGGCAAGCGCGGTCAGTCGCTGGTCAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGgtt                                                                                          >  2:1412344/3‑149 (MQ=255)
                                                            ccATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTc                                                                                         <  2:1129677/149‑1 (MQ=255)
                                                                  tATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATaa                                                                                   >  2:356365/1‑149 (MQ=255)
                                                                    tGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAAAGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAAcg                                                                                 >  2:1193344/1‑149 (MQ=255)
                                                                        ggATGTCAAAGCACTGCATCCGCTGCACTTCCGGCACTTTCACCTCGCTGGCAAGCGCGGTCAGGCGCTGGTGAACGGGAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGcc                                                                             <  2:1361432/148‑1 (MQ=255)
                                                                         gATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGccc                                                                            >  1:966245/1‑149 (MQ=255)
                                                                            gTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTATATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGACCTAt                                                                         >  2:365654/1‑149 (MQ=255)
                                                                             tCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATc                                                                        >  1:1277933/1‑149 (MQ=255)
                                                                               aaaGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATcgc                                                                      >  1:434830/1‑149 (MQ=255)
                                                                                 aGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGACAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGCTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATcgccg                                                                    >  2:554718/1‑149 (MQ=255)
                                                                                  gCACTCCATCCGCTTCACTTCCGGCAAATTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATcgccgc                                                                   <  1:311327/149‑1 (MQ=255)
                                                                                    aCTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATcgccgcta                                                                 <  1:172944/149‑3 (MQ=255)
                                                                                    aCTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGGCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGa                                                                 <  2:937511/149‑1 (MQ=255)
                                                                                        cATCCCCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCCTCGCCGCATTGTTGAGCGCGAGTGTCAGATAACGCTCCCTATCGCCGCGAggtt                                                             >  2:160502/1‑149 (MQ=255)
                                                                                          tCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGGCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAggtttg                                                           <  1:642412/149‑1 (MQ=255)
                                                                                            cGCTTCACATCCGGCAACTTCAACACGCTGCCAAGCGCGGTCAGTCCCTGGTGAACGGACGATTGCTGCGAAAGTTTGCTGGTTACGGCCGGCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAggtttggt                                                         <  2:1483439/149‑1 (MQ=255)
                                                                                             gCTTCACTTCCGGCAATTTCAACACGCTGGCATGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTAGATGGTTAAGGCCGTCGCAGCATTGTTGCGCGCGAGTTTCAGATAACGCGACGTATCGACGCGAggtttggtt                                                        >  1:107492/1‑149 (MQ=255)
                                                                                                tCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGACGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGa                                                     >  2:1068910/1‑149 (MQ=255)
                                                                                                     tCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACAtt                                                <  1:548814/149‑1 (MQ=255)
                                                                                                           actttCACCACGCAGGCAAGCGCGGCCAGTCGCTGGGGAACCGTAGATTGCTGCGAAAGTTAGCTGGTTACGGCCGGCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTACTCtt                                          <  2:887612/147‑1 (MQ=255)
                                                                                                           aaTTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCtt                                          <  1:1186392/149‑1 (MQ=255)
                                                                                                           aaTTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCtt                                          <  2:437478/149‑1 (MQ=255)
                                                                                                             tttCAACACGCTGGCAAGCGCGGGCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGc                                        <  1:1406084/149‑1 (MQ=255)
                                                                                                              ttCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCg                                       >  2:820518/1‑149 (MQ=255)
                                                                                                              ttCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCg                                       >  1:1372682/1‑149 (MQ=255)
                                                                                                              ttCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCg                                       >  1:984637/1‑149 (MQ=255)
                                                                                                                cAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTc                                     >  1:1068268/1‑149 (MQ=255)
                                                                                                                    acGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGc                                 <  1:1370178/149‑1 (MQ=255)
                                                                                                                    acGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGGTAGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGAGCCCTATCGCCGGGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGc                                 >  2:271599/1‑149 (MQ=255)
                                                                                                                          gCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTc                           >  1:1548298/1‑149 (MQ=255)
                                                                                                                           cAAGCGCGGTCAGGCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCGGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGGCCCGCCAGTTCt                          <  2:1411910/149‑1 (MQ=255)
                                                                                                                            aaGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTg                         <  1:18767/149‑1 (MQ=255)
                                                                                                                             aGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGGCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGa                        <  2:715707/149‑1 (MQ=255)
                                                                                                                               cgcgGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGaaa                      <  2:69707/149‑1 (MQ=255)
                                                                                                                                 cgGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAgg                    <  1:355136/149‑1 (MQ=255)
                                                                                                                                  ggTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAggg                   >  2:1025308/1‑149 (MQ=255)
                                                                                                                                    tCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTGTGAAAGGGaa                 >  2:36732/1‑149 (MQ=255)
                                                                                                                                      aGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATAAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATc               >  1:451838/1‑149 (MQ=255)
                                                                                                                                        tCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCgg             >  1:867205/1‑149 (MQ=255)
                                                                                                                                          gCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCTCCGCATTGGTGCGCGCGAGTTTCAGATAACGAGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGACAGTTCTGAAAGGGAATCGGCg           >  2:1019878/1‑149 (MQ=255)
                                                                                                                                              gTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTGTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGagcc       >  1:1053606/1‑148 (MQ=255)
                                                                                                                                                gAACGGTAGATTGCTCCGAAAGTTTGCTTGTTAAGGCCGTCGCCGCAATGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGagcagc     <  2:644784/149‑1 (MQ=255)
                                                                                                                                                gAACAGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGagcagc     <  2:227849/149‑1 (MQ=255)
                                                                                                                                                gAACAGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGagcagc     <  2:221885/149‑1 (MQ=255)
                                                                                                                                                  aCGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGAGCAGCGt   <  2:143410/149‑1 (MQ=255)
                                                                                                                                                   cGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGAGCAGCGtt  >  1:148080/1‑149 (MQ=255)
                                                                                                                                                   cGGGAGATTGCTGCGAAAGTTTGCTGGTAAAGGCCGGCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGAGCAGCGtt  <  2:1241083/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACTCCGCACGCAGCGGGCCGTTAGCATCAAACACCACACAGGAAGCGACGGTTTGTTCGCCCATGGTATGGCTGATGTCAAAGCACTCCATCCGCTTCACTTCCGGCAATTTCAACACGCTGGCAAGCGCGGTCAGTCGCTGGTGAACGGTAGATTGCTGCGAAAGTTTGCTGGTTAAGGCCGTCGCCGCATTGGTGCGCGCGAGTTTCAGATAACGCGCCCTATCGCCGCGAGGTTTGGTTTGAACATTAATCTTGCGTCCCGCCAGTTCTGAAAGGGAATCGGCGAGCAGCGTT  >  NC_000913/1993434‑1993729

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: