Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 572,058 A→T 8.6% E202V (GAG→GTG)  ybcM → DLP12 prophage; putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913572,0580AT 8.6% 1000.0 / 16.6 328E202V (GAG→GTG) ybcMDLP12 prophage; putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (153/130);  new base T (3/34);  total (157/171)
Fisher's exact test for biased strand distribution p-value = 3.30e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

TCATTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGT  >  NC_000913/571909‑572206
                                                                                                                                                     |                                                                                                                                                    
tCATTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGa                                                                                                                                                      >  2:2774355/1‑150 (MQ=255)
  aTTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACACGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                                                                                                                                    <  1:4682363/150‑1 (MQ=255)
  aTTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                                                                                                                                    >  1:4301192/1‑150 (MQ=255)
  aTTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                                                                                                                                    <  1:1732729/150‑1 (MQ=255)
  aTTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                                                                                                                                    <  2:3425733/150‑1 (MQ=255)
  aTTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAgg                                                                                                                                                    <  2:4624738/150‑1 (MQ=255)
   tttACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATGTCTAGATTTCTGTTTGTTCATTACATCGCCGTCTAGCAAGCGAggg                                                                                                                                                   >  1:4825180/1‑150 (MQ=255)
   tttACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAggg                                                                                                                                                   >  1:1819275/1‑150 (MQ=255)
        ttttATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTaa                                                                                                                                              >  2:4530206/1‑150 (MQ=255)
        ttttATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTaa                                                                                                                                              <  1:3628047/150‑1 (MQ=255)
        ttttATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTaa                                                                                                                                              <  1:4490174/150‑1 (MQ=255)
               atatAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCa                                                                                                                                       >  1:2586901/1‑150 (MQ=255)
                 ataGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCaaa                                                                                                                                     <  2:1983946/150‑1 (MQ=255)
                 ataGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCaaa                                                                                                                                     >  1:633470/1‑150 (MQ=255)
                 ataGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCaaa                                                                                                                                     <  2:621893/150‑1 (MQ=255)
                   aGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGt                                                                                                                                   >  2:1152352/1‑150 (MQ=255)
                     tAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTAGATTTCTGTTTCTAGATTAAATCGCCGTCTAGCAAGGGTGGGGGTAAGTTTCCAAAGtat                                                                                                                                 >  1:4579978/1‑150 (MQ=255)
                      aaaCAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGtata                                                                                                                                <  2:412453/150‑1 (MQ=255)
                      aaaCAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGtata                                                                                                                                <  1:742544/150‑1 (MQ=255)
                          aaTCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTgg                                                                                                                            >  1:4132983/1‑150 (MQ=255)
                           aTCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGa                                                                                                                           <  1:811097/150‑1 (MQ=255)
                           aTCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGa                                                                                                                           <  1:869424/150‑1 (MQ=255)
                            tCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGAc                                                                                                                          <  2:2027468/150‑1 (MQ=255)
                               ttGATAAGTCGTTGTTATAACCTGCTGATATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGAt                                                                                                                       <  2:5116360/150‑1 (MQ=255)
                                tGATAAGTCGTTGTTATAACCAGCAGCTATCCGAACCCGGCACAAAATGGAAAGCAAACAAGGTTTCTCGATATATCAACATTTCTGTTTCTACATTAAATAGCCGTCTATCAAGCGAGGGGGTAAGTTTACAAAGTATACTGGACGAtg                                                                                                                      >  1:4719314/1‑150 (MQ=255)
                                  aTAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGAtgtg                                                                                                                    <  1:1191185/150‑1 (MQ=255)
                                   tAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGa                                                                                                                   <  2:2162025/150‑1 (MQ=255)
                                      gTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGt                                                                                                                <  2:819455/150‑1 (MQ=255)
                                       tCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGtt                                                                                                               <  1:553335/150‑1 (MQ=255)
                                       tCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGtt                                                                                                               <  2:322963/150‑1 (MQ=255)
                                       tCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGtt                                                                                                               >  2:1246290/1‑150 (MQ=255)
                                         gttgttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAggg                                                                                                                                                   >  2:49477/1‑112 (MQ=255)
                                         gttgttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTaa                                                                                                             >  2:3189907/1‑150 (MQ=255)
                                         gttgttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTaa                                                                                                             <  1:537158/150‑1 (MQ=255)
                                          ttgttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTaaa                                                                                                            >  1:1450982/1‑150 (MQ=255)
                                            gttATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAaata                                                                                                          <  1:3007463/150‑1 (MQ=255)
                                                taACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGc                                                                                                      <  2:5015303/150‑1 (MQ=255)
                                                taACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGc                                                                                                      <  1:1036487/150‑1 (MQ=255)
                                                taACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGc                                                                                                      >  2:420904/1‑150 (MQ=255)
                                                   ccTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGtt                                                                                                   <  2:1397999/150‑1 (MQ=255)
                                                        tgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTg                                                                                              <  2:2641421/150‑1 (MQ=255)
                                                        tgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTg                                                                                              >  1:1275462/1‑150 (MQ=255)
                                                        tgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTg                                                                                              <  2:5428090/150‑1 (MQ=255)
                                                        tgctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTg                                                                                              >  1:597179/1‑150 (MQ=255)
                                                         gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTTCATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                             <  1:5061515/150‑1 (MQ=255)
                                                         gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGGGGTTAAATAATGCGTTGTCTGc                                                                                             <  2:2425643/150‑1 (MQ=255)
                                                         gctATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGc                                                                                             >  2:4926326/1‑150 (MQ=255)
                                                              ccGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATAc                                                                                        >  2:1712359/1‑150 (MQ=255)
                                                               cGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACa                                                                                       >  1:1592832/1‑150 (MQ=255)
                                                                gAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACaa                                                                                      >  1:2174124/1‑150 (MQ=255)
                                                                    ccGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAacga                                                                                  >  1:5159109/1‑150 (MQ=255)
                                                                      ggCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCTTCTAGCAAGCGAGGGGGTAAGTTTACAAAGTATACTGGACGATTTGAGGTTAAATAATGAGTTGTCTGCTATACAAacgacg                                                                                >  1:1922748/1‑150 (MQ=255)
                                                                      ggCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGAGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAacgacg                                                                                <  2:2202224/150‑1 (MQ=255)
                                                                       gCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACgg                                                                               >  1:2727027/1‑150 (MQ=255)
                                                                         acaAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTa                                                                             >  1:5009160/1‑150 (MQ=255)
                                                                          caAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTaa                                                                            >  2:4360582/1‑150 (MQ=255)
                                                                          aaaaaTGGACAGAAAAAAAGGTTGCTCGATATCTCAACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAATTTTCAAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTaa                                                                            <  2:5353047/149‑1 (MQ=255)
                                                                           aaaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTaaa                                                                           <  2:4719314/150‑1 (MQ=255)
                                                                           aaaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTaaa                                                                           >  1:2083775/1‑150 (MQ=255)
                                                                            aaaTGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTaaaa                                                                          >  2:2243506/1‑150 (MQ=255)
                                                                               tGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCt                                                                       >  2:3106692/1‑150 (MQ=255)
                                                                                ggACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCta                                                                      >  1:4170507/1‑150 (MQ=255)
                                                                                   cAGAAAAAACGGTTGCTCGATATCTCTACATTTGTTTTTCTAAATTACTTCGCCGTCTCGCACGCGCGGGGTTATTTTTCAAAAGTATAGGGGACTTTGTGTGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATaa                                                                   <  2:1856386/150‑1 (MQ=255)
                                                                                        aaCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCgag                                                              >  1:5102355/1‑150 (MQ=255)
                                                                                         aCAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCgaga                                                             >  2:4407546/1‑150 (MQ=255)
                                                                                              gTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATATTGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGcc                                                        <  2:1738083/150‑1 (MQ=255)
                                                                                               ttGCTCGCTATCTCTCCCTTTCTGTTTCTACTTTACCTCGCCGTCTCGCACGGGTGGGGGTATGTTTCCAATGTATAGGGGCGGTTGTGTGTTTAATTATTGGGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCa                                                       <  1:3570393/150‑1 (MQ=255)
                                                                                                   tCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAgaga                                                   >  2:3096925/1‑150 (MQ=255)
                                                                                                     gATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGaaa                                                 >  1:949697/1‑150 (MQ=255)
                                                                                                      ctatCTCTACATTTCTGTTTCTACATTACAGCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGaaaa                                                <  2:2456508/149‑1 (MQ=255)
                                                                                                      atatCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGCGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGaaaa                                                <  2:454408/150‑1 (MQ=255)
                                                                                                         tctcTACATTTCTTTTTCTACATTACATCGCCGTCTAGCAAGCGTGGGGGTAAGTTTACAAAGTATACTGGACGTTGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATgg                                             <  2:690604/150‑1 (MQ=255)
                                                                                                         tctcTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATgg                                             >  2:3655825/1‑150 (MQ=255)
                                                                                                           tctACATTTGTTTTTCTAAATTACATCGCGGTCTAGCAAGGGTGGGGGTATGTTTCAAATGTATACTGGACGTTGTGTGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGtt                                           <  1:5416497/150‑1 (MQ=255)
                                                                                                            ctACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGCGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTta                                          >  1:1674253/1‑150 (MQ=255)
                                                                                                             tACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTtat                                         >  1:3731621/1‑150 (MQ=255)
                                                                                                                 tttCTGTTTCTACATTACATCGCCGTCTAGCAAGCGTGGGGTAAAGTTTCCAAAGTATACTGGACGATGTGTGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAgt                                     <  1:3502979/150‑1 (MQ=255)
                                                                                                                 tttCTGTTTCTACATTACATCGCCGTCTAGCAAGCGTGGGGGTACGTTCCCAACGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAgt                                     <  2:3875226/150‑1 (MQ=255)
                                                                                                                 tttCTGTTTATAAATTACATCGCCGTCTAGAAAGGGAGGGGTTAAGTTTCCAAATTATACTGTACGATGTTTGGTTAAATAATGGGTTGTCTGATATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAgt                                     <  1:1748640/150‑1 (MQ=255)
                                                                                                                        ttCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTtc                              >  2:2604300/1‑150 (MQ=255)
                                                                                                                         tCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTCTACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTtct                             >  2:3576455/1‑150 (MQ=255)
                                                                                                                           tacattacatCGCCTTCTAGCAAGCGTGGGGGTAATTTTCCAAAGTATACTGGACTATGTGCGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCg                           <  1:2601746/150‑1 (MQ=255)
                                                                                                                              attacatCGCCGTCTAGCATGGGTGGGGGTATGTTTCAAATGTTTAGTGGACGATGTGAGTTTAAATTATGCGTTTTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                        <  2:4263660/150‑1 (MQ=255)
                                                                                                                              attacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACTATGTGAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                        <  2:145284/150‑1 (MQ=255)
                                                                                                                              attacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGttt                        <  1:542536/150‑1 (MQ=255)
                                                                                                                               ttacatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGtttt                       >  1:1607463/1‑150 (MQ=255)
                                                                                                                                tacatCGCCGTCTAGCAAGCGTGGGGGTAATTTTCCAAAGTATACTGGACGATTTGAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTa                      <  2:1086883/150‑1 (MQ=255)
                                                                                                                                 acatCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCAAACTTATACTGGACGATGTGAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTAc                     <  1:3932163/150‑1 (MQ=255)
                                                                                                                                 acatCGCCGTATAGCAAGCGCGGGTGTATGTTTCCAATGTATTCTTGACTTTGGTTTGTAAATTAATGCGTTGTCTGCTATAAAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTAc                     <  1:4916128/150‑1 (MQ=255)
                                                                                                                                 acagCGCGGTCTAGCTAGTGAGGGGTTCATTTACAAACGTTTACTGGAGTCTGTGCGGTTAATTACTGCGTTGTCTGCTATACCAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTAc                     <  1:2431613/150‑1 (MQ=255)
                                                                                                                                    tCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGa                  >  2:862327/1‑150 (MQ=255)
                                                                                                                                    acGCCGTCTAGCAAGCGTGGGGTTAAGTTTCCAAATTATACTGGACGATGTGAGTTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGa                  <  1:2924206/149‑1 (MQ=255)
                                                                                                                                        cGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGa              >  1:3188097/1‑150 (MQ=255)
                                                                                                                                            ttgCACGCGTGGGGGTAAGTTTCCACAGTATACTGGACGAGGTGTGGTTATATAATGCGTTGTCTGCTATACACACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTc          <  2:940526/148‑1 (MQ=255)
                                                                                                                                              gCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCAt        >  1:5151157/1‑150 (MQ=255)
                                                                                                                                              gCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCAt        >  1:3797388/1‑150 (MQ=255)
                                                                                                                                              gCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCAt        >  2:861258/1‑150 (MQ=255)
                                                                                                                                                  gCGTGGGGGTAAGTTTCCAAAGTATACTGGACGATGTTAGGTTAATTAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATc    <  1:3127179/150‑1 (MQ=255)
                                                                                                                                                  gCGAGGGGGTAAGTTTCAAAAGTATACTGGACGATGTGAGTTTAAATAATGCTTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATc    <  2:2902745/150‑1 (MQ=255)
                                                                                                                                                    gAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGt  >  1:2331261/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCATTTACTTTTATCATATAGTAAACAATCATTGATAAGTCGTTGTTATAACCTGCTGCTATCCGAACCCGGCACAAAATGGACAGCAAACAAGGTTGCTCGATATCTCTACATTTCTGTTTCTACATTACATCGCCGTCTAGCAAGCGAGGGGGTAAGTTTCCAAAGTATACTGGACGATGTGAGGTTAAATAATGCGTTGTCTGCTATACAAACGACGGTAAAACCTATAAGCGAGATTGCCAGAGAAAATGGTTATAAGTGTCCTTCTCGTTTTACTGAAAGATTTCATAATCGT  >  NC_000913/571909‑572206

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