Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,227,052 A→T 10.3% intergenic (+84/+19) 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,0520AT10.3% 272.5 / 18.6 97intergenic (+84/+19)ycgJ/pliGPF05666 family protein YcgJ/inhibitor of g‑type lysozyme
Reads supporting (aligned to +/- strand):  ref base A (38/49);  new base T (7/3);  total (45/52)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.11e-01

GATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCT  >  NC_000913/1226903‑1227200
                                                                                                                                                     |                                                                                                                                                    
caTGTTAAGAGTAGTGGATAAAACAATCAAACCAAAACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGgaga                                                                                                                                                      >  2:1093926/2‑150 (MQ=255)
  tGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGt                                                                                                                                                    <  2:1414704/150‑1 (MQ=255)
    gTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTcc                                                                                                                                                  >  1:41029/1‑150 (MQ=255)
     tAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTccc                                                                                                                                                 <  1:9926/150‑1 (MQ=255)
     tAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTccc                                                                                                                                                 >  1:46642/1‑150 (MQ=255)
       aGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCg                                                                                                                                               <  2:1034618/150‑1 (MQ=255)
       aGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCg                                                                                                                                               >  2:423894/1‑150 (MQ=255)
          ggAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCg                                                                                                                                            >  2:592453/1‑150 (MQ=255)
          ggAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCGTca                                                                                                                                            <  2:1231789/150‑2 (MQ=255)
           gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCgg                                                                                                                                           <  1:286922/150‑1 (MQ=255)
           gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGACTCTCCCCGTcat                                                                                                                                           <  1:423894/150‑3 (MQ=255)
           gAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGGGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCgg                                                                                                                                           <  2:1005495/150‑1 (MQ=255)
               ggAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGa                                                                                                                                       >  1:1153730/1‑150 (MQ=255)
               ggAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGa                                                                                                                                       <  1:777771/150‑1 (MQ=255)
                gAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGAc                                                                                                                                      >  2:1525843/1‑150 (MQ=255)
                gAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGAc                                                                                                                                      <  1:592453/150‑1 (MQ=255)
                    aatcaatcaatcCAGCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTGTATTCGTCATTAACCTCTTGAGGGACAGAGATGAGTGTGGCTTTAAGGTTTATGACGGGGAGAGTCCCCGTCGGTTGAGAttt                                                                                                                                  <  2:480249/150‑1 (MQ=255)
                    aatcaatcAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACAttt                                                                                                                                  >  1:132958/1‑150 (MQ=255)
                    aatcaatcAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACAttt                                                                                                                                  >  1:1492889/1‑150 (MQ=255)
                      tcaatcaaACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTAc                                                                                                                                <  2:1492889/150‑1 (MQ=255)
                       caatcaaACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACt                                                                                                                               <  2:829881/150‑1 (MQ=255)
                         atcaaACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTa                                                                                                                             <  2:1505274/150‑1 (MQ=255)
                          tcaaACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTaa                                                                                                                            >  2:1228834/1‑150 (MQ=255)
                            aaaCCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAAtt                                                                                                                          <  2:775474/150‑1 (MQ=255)
                                 accacAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGAGAAAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCACGTAGGTTGACATTTACTTAATTTGAat                                                                                                                     >  2:267559/1‑150 (MQ=255)
                                  ccacAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAata                                                                                                                    <  2:237704/150‑1 (MQ=255)
                                    acaAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATc                                                                                                                  <  2:1173327/150‑1 (MQ=255)
                                     caAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCg                                                                                                                 >  2:87996/1‑150 (MQ=255)
                                      aaaGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGa                                                                                                                <  1:345156/150‑1 (MQ=255)
                                        aGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGGGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACg                                                                                                              <  2:1219575/150‑1 (MQ=255)
                                           tGTTAAATATCTGTCGTGAGAGAAGTCTTGTTTAAATTATGTTCTATTCGTAATTAACCTCTTGAGCCACATAGAGGAATAGGACATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTta                                                                                                           >  2:690488/1‑150 (MQ=255)
                                               atttatGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAtt                                                                                                       >  2:79289/1‑150 (MQ=255)
                                                  tatGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATAttttt                                                                                                    >  1:1137963/1‑150 (MQ=255)
                                                    tgtgTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTa                                                                                                  <  2:552879/150‑1 (MQ=255)
                                                      tgtCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgt                                                                                                <  1:1356435/150‑1 (MQ=255)
                                                       gtCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAat                                                                                                                     <  2:15319/128‑1 (MQ=255)
                                                       gtCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAat                                                                                                                     >  1:15319/1‑128 (MQ=255)
                                                         catgagtgaTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgtttt                                                                                             >  1:1447881/3‑150 (MQ=255)
                                                         cGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgtttt                                                                                             >  1:811739/1‑150 (MQ=255)
                                                         cGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgtttt                                                                                             >  2:473234/1‑150 (MQ=255)
                                                            gagtgaTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAgttttgtt                                                                                          >  2:732129/1‑150 (MQ=255)
                                                                gaTGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTa                                                                                      <  1:935682/150‑1 (MQ=255)
                                                                  tGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACg                                                                                    <  2:173796/150‑1 (MQ=255)
                                                                   gTCTTGTTTACATAATTTTCTATTCGTCATTAACCTCTTGAGCGACAGAGATGAGTGGGGCATTAACTTTTATGACGGTGATAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACgg                                                                                   <  2:542568/150‑1 (MQ=255)
                                                                        gTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTtat                                                                                                          <  2:1358312/122‑1 (MQ=255)
                                                                        gTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTtat                                                                                                          >  1:1358312/1‑122 (MQ=255)
                                                                         tttACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGcgtcg                                                                             >  2:773406/1‑150 (MQ=255)
                                                                           tACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGtt                                                                           >  2:285599/1‑150 (MQ=255)
                                                                           tACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGtt                                                                           >  1:635415/1‑150 (MQ=255)
                                                                            aCATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGttt                                                                          <  1:1273820/150‑1 (MQ=255)
                                                                            aCATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGttt                                                                          <  1:1268647/150‑1 (MQ=255)
                                                                              aTTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCg                                                                        <  2:1120382/150‑1 (MQ=255)
                                                                                tATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCggg                                                                      <  1:160830/150‑1 (MQ=255)
                                                                                tATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGCGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCggg                                                                      <  1:149390/150‑1 (MQ=255)
                                                                                   gTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGttt                                                                   >  1:1085145/1‑150 (MQ=255)
                                                                                      cTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCa                                                                >  2:1264553/1‑150 (MQ=255)
                                                                                        aTTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAac                                                              <  1:473234/150‑1 (MQ=255)
                                                                                           cGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGAGTCGTTTCGGGTTTGCAACACt                                                           >  2:59829/1‑150 (MQ=255)
                                                                                              cATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGc                                                        <  1:577976/150‑1 (MQ=255)
                                                                                              cATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGc                                                        <  1:274730/150‑1 (MQ=255)
                                                                                                 tAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGc                                                     <  1:1133886/150‑1 (MQ=255)
                                                                                                  aaCCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCt                                                    >  2:428778/1‑150 (MQ=255)
                                                                                                    ccTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCa                                                  <  1:796058/150‑1 (MQ=255)
                                                                                                             cGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCa                                         >  1:552244/1‑150 (MQ=255)
                                                                                                                cAGAGAGGAGTGGTGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGaa                                      <  2:1430540/150‑1 (MQ=255)
                                                                                                                cAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGaa                                      <  1:973219/150‑1 (MQ=255)
                                                                                                                  gagagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGc                                    >  1:1203469/1‑150 (MQ=255)
                                                                                                                   agagGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCTACACTCGCAGCTCATATTTGCCAGAAGCa                                   >  2:1219747/1‑150 (MQ=255)
                                                                                                                    gaggagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACggg                                                                                  <  2:237993/102‑1 (MQ=255)
                                                                                                                    gaggagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACggg                                                                                  >  1:237993/1‑102 (MQ=255)
                                                                                                                      ggagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAggg                                <  2:244343/150‑1 (MQ=255)
                                                                                                                      ggagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAggg                                <  1:854289/150‑1 (MQ=255)
                                                                                                                       gagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGga                               <  1:79289/150‑1 (MQ=255)
                                                                                                                       gagTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGga                               <  2:811739/150‑1 (MQ=255)
                                                                                                                        agTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgag                              <  1:227904/150‑1 (MQ=255)
                                                                                                                         gTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagt                             >  2:428628/1‑150 (MQ=255)
                                                                                                                            gggCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtgag                          <  1:318558/150‑1 (MQ=255)
                                                                                                                             ggCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtgagt                         <  1:773406/150‑1 (MQ=255)
                                                                                                                             ggCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGgagtgagt                         >  2:1220227/1‑150 (MQ=255)
                                                                                                                               cATTAAAGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTAATTAATTTTAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTAt                       <  1:267559/150‑1 (MQ=255)
                                                                                                                                 ttAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTg                     <  2:41029/150‑1 (MQ=255)
                                                                                                                                   aaCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTc                   >  2:1459257/1‑150 (MQ=255)
                                                                                                                                    aCGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTcc                  <  1:381347/150‑1 (MQ=255)
                                                                                                                                     cGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTTTATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGATGCAGGGAGTGAGTATTGTCCg                 >  2:1395009/1‑150 (MQ=255)
                                                                                                                                       caaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGcc                                          >  2:1174531/6‑123 (MQ=255)
                                                                                                                                       caaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGcc                                          <  1:1174531/118‑1 (MQ=255)
                                                                                                                                       caaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTg               >  2:724453/6‑150 (MQ=255)
                                                                                                                                       caaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTg               >  2:729827/6‑150 (MQ=255)
                                                                                                                                        aaccgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGa              >  2:676694/5‑150 (MQ=255)
                                                                                                                                         tttGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc             >  2:944197/3‑150 (MQ=255)
                                                                                                                                         accgACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGAc             >  1:829283/4‑150 (MQ=255)
                                                                                                                                          aTGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACt            <  1:1228834/150‑1 (MQ=255)
                                                                                                                                           tGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTa           <  2:1447881/150‑1 (MQ=255)
                                                                                                                                            gACGGGGACTCTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGAATTGTCCGTGACTAt          >  1:1260286/1‑150 (MQ=255)
                                                                                                                                              cGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCg        >  2:562525/1‑150 (MQ=255)
                                                                                                                                              cGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCg        >  2:138014/1‑150 (MQ=255)
                                                                                                                                              cGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCg        <  2:132958/150‑1 (MQ=255)
                                                                                                                                                gggAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAg      >  2:445284/1‑150 (MQ=255)
                                                                                                                                                    gagTGCCCGTCGGTTGAAATTTACATAATTTTATTATCGACCGGATATTATTTAGTTTTGTTTTTAAGGGCGGGGGTTCGGGTTTGCCACACTCGAACGTCATATTTACCAGTAGCAGGGAGTGAGTATTGTCCGAGACTAGCGCGATCt  >  2:880916/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATGGTAAGCGGAGTGGAAAAATCAATCAAACCACCACAAAGATGTTATTTATGTGTCGTGAGTGATGTCTTGTTTACATTATGTTCTATTCGTCATTAACCTCTTGAGCGACAGAGAGGAGTGGGGCATTAACGTTTATGACGGGGAGAGTCCCCGTCGGTTGACATTTACTTAATTTGAATATCGACGTTATATTTTTTAGTTTTGTTTTTACGGGCGTCGTTTCGGGTTTGCAACACTCGCAGCTCATATTTGCCAGAAGCAGGGAGTGAGTATTGTCCGTGACTATCGAGTTCT  >  NC_000913/1226903‑1227200

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: