Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,966,862 G→A 24.2% R149H (CGT→CAT)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,8620GA24.2% 253.6 / 76.0 132R149H (CGT→CAT) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (43/57);  new base A (13/19);  total (56/76)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 132 total aligned reads displayed.

GATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAAC  >  NC_000913/3966716‑3967010
                                                                                                                                                  |                                                                                                                                                    
taTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTTAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                                                                                                                                   <  1:2096603/149‑1 (MQ=255)
taTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                                                                                                                                   <  1:1594512/149‑1 (MQ=255)
taTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                                                                                                                                   <  1:676679/149‑1 (MQ=255)
gATTCGCCCGCCGAAAGAAGTTGAACGCTATTTTGCGCTGCTTAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAACCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                                                                                                                                   <  1:638359/150‑1 (MQ=255)
gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTTAAAGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAAACCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTgg                                                                                                                                                   <  1:629787/150‑1 (MQ=255)
gATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATgg                                                                                                                                                   >  1:816647/1‑150 (MQ=255)
 aatCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGTAGTTAACTTAGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt                                                                                                                                                  <  1:334294/148‑1 (MQ=255)
 aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt                                                                                                                                                  <  1:2023600/150‑1 (MQ=255)
 aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGt                                                                                                                                                  <  1:1583180/150‑1 (MQ=255)
 aTTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGt                                                                                                                                                  <  1:2028849/150‑1 (MQ=255)
   tCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTaa                                                                                                                                                <  1:1917109/150‑1 (MQ=255)
    cgcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGCTAAc                                                                                                                                               >  1:698746/1‑150 (MQ=255)
     gcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACg                                                                                                                                              >  1:447719/1‑150 (MQ=255)
     gcccgccGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCCGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACg                                                                                                                                              >  2:1590613/1‑150 (MQ=255)
       ccgccgAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAATTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGt                                                                                                                                            <  1:586190/150‑1 (MQ=255)
       acgccgAAAGAAGGTGAACTCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAATATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGt                                                                                                                                            <  1:307044/149‑1 (MQ=255)
          ccgAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCt                                                                                                                                         <  2:1655532/150‑1 (MQ=255)
             aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                      >  2:462756/1‑150 (MQ=255)
             aaaGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACt                                                                                                                                      <  1:2006312/150‑1 (MQ=255)
              aagaagGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTg                                                                                                                                     >  1:1250738/1‑150 (MQ=255)
                  agGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGa                                                                                                                                 <  2:107341/150‑1 (MQ=255)
                   ggTGAACGCTAATTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAACGCTCGTCTGCGTATGTAACGTGGTAACGGTTCTACTGAAGAt                                                                                                                                <  1:1210404/150‑1 (MQ=255)
                        aCGCTATTTTGCGCTGCTGAAAGTTAAGGAAGTTAACTTCGACAAACCTGAAAACGCCCTCAACAAAATCCTCTTTGAGAACTTAACCCCGCGGCACGCAAACTCTCGTCTGCGTATGGAACGTGGAAACGGTTCTACTGAAGATTTAAc                                                                                                                           <  1:1743746/150‑1 (MQ=255)
                        aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                           <  1:1907112/150‑1 (MQ=255)
                        aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                           >  1:802676/1‑150 (MQ=255)
                        aCGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                           <  1:1338587/150‑1 (MQ=255)
                        aCGCTATTTTGAGATGCTGAAAGTTAACTAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAAc                                                                                                                           <  1:671643/150‑1 (MQ=255)
                                 tGCGCTGCTGAAAGTTAACGAAGTTAAATTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCAGCACGCAAACTCTCGTCAGCGTATGGAACATGGTAACGCTTCTACTGAAGAATTAACTGCTCGCGt                                                                                                                  >  1:1318363/1‑150 (MQ=255)
                                  gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTTACGGTTCTACTGAAGATTTAACTGCTCGCGTa                                                                                                                 >  2:223883/1‑150 (MQ=255)
                                  gcgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTa                                                                                                                 <  2:1079092/150‑1 (MQ=255)
                                   cgcTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTAc                                                                                                                >  2:588934/1‑150 (MQ=255)
                                      tgctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTgg                                                                                                             <  1:1446346/150‑1 (MQ=255)
                                       gctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAAAGGTTCTAATGAAGATTTAACTGCTCGCGTAATGGa                                                                                                            >  1:1177094/1‑150 (MQ=255)
                                        ctgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGAt                                                                                                           >  1:1140623/1‑150 (MQ=255)
                                         tgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGGGGTACAGGTTCTATTGAAGATTTATCTGCTCGCGTACTGGATc                                                                                                          >  2:1335875/1‑150 (MQ=255)
                                         tgAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAAAGCCCGCAACAAAATCCTCTATGAGAACTTAACACCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTACTACTGAAGATTTAACTGCTCGCGTACTGTATc                                                                                                          >  1:780108/1‑150 (MQ=255)
                                         tcgaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAACCTCTTTGAGAACTTAACCCCGCTGCACGCACACTCGCGTCTGCGAATGGAACATGGATACGGTTCTACTGAAGATTTAACTGCTCGCGTCATGGATc                                                                                                          >  2:1831960/4‑150 (MQ=255)
                                          gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                         >  2:1161611/1‑150 (MQ=255)
                                          gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                         >  2:1922757/1‑150 (MQ=255)
                                          gAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                         <  2:674647/150‑1 (MQ=255)
                                          cgaaGTTAACTTCGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCt                                                                                                         >  1:479532/3‑150 (MQ=255)
                                            aaGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTgg                                                                                                       >  1:1853772/1‑150 (MQ=255)
                                              gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                                                                     <  1:835007/150‑1 (MQ=255)
                                              gTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCa                                                                                                     <  2:251382/150‑1 (MQ=255)
                                                tAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTACCCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATc                                                                                                   <  1:174161/150‑1 (MQ=255)
                                                     aaGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGAATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATAACCTa                                                                                              >  2:942062/1‑150 (MQ=255)
                                                       gTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATc                                                                                            <  1:225283/150‑1 (MQ=255)
                                                         tAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCgg                                                                                          <  1:882169/150‑1 (MQ=255)
                                                         tAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCgg                                                                                          >  1:1059702/1‑150 (MQ=255)
                                                          aaCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCggg                                                                                         >  2:1883396/1‑149 (MQ=255)
                                                           aCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                        >  1:445359/1‑150 (MQ=255)
                                                           aCTTCGACAAACCTGAAAACGCCCGCAACAAAATCATCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTATGCGTATGGAACGTGGTAACGGTTTTAATGAAGATTTAAGTGATCGCGTACTGGATCTGGCATCACCTATCGGTc                                                                                        <  2:399416/150‑1 (MQ=255)
                                                            cTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCg                                                                                       <  1:658901/150‑1 (MQ=255)
                                                             ttCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACACAAACTCTCGTCTGCGTATGGAACGGGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCgg                                                                                      >  2:1040400/1‑149 (MQ=255)
                                                               cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                                                                    <  2:1345527/150‑1 (MQ=255)
                                                               cGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTAATGAAGTTTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTgg                                                                                    <  1:31073/150‑1 (MQ=255)
                                                                 aCAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTc                                                                                  <  2:526567/150‑1 (MQ=255)
                                                                       cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTg                                                                            >  1:1941202/1‑150 (MQ=255)
                                                                       cTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCGATCGGTCGTGGGCAGCGTg                                                                            >  1:677344/1‑150 (MQ=255)
                                                                          aaaaCGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTc                                                                         <  1:1416344/150‑1 (MQ=255)
                                                                             aCGCCAGCAACAAAATGATCTTTGAGAACGTAACCCCGCTGCACGCAAACTCTCGTCTGCGTAGTGAACGTGGTAACGGTTCTACTGAAGATTTAACTTCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGa                                                                      <  2:341795/150‑1 (MQ=255)
                                                                                 ccGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCCTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATtgt                                                                  <  2:1653366/150‑1 (MQ=255)
                                                                                      acaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAc                                                             >  2:952215/1‑150 (MQ=255)
                                                                                      acaaAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAc                                                             >  1:1407335/1‑150 (MQ=255)
                                                                                         aaaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCAGGCAACACCTATCGGTCGTGGTCAGCGTGGACAGATTGTGGCAacgc                                                          >  1:1044036/1‑150 (MQ=255)
                                                                                          aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                         >  1:190117/1‑150 (MQ=255)
                                                                                          aaTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgcc                                                         <  1:1200636/150‑1 (MQ=255)
                                                                                           aTCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAccgccg                                                        <  2:675271/150‑1 (MQ=255)
                                                                                           aTCCTCTTTGAGAACTTAACCCCGCCGCTCGCAATCTCTGGTCTCCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGAGCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCAcctccg                                                        <  2:292883/150‑1 (MQ=255)
                                                                                             ccTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaa                                                      <  2:1140623/150‑1 (MQ=255)
                                                                                              ctctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGaaa                                                     >  1:1599780/1‑150 (MQ=255)
                                                                                                ctTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGc                                                   >  2:740119/1‑150 (MQ=255)
                                                                                                  ttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTAGAGTACTGGATCTGGCATAACCTATAGGTCGTGGTCAGCTTTGTCTGATTTTGGCACCGCCGAAATCCg                                                 >  2:1100856/1‑150 (MQ=255)
                                                                                                  ttGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCg                                                 <  1:1486048/150‑1 (MQ=255)
                                                                                                   tGATAACTTAACCCCGCTGCACGCTAACACTCGTCTGCGTATGGAACAAGTTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                <  2:2096970/150‑1 (MQ=255)
                                                                                                   tGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCgg                                                <  1:770080/150‑1 (MQ=255)
                                                                                                    gagaACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAATATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGTTCAGCGTTGTCTGATTGTGGCACCGCCGAAAGCCGGt                                               <  1:1475652/150‑1 (MQ=255)
                                                                                                     agaACTTCACCCCGCTGCACGCACACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTa                                              <  1:92041/150‑1 (MQ=255)
                                                                                                       aaCTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTaaa                                            <  1:1728754/150‑1 (MQ=255)
                                                                                                         cTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAAc                                          >  1:1222463/1‑150 (MQ=255)
                                                                                                              ccccGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTTTCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgg                                     >  1:440046/1‑149 (MQ=255)
                                                                                                                 cGCTGCACGCAAACTCTCGTCTGCGTATGGACCGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGGGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgc                                  <  2:1436663/150‑1 (MQ=255)
                                                                                                                  gCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgct                                 <  2:649354/150‑1 (MQ=255)
                                                                                                                  gCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCGCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgct                                 <  1:1114502/150‑1 (MQ=255)
                                                                                                                   cTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctg                                <  2:479532/150‑1 (MQ=255)
                                                                                                                    tGCACGCAATCTCTCGTCTGCGTATGGAACGGGGTAACGGTTCTTCTGAAGATTTAAGTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTTGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCAtgctgctgc                               <  2:770295/150‑1 (MQ=255)
                                                                                                                     gCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCa                              <  2:816647/150‑1 (MQ=255)
                                                                                                                     gCACGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCa                              <  1:1631931/150‑1 (MQ=255)
                                                                                                                       aCGCAAACTCTCGTCTGCGTATGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGa                            >  1:1567305/1‑150 (MQ=255)
                                                                                                                          cAAACTCTCGTCTGCGAATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTAGCGTACTGGATCTGGCATCACCTATCGGTATTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCATAACa                         >  1:264852/1‑150 (MQ=255)
                                                                                                                             aCTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCATCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTTCAGAACAttt                      >  1:60154/1‑149 (MQ=255)
                                                                                                                              ctctCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGc                     >  2:1969532/1‑150 (MQ=255)
                                                                                                                              ctctCGGCTGCGTATGGAACGTGGTAACGGTGCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGc                     <  1:1843808/150‑1 (MQ=255)
                                                                                                                                  cGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCag                 >  2:106636/1‑150 (MQ=255)
                                                                                                                                      tGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCa             >  2:57221/1‑150 (MQ=255)
                                                                                                                                          tATGGAACGTGGTAACGGTTCTAATGAAGATTAAACTGCCCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGc         <  2:949592/150‑1 (MQ=255)
                                                                                                                                           aTGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt        >  1:1644013/1‑150 (MQ=255)
                                                                                                                                           aTGGAACATGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCt        <  1:1936229/150‑1 (MQ=255)
                                                                                                                                            tGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCtt       <  2:447719/150‑1 (MQ=255)
                                                                                                                                               aaCGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACa    <  1:496463/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGAGAACTTAACCCCGCTGCACGCAAACTCTCGTCTGCGTATGGAACGTGGTAACGGTTCTACTGAAGATTTAACTGCTCGCGTACTGGATCTGGCATCACCTATCGGTCGTGGTCAGCGTGGTCTGATTGTGGCACCGCCGAAAGCCGGTAAAACCATGCTGCTGCAGAACATTGCTCAGAGCATTGCTTACAAC  >  NC_000913/3966716‑3967010

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: