Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,649,718 A→G E37G (GAA→GGA)  dicB → Qin prophage; cell division inhibition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,649,7180AG99.0% 367.6 / ‑4.6 104E37G (GAA→GGA) dicBQin prophage; cell division inhibition protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (53/50);  total (53/51)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 104 total aligned reads displayed.

CTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCA  >  NC_000913/1649575‑1649863
                                                                                                                                               |                                                                                                                                                 
cTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATACCATTAACAAGAGAAAACAAGAACggg                                                                                                                                              <  8:337315/149‑1 (MQ=255)
 tGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACgg                                                                                                                                               >  4:130150/1‑147 (MQ=255)
    aTCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGc                                                                                                                                           >  3:374661/1‑148 (MQ=255)
     tCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTa                                                                                                                                         <  2:219152/149‑1 (MQ=255)
       acaAAATTCATTGTTGAGGACGCGATAATGAATACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTTTCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTAtt                                                                                                                                       <  3:187759/149‑1 (MQ=255)
       acaAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTAt                                                                                                                                        >  4:117457/1‑148 (MQ=255)
         aaaaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATtaa                                                                                                                                     <  3:216409/149‑1 (MQ=255)
         aaaaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATta                                                                                                                                      >  2:47034/1‑148 (MQ=255)
         aaaaTTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATta                                                                                                                                      >  7:305187/1‑148 (MQ=255)
               cATTGTTGAGGACGCGATAATGAAAACTTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCAAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATtaaataaa                                                                                                                                <  5:137811/148‑1 (MQ=255)
                aTTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTaaataaaa                                                                                                                               <  1:47034/148‑1 (MQ=255)
                aTTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTaaataaaa                                                                                                                               >  3:67718/1‑148 (MQ=255)
                 ttgttgAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTaaataaaata                                                                                                                             <  3:117457/149‑1 (MQ=255)
                         gACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGt                                                                                                                     >  8:167440/1‑149 (MQ=255)
                          aCGCGGTAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGtt                                                                                                                    <  1:36907/149‑1 (MQ=255)
                           cgcgATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGttt                                                                                                                   <  4:374661/149‑1 (MQ=255)
                              gATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCaa                                                                                                                >  8:405442/1‑149 (MQ=255)
                               aTAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAAt                                                                                                               >  8:350734/1‑149 (MQ=255)
                                 aaTGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGAGATTAAATAAAATATGCATTGTTTCAATGc                                                                                                             >  4:287921/1‑149 (MQ=255)
                                  aTGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCt                                                                                                            <  1:225630/149‑1 (MQ=255)
                                   tGAAAACGTTATTACGAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCTGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTg                                                                                                           <  5:106567/149‑1 (MQ=255)
                                    gAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTTAATAAAATATGCATTGTTTCAATGCTgg                                                                                                          >  5:37054/1‑149 (MQ=255)
                                    gAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTgg                                                                                                          >  4:260254/1‑149 (MQ=255)
                                       aaCGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTc                                                                                                       <  8:73110/149‑1 (MQ=255)
                                          gTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGtt                                                                                                    <  8:18399/149‑1 (MQ=255)
                                          gTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGtt                                                                                                    <  3:212840/149‑1 (MQ=255)
                                           ttattaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGtt                                                                                                                    <  1:143073/132‑1 (MQ=255)
                                           ttattaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGtt                                                                                                                    >  2:143073/1‑132 (MQ=255)
                                            tattaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTa                                                                                                  >  8:214435/1‑149 (MQ=255)
                                             attaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGa                                                                                                                                             >  8:43036/1‑105 (MQ=255)
                                             attaCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGa                                                                                                                                             <  7:43036/105‑1 (MQ=255)
                                               taCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGt                                                                                               >  4:249763/1‑149 (MQ=255)
                                                aCCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTc                                                                                              >  5:295390/1‑149 (MQ=255)
                                                 ccAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATg                                                                                                                           >  8:56687/1‑119 (MQ=255)
                                                 ccAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATg                                                                                                                           <  7:56687/119‑1 (MQ=255)
                                                 ccAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCt                                                                                             >  4:51572/1‑149 (MQ=255)
                                                 ccAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCt                                                                                             >  2:83170/1‑149 (MQ=255)
                                                  cAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACTAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGATCGTTTACGTCt                                                                                             >  3:255627/1‑148 (MQ=255)
                                                  cAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCt                                                                                             >  1:394949/1‑148 (MQ=255)
                                                   aaaCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGa                                                                                           <  4:249988/149‑1 (MQ=255)
                                                    aaCTTTAATACGTATGAAGGTTGTTTTGAAATTTGTGTAACTATCAGTAACACAGTATTTACTTAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                                          <  3:194114/149‑1 (MQ=255)
                                                    aaCGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGAt                                                                                          <  8:305187/149‑1 (MQ=255)
                                                      cGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGc                                                                                        <  1:1663/149‑1 (MQ=255)
                                                       gTTAATACGTCTGAAGGATGTTTTGAAATTGGTCACACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGATATTAAAAAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGc                                                                                        >  4:321288/1‑148 (MQ=255)
                                                           aTACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCaaaa                                                                                   <  1:83170/149‑1 (MQ=255)
                                                             aCGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAg                                                                                  >  7:246533/1‑148 (MQ=255)
                                                                tCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACCAGGACGGGAGCTATTAAAAAAAATATGCATTGTTTCAATGCTCGCTCGTTTACGTCTGATGCCAAAAGGAtg                                                                              >  4:77902/1‑149 (MQ=255)
                                                                 cTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGt                                                                                               >  8:197963/1‑131 (MQ=255)
                                                                 cTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGt                                                                                               <  7:197963/131‑1 (MQ=255)
                                                                  tGAATGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCaaaa                                                                                   <  3:5576/142‑1 (MQ=255)
                                                                  tGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGTCTCGTTTACGTCTGATGCCAAAAGGAtttg                                                                            >  6:377163/1‑149 (MQ=255)
                                                                  tGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCaaaa                                                                                   >  4:5576/1‑142 (MQ=255)
                                                                    aaGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGc                                                                           <  8:208377/148‑1 (MQ=255)
                                                                       gttgttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACaa                                                                       >  4:156404/1‑149 (MQ=255)
                                                                          gttTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATg                                                                     >  4:360692/1‑148 (MQ=255)
                                                                           ttttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGaa                                                                   >  5:137821/1‑149 (MQ=255)
                                                                           ttttGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGaa                                                                   >  2:144518/1‑149 (MQ=255)
                                                                              tGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAAtt                                                                 <  8:246533/148‑1 (MQ=255)
                                                                               gAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                                                <  7:350734/148‑1 (MQ=255)
                                                                                aaaTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCa                                                               <  7:167440/148‑1 (MQ=255)
                                                                                aaaTTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCa                                                               <  7:405442/148‑1 (MQ=255)
                                                                                      gtgtCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTtg                                                        >  3:196698/1‑149 (MQ=255)
                                                                                        gtCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGTATGTGCACAATGAATTCAGCATTtgtg                                                      >  3:202743/1‑149 (MQ=255)
                                                                                         tCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATACGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGc                                                     <  4:280294/149‑1 (MQ=255)
                                                                                          cACTATCAGTTACCCAGTATTTACTGAATATGCGATTAACAAGAGAAAACTAGGACGGGAGCTATTAAATAAATTATGCATTGTTTCAATGCTGGCTCGTTTACGTGTGATGCCAAAATGATGTGCACTATGAATTCAGCATTTGTGc                                                     <  3:321288/148‑1 (MQ=255)
                                                                                          cACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCt                                                    >  3:70191/1‑149 (MQ=255)
                                                                                                 aGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCt                                              <  1:144518/148‑1 (MQ=255)
                                                                                                 aGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCt                                              <  3:104058/148‑1 (MQ=255)
                                                                                                 aGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCt                                              <  3:27480/148‑1 (MQ=255)
                                                                                                 aGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCt                                              <  4:196698/148‑1 (MQ=255)
                                                                                                  gTACCCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGATTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGa                                            >  8:131767/1‑149 (MQ=255)
                                                                                                  gTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTg                                             >  3:303540/1‑148 (MQ=255)
                                                                                                     aCCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAg                                         >  1:407782/1‑149 (MQ=255)
                                                                                                        cAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGttt                                      >  7:233343/1‑149 (MQ=255)
                                                                                                         aGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGtttt                                     <  3:130150/149‑1 (MQ=255)
                                                                                                            aTTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCtt                                                   <  1:111821/132‑1 (MQ=255)
                                                                                                            aTTTACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCtt                                                   >  2:111821/1‑132 (MQ=255)
                                                                                                             tttACTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCtt                                 <  7:255129/149‑1 (MQ=255)
                                                                                                                aCTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGt                               <  2:394949/148‑1 (MQ=255)
                                                                                                                 cTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGttt                                                                                                   >  6:365657/1‑79 (MQ=255)
                                                                                                                 cTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGttt                                                                                                   <  5:365657/79‑1 (MQ=255)
                                                                                                                 cTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGttt                             >  6:210342/1‑149 (MQ=255)
                                                                                                                 cTGAAGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGttt                             >  5:280610/1‑149 (MQ=255)
                                                                                                                    aaGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTcc                           <  6:137821/148‑1 (MQ=255)
                                                                                                                    aaGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCg                          <  6:37054/149‑1 (MQ=255)
                                                                                                                    aaGATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCg                          >  3:199024/1‑149 (MQ=255)
                                                                                                                      gATGCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGga                        <  3:260254/149‑1 (MQ=255)
                                                                                                                         gCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                                                <  1:93931/106‑1 (MQ=255)
                                                                                                                         gCCATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTc                                                                >  2:93931/1‑106 (MQ=255)
                                                                                                                           cATTAACAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGcc                   <  5:210342/149‑1 (MQ=255)
                                                                                                                               aaCAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGtt               >  1:151598/1‑149 (MQ=255)
                                                                                                                                 cAAGAGAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGa             <  4:70191/149‑1 (MQ=255)
                                                                                                                                   agagAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGata           <  4:67718/149‑1 (MQ=255)
                                                                                                                                   agagAAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGata           <  6:280610/149‑1 (MQ=255)
                                                                                                                                    gagaAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGga                        >  1:289303/1‑135 (MQ=255)
                                                                                                                                    gagaAAACAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGga                        <  2:289303/135‑1 (MQ=255)
                                                                                                                                         aaCAAGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCCGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGcag     >  4:78307/1‑149 (MQ=255)
                                                                                                                                           cAAGGACGGGAGCTATTAAATAAAATATGGATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGcagt    >  7:270605/1‑148 (MQ=255)
                                                                                                                                            aaGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGcagtca  <  2:407782/149‑1 (MQ=255)
                                                                                                                                            aaGGACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGcagt    >  4:114707/1‑147 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
CTGCATCACAAAATTCATTGTTGAGGACGCGATAATGAAAACGTTATTACCAAACGTTAATACGTCTGAAGGTTGTTTTGAAATTGGTGTCACTATCAGTAACCCAGTATTTACTGAAGATGCCATTAACAAGAGAAAACAAGAACGGGAGCTATTAAATAAAATATGCATTGTTTCAATGCTGGCTCGTTTACGTCTGATGCCAAAAGGATGTGCACAATGAATTCAGCATTTGTGCTTGTTCTGACAGTTTTTCTTGTTTCCGGAGAGCCAGTTGATATTGCAGTCA  >  NC_000913/1649575‑1649863

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