Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 843,749 C→A D578Y (GAT→TAT)  ybiO ← mechanosensitive channel protein, intermediate conductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913843,7490CA98.9% 334.7 / NA 93D578Y (GAT→TAT) ybiOmechanosensitive channel protein, intermediate conductance
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (52/40);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

CCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGACAACA  >  NC_000913/843602‑843895
                                                                                                                                                   |                                                                                                                                                  
ccGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATcc                                                                                                                                                 <  1:1172370/151‑1 (MQ=255)
ccGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATcc                                                                                                                                                 <  2:1095881/151‑1 (MQ=255)
ccGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATcc                                                                                                                                                 <  1:332844/151‑1 (MQ=255)
  gAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAg                                                                                                                                               >  2:108294/1‑151 (MQ=255)
  caaCCAATGCCGCGGACAATGTTTGCAAAGGTGGTTATCGAAGACCACTGGAGGATGTGATACGCCCCGGTATCCTGTCGTACGCCCCCGGAGCTACGCGAGATCCGTTACACTGTGCCGGTCATCGGCCCGATAGTCACAAAAGATCCat                                                                                                                                               <  2:998458/150‑2 (MQ=255)
   aaCCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgt                                                                                                                                              >  2:766834/1‑151 (MQ=255)
     ccAATGCCGCTGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACTCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                                                                            <  2:818238/151‑1 (MQ=255)
     ccAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAgtgt                                                                                                                                            <  1:526147/151‑1 (MQ=255)
      cAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGtt                                                                                                                                           >  2:1104379/1‑151 (MQ=255)
            cgcgGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGcc                                                                                                                                     <  1:647573/151‑1 (MQ=255)
               ggACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGtt                                                                                                                                  >  1:786017/1‑151 (MQ=255)
                    aaGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCaa                                                                                                                             >  2:162769/1‑151 (MQ=255)
                       ttAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACt                                                                                                                          <  2:135733/151‑1 (MQ=255)
                           cAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaat                                                                                                                      >  2:258756/1‑151 (MQ=255)
                               ggtggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                                                                                  >  2:409912/1‑151 (MQ=255)
                               ggtggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                                                                                  >  1:748597/1‑151 (MQ=255)
                               ggtggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaat                                                                                                                  >  1:1108531/1‑151 (MQ=255)
                                 tggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGaataaata                                                                                                                 >  2:1013387/1‑150 (MQ=255)
                                 tggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                                                                                                                <  2:450294/151‑1 (MQ=255)
                                 tggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGCGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                                                                                                                <  2:283239/151‑1 (MQ=255)
                                 tggtTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGGCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGTCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAc                                                                                                                <  2:653821/151‑1 (MQ=255)
                                    ttATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATAccc                                                                                                              <  1:281566/150‑1 (MQ=255)
                                             aCCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTGCGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTAAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataat                                                                                                    <  2:532510/151‑1 (MQ=255)
                                             aCCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataat                                                                                                    <  1:650180/151‑1 (MQ=255)
                                              ccACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                                                                   <  2:939262/151‑1 (MQ=255)
                                              ccACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGataata                                                                                                   <  2:231803/151‑1 (MQ=255)
                                               cACGGAATGATGTGATACGCCCCGGTATCGTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAGCGGCCGGATAGTCCCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAatat                                                                                                  <  2:128866/151‑1 (MQ=255)
                                               cACGGAATGATGTGATACGCCACGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAatat                                                                                                  <  2:955603/151‑1 (MQ=255)
                                                ccGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATc                                                                                                 <  2:384908/150‑1 (MQ=255)
                                                aCGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATc                                                                                                 <  1:301991/151‑1 (MQ=255)
                                                 cGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCt                                                                                                <  1:140037/151‑1 (MQ=255)
                                                   gAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCttt                                                                                              >  1:952456/1‑151 (MQ=255)
                                                   gAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCttt                                                                                              >  1:447835/1‑151 (MQ=255)
                                                   gAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCttt                                                                                              >  1:524827/1‑151 (MQ=255)
                                                   gAATGATGTGATACGCCCCGGTATCCTGTCGTACCCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGCTAATATCttt                                                                                              >  2:338428/1‑151 (MQ=255)
                                                     atgatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCa                                                                                            >  1:966012/1‑151 (MQ=255)
                                                       gatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                                                                           >  2:170474/1‑150 (MQ=255)
                                                       gatgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAc                                                                                           >  1:821273/1‑150 (MQ=255)
                                                        atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGGGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCAcc                                                                                          <  2:561039/150‑1 (MQ=255)
                                                        atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCa                                                                                         >  1:568587/1‑151 (MQ=255)
                                                        atgTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCa                                                                                         >  1:1138706/1‑151 (MQ=255)
                                                          gtgATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAg                                                                                        >  1:1072268/1‑150 (MQ=255)
                                                                  ccccGGTATCCTGTCGTACGCTCGCTGTGCGAATCGACATCCGTTCCGCTGTGCCGTTCAACGGCACGATAGTCGCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATGTCTTTCACCAGCGTTTGCGa                                                                               <  2:1193039/151‑1 (MQ=255)
                                                                  ccccGGCATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGACGGTCAACGGCCCGATTGTCCCCAAATCTTCAGTGTTTATGCCGTTTTCAAACTGAATAAATACCCCGGTGATATTATCTTTCACCAGCGTTCGCGa                                                                               >  2:1082234/1‑151 (MQ=255)
                                                                    ccGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGaa                                                                              >  1:992453/1‑150 (MQ=255)
                                                                      ggTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCa                                                                           >  1:509181/1‑151 (MQ=255)
                                                                       gTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCaa                                                                          >  2:929716/1‑151 (MQ=255)
                                                                       gTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCaa                                                                          >  1:12370/1‑151 (MQ=255)
                                                                         aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCaaa                                                                         >  2:371755/1‑150 (MQ=255)
                                                                         aTCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCaaa                                                                         >  1:577958/1‑150 (MQ=255)
                                                                              gTCGTACGCCCACGGCGCGCATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGACGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                                                                   >  2:1012787/1‑151 (MQ=255)
                                                                              gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                                                                   >  1:604736/1‑151 (MQ=255)
                                                                              gTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                                                                   >  1:1189664/1‑151 (MQ=255)
                                                                               tCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGAt                                                                   <  1:333784/150‑1 (MQ=255)
                                                                                  tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc                                                               <  1:849383/151‑1 (MQ=255)
                                                                                  tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc                                                               <  1:428844/151‑1 (MQ=255)
                                                                                  tACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGcc                                                               <  1:603480/151‑1 (MQ=255)
                                                                                      cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGa                                                            >  2:768953/1‑150 (MQ=255)
                                                                                      cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGa                                                            >  1:859145/1‑150 (MQ=255)
                                                                                      cccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAc                                                           <  1:331483/151‑1 (MQ=255)
                                                                                       ccACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGAcc                                                          <  2:545828/151‑1 (MQ=255)
                                                                                        cccGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTTACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt                                                         <  1:519267/149‑1 (MQ=255)
                                                                                        cACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCt                                                         <  1:979736/151‑1 (MQ=255)
                                                                                          cGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCCCCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTaa                                                       <  2:579284/151‑1 (MQ=255)
                                                                                             aGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATg                                                     <  1:970042/150‑1 (MQ=255)
                                                                                               cGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGcc                                                   <  1:454653/150‑1 (MQ=255)
                                                                                                     gACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCa                                             >  1:611889/1‑150 (MQ=255)
                                                                                                      aCATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCAc                                            <  1:769292/150‑1 (MQ=255)
                                                                                                         tCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTg                                         >  1:197760/1‑150 (MQ=255)
                                                                                                          ccGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGc                                        >  1:752151/1‑150 (MQ=255)
                                                                                                            gTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAg                                     >  1:817267/1‑151 (MQ=255)
                                                                                                               cctcTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATCCCCCGGTGATAATATCTTTCCCCACCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCaa                                  <  2:305061/148‑1 (MQ=255)
                                                                                                               ccACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCaa                                  >  2:921867/1‑151 (MQ=255)
                                                                                                               ccACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCaa                                  >  2:218345/1‑151 (MQ=255)
                                                                                                               ccACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCaa                                  >  2:680524/1‑151 (MQ=255)
                                                                                                               ccACTGTGCCGGTCAACGGCCCGATAGACACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGCGATAATATCTTTCATCAGCGCTTGCGAAACAAACGAGACAGCAACACCTAATGCCCCGGAACATGCCAGCa                                   >  2:763272/1‑150 (MQ=255)
                                                                                                                cACTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAAt                                 <  2:453975/151‑1 (MQ=255)
                                                                                                                  cTGTGCCGGTCAACTGCCCGATAGTCACCAAATATCCGGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGCCCTAATGCCCCGGCACATGCCAGCAATgg                               >  2:1064988/1‑151 (MQ=255)
                                                                                                                  cTGTGCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTATTCATGCCGCTTTTAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCTCTTGCGAACCAAACGAGTTAGCCAGACCTAATGCCCCGGCACCTGCCATCAATgg                               >  2:100691/1‑151 (MQ=255)
                                                                                                                    gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGgcg                             >  2:985705/1‑151 (MQ=255)
                                                                                                                    gtgCCGGTCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGgcg                             >  1:17316/1‑151 (MQ=255)
                                                                                                                           tCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCNatat                       <  2:394802/150‑1 (MQ=255)
                                                                                                                           tCAACGGCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGatat                       <  1:379270/150‑1 (MQ=255)
                                                                                                                                 gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATATTATCTTTCACCAGCGCTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGTGCGATATTGACGc                 >  2:91549/1‑150 (MQ=255)
                                                                                                                                 gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGc                 >  2:534426/1‑150 (MQ=255)
                                                                                                                                 gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGc                 >  1:600843/1‑150 (MQ=255)
                                                                                                                                 gCCCGATAGTCACCAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGc                 >  1:720198/1‑150 (MQ=255)
                                                                                                                                             ccATATATCCAGTGGCCATGCCGTTTTCACCCTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGCTTGCGAACCAAACGAGATAGCCCGACCTAACGCCCCCGCCCCCGCCAGCAATGGCGCGGTATTGTCGCCGCTTTCCGagaa    >  2:868757/1‑151 (MQ=255)
                                                                                                                                             ccAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGacaa    >  1:993303/1‑151 (MQ=255)
                                                                                                                                             ccAAATATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGacaa    >  1:1116349/1‑151 (MQ=255)
                                                                                                                                              cAAATATACAGTGTTCATGCCGTTTTCAAGCGGAAGAAATCCCCCGCTGCTAATATCTTTCCCCAGCGTGTGCGAACCAAACCAGATAGCCAGCCCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGacaac   <  2:902354/151‑1 (MQ=255)
                                                                                                                                                aaTATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGacaaca  >  1:96155/1‑150 (MQ=255)
                                                                                                                                                aaTATCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGacaaca  >  1:908205/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGAACCAATGCCGCGGACAAAGTTAGCAAAGGTGGTTATCGAAGACCACGGAATGATGTGATACGCCCCGGTATCCTGTCGTACGCCCACGGAGCGAATCGACATCCGTTCCACTGTGCCGGTCAACGGCCCGATAGTCACCAAATCTCCAGTGTTCATGCCGTTTTCAAACTGAATAAATACCCCGGTGATAATATCTTTCACCAGCGTTTGCGAACCAAACGAGATAGCCAGACCTAATGCCCCGGCACCTGCCAGCAATGGCGCGATATTGACGCCGATTTCCGACAACA  >  NC_000913/843602‑843895

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