Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,227,053 G→C 6.8% intergenic (+85/+18) ycgJ → / ← pliG PF05666 family protein YcgJ/inhibitor of g‑type lysozyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,227,0530GC 6.8% 350.6 / 12.3 103intergenic (+85/+18)ycgJ/pliGPF05666 family protein YcgJ/inhibitor of g‑type lysozyme
Reads supporting (aligned to +/- strand):  ref base G (49/47);  new base C (3/4);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 106 total aligned reads displayed.

GTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGG  >  NC_000913/1226907‑1227202
                                                                                                                                                  |                                                                                                                                                     
gTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGCGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTcc                                                                                                                                                    >  1:849039/1‑150 (MQ=255)
gTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTcc                                                                                                                                                    >  2:1031951/1‑150 (MQ=255)
 tAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTccc                                                                                                                                                   >  1:158773/1‑150 (MQ=255)
 tAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTccc                                                                                                                                                   >  1:1902104/1‑150 (MQ=255)
       gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTTAGTGATGTCTTGTTTACATTATGTTCTATTCTTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCGTcat                                                                                                                                             <  2:77108/150‑3 (MQ=255)
       gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCgg                                                                                                                                             <  1:1985274/150‑1 (MQ=255)
       gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCGTcat                                                                                                                                             <  1:498748/150‑3 (MQ=255)
       gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCGTcat                                                                                                                                             <  1:1802744/150‑3 (MQ=255)
        aGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACTTTTATGACGGGGAGAGTCCCCGTCGGt                                                                                                                                            <  1:131724/150‑1 (MQ=255)
        aGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCGTcata                                                                                                                                            <  1:1675351/150‑4 (MQ=255)
         gTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGCATTAACGTTTATGACTGGGAGAGTCCCCGTCGGTTg                                                                                                                                          <  2:1951195/150‑1 (MQ=255)
           ggAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGa                                                                                                                                         >  1:1516002/1‑150 (MQ=255)
           ggAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGa                                                                                                                                         >  2:1315093/1‑150 (MQ=255)
            gAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGAc                                                                                                                                        >  2:695687/1‑150 (MQ=255)
             aaaaaTCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACa                                                                                                                                       <  2:748723/150‑1 (MQ=255)
                   caatcaaACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACt                                                                                                                                 >  1:1748707/1‑150 (MQ=255)
                    aaacaaACAACAAAAAAGAAATAATTTATCTGTCTTTAGTTATCTCTTGTATACATTAAATTCTATAATACATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGAGAGTCCCCGTCGGTTGACATTTACTTa                                                                                                                               >  2:165904/4‑150 (MQ=255)
                        aaaCCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAAtt                                                                                                                            <  2:1798095/150‑1 (MQ=255)
                        aaaCCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAAtt                                                                                                                            <  1:384821/150‑1 (MQ=255)
                         aaCCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAAttt                                                                                                                           <  1:1384982/150‑1 (MQ=255)
                         aaCCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAAttt                                                                                                                           >  1:176368/1‑150 (MQ=255)
                            caccacAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGaa                                                                                                                        <  2:953593/150‑1 (MQ=255)
                             accacAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAat                                                                                                                       <  2:1904375/150‑1 (MQ=255)
                              ccacAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAata                                                                                                                      <  1:1276534/150‑1 (MQ=255)
                              ccacAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAata                                                                                                                      <  1:1962907/150‑1 (MQ=255)
                                 caAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCg                                                                                                                   >  2:1353379/1‑150 (MQ=255)
                                  aaaGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGa                                                                                                                  <  2:1902104/150‑1 (MQ=255)
                                    aGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACg                                                                                                                <  1:575394/150‑1 (MQ=255)
                                    aGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACg                                                                                                                <  1:17893/150‑1 (MQ=255)
                                    aGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACg                                                                                                                <  1:1216292/150‑1 (MQ=255)
                                       tGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTta                                                                                                             >  1:1813413/1‑150 (MQ=255)
                                       tGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTta                                                                                                             >  1:374210/1‑150 (MQ=255)
                                           atttatGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAtt                                                                                                         >  1:1189696/1‑150 (MQ=255)
                                              tatGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAttttt                                                                                                      <  1:695687/150‑1 (MQ=255)
                                               atGTGTCGTGAGTGTTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGTGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAtttttt                                                                                                     <  1:996007/150‑1 (MQ=255)
                                               atGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAtttttt                                                                                                     <  1:1315093/150‑1 (MQ=255)
                                               atGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAtttttt                                                                                                     <  1:1666840/150‑1 (MQ=255)
                                                 gtgtCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAg                                                                                                   <  1:1874657/150‑1 (MQ=255)
                                                  tgtCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgt                                                                                                  >  2:95278/1‑150 (MQ=255)
                                                   gtCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGTGCATTAACGTTTATGACGGGGAGAGTCACCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgtt                                                                                                 >  1:1895711/1‑150 (MQ=255)
                                                     cGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgtttt                                                                                               >  2:1841676/1‑150 (MQ=255)
                                                        gagtgaTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgttttgtt                                                                                            >  2:183663/1‑150 (MQ=255)
                                                          gtgaTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgttttgtttt                                                                                          <  1:1659602/150‑1 (MQ=255)
                                                             aTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTAc                                                                                       <  2:541201/150‑1 (MQ=255)
                                                                   tGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGcgt                                                                                 >  2:763390/1‑150 (MQ=255)
                                                                         cATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTc                                                                           >  1:1332011/1‑150 (MQ=255)
                                                                          aTTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCg                                                                          >  1:1455116/1‑150 (MQ=255)
                                                                            tATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCggg                                                                        <  2:374210/150‑1 (MQ=255)
                                                                               gTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGttt                                                                     >  1:1394817/1‑150 (MQ=255)
                                                                                ttCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGTGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTg                                                                    <  2:1338502/150‑1 (MQ=255)
                                                                                ttCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTg                                                                    >  1:1407565/1‑150 (MQ=255)
                                                                                 tCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGc                                                                   <  2:1748707/150‑1 (MQ=255)
                                                                                 tCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGc                                                                   >  1:845935/1‑150 (MQ=255)
                                                                                  cTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCa                                                                  <  2:1813413/150‑1 (MQ=255)
                                                                                  cTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCa                                                                  <  1:1353379/150‑1 (MQ=255)
                                                                                    aTTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAac                                                                >  1:1671717/1‑150 (MQ=255)
                                                                                          cATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGc                                                          <  2:808673/150‑1 (MQ=255)
                                                                                          cATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGc                                                          <  2:1455116/150‑1 (MQ=255)
                                                                                          cATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGc                                                          <  2:138576/150‑1 (MQ=255)
                                                                                               aCCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTc                                                     <  1:1841676/150‑1 (MQ=255)
                                                                                               aCCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCa                                                    <  1:165904/150‑1 (MQ=255)
                                                                                                 ctctTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCat                                                   <  1:866680/150‑1 (MQ=255)
                                                                                                     tGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATAttt                                               >  2:1058012/1‑150 (MQ=255)
                                                                                                          gACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAg                                          <  1:183663/150‑1 (MQ=255)
                                                                                                           aCAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGa                                         <  2:1407565/150‑1 (MQ=255)
                                                                                                              gagagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGc                                      >  1:656036/1‑150 (MQ=255)
                                                                                                              gagagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGc                                      >  2:1495324/1‑150 (MQ=255)
                                                                                                              gagagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGc                                      >  1:1102537/1‑150 (MQ=255)
                                                                                                               agagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCa                                     >  1:1940598/1‑150 (MQ=255)
                                                                                                               agagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCa                                     <  2:158773/150‑1 (MQ=255)
                                                                                                                   gagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGga                                 >  2:398423/1‑150 (MQ=255)
                                                                                                                     gTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagt                               <  2:1671717/150‑1 (MQ=255)
                                                                                                                     gTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagt                               <  1:1783171/150‑1 (MQ=255)
                                                                                                                     gTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagt                               >  2:1144699/1‑150 (MQ=255)
                                                                                                                      tGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtg                              <  2:158138/150‑1 (MQ=255)
                                                                                                                       ggggCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTCAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtga                             >  2:542247/1‑150 (MQ=255)
                                                                                                                          gCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTa                          >  1:393156/1‑150 (MQ=255)
                                                                                                                           cATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGTGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTAt                         >  2:1012363/1‑150 (MQ=255)
                                                                                                                              tAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGt                      <  2:1332011/150‑1 (MQ=255)
                                                                                                                               aaCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTc                     >  2:1253890/1‑150 (MQ=255)
                                                                                                                                 gtcaaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCg                   >  1:1661319/8‑150 (MQ=255)
                                                                                                                                 cGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCg                   >  1:288406/1‑150 (MQ=255)
                                                                                                                                  tcaaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGt                  >  2:175635/7‑150 (MQ=255)
                                                                                                                                   tttATGACGGGGAGAGTCCCCGTCTGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGGAACACTCGCAGCTCATATTTGCGAGAAGCAGTGAGTGAGTATTGTCCGTg                 <  2:1516002/150‑1 (MQ=255)
                                                                                                                                   tttATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTg                 <  2:52290/150‑1 (MQ=255)
                                                                                                                                     tATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc               <  1:640302/150‑1 (MQ=255)
                                                                                                                                     tATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc               <  1:1624522/150‑1 (MQ=255)
                                                                                                                                      aTGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACt              >  2:1815145/1‑150 (MQ=255)
                                                                                                                                      aTGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACt              >  2:1434453/1‑150 (MQ=255)
                                                                                                                                       tGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgag                                <  1:414262/131‑1 (MQ=255)
                                                                                                                                       tGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgag                                >  2:414262/1‑131 (MQ=255)
                                                                                                                                        gACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACTCTCGCAGCTCATATTTGCCAGAAGCAGGTAGTGAGTATTGTCCGTGACTAt            >  2:582759/1‑150 (MQ=255)
                                                                                                                                           tgggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGa         >  2:1418500/2‑150 (MQ=255)
                                                                                                                                           ggggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGa         >  1:743074/1‑150 (MQ=255)
                                                                                                                                            gggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAg        >  1:1695824/1‑150 (MQ=255)
                                                                                                                                            gggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAg        >  2:607919/1‑150 (MQ=255)
                                                                                                                                             ggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGt       >  1:1016581/1‑150 (MQ=255)
                                                                                                                                                 agTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTg   <  1:1281756/150‑1 (MQ=255)
                                                                                                                                                 agTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTg   <  1:1012363/150‑1 (MQ=255)
                                                                                                                                                  gTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTgg  >  2:669729/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCTGG  >  NC_000913/1226907‑1227202

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