Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 487,092 T→C L186P (CTG→CCG)  mscK → mechanosensitive channel protein, intermediate conductance, K+ regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913487,0920TC99.0% 355.6 / ‑4.0 97L186P (CTG→CCG) mscKmechanosensitive channel protein, intermediate conductance, K+ regulated
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (49/47);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.38e-01

GCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAGC  >  NC_000913/486947‑487238
                                                                                                                                                 |                                                                                                                                                  
gCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATg                                                                                                                                               >  2:916148/1‑151 (MQ=255)
gCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATg                                                                                                                                               <  2:556004/151‑1 (MQ=255)
   cAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATggg                                                                                                                                             >  1:919340/1‑150 (MQ=255)
        gcTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGGTGGGACTGAt                                                                                                                                       <  2:237095/151‑1 (MQ=255)
        gcTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGAt                                                                                                                                       >  1:45592/1‑151 (MQ=255)
        gcTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGAt                                                                                                                                       >  2:978899/1‑151 (MQ=255)
         cTGGACGATTTGAAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTCCGTTACAGACGCAGCAAGAACGCGGGCAAAATGCGATGTATAACGCTTCGCAGCAGCTTCAACAAATTCGCAGTCGTCCGGATGGGACTGAt                                                                                                                                       <  2:1245491/150‑1 (MQ=255)
           ggACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTc                                                                                                                                    <  2:230635/151‑1 (MQ=255)
             aCGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCCTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCgg                                                                                                                                  <  1:789560/151‑1 (MQ=255)
              cGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGc                                                                                                                                 >  1:892282/1‑151 (MQ=255)
              cGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGc                                                                                                                                 >  1:110113/1‑151 (MQ=255)
               gATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCg                                                                                                                                >  1:1016575/1‑151 (MQ=255)
                aTTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCg                                                                                                                                >  2:1180575/1‑150 (MQ=255)
                  ttGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCgag                                                                                                                              <  1:1107277/150‑1 (MQ=255)
                   tGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCgaga                                                                                                                             <  2:276214/150‑1 (MQ=255)
                   tGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGAc                                                                                                                            >  1:705660/1‑151 (MQ=255)
                    gCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGAc                                                                                                                            >  2:1182551/1‑150 (MQ=255)
                     cAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGAc                                                                                                                            >  2:747306/1‑149 (MQ=255)
                      aaaaCGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGc                                                                                                                         >  1:296853/1‑151 (MQ=255)
                       aaaCGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGcc                                                                                                                        >  2:8048/1‑151 (MQ=255)
                         aCGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCtt                                                                                                                      <  2:1173670/151‑1 (MQ=255)
                         aCGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCtt                                                                                                                      <  2:899372/151‑1 (MQ=255)
                         aCGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCtt                                                                                                                      >  2:994892/1‑151 (MQ=255)
                           gCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTAc                                                                                                                    <  2:247789/151‑1 (MQ=255)
                             acaAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACg                                                                                                                   >  1:1128611/1‑150 (MQ=255)
                             acaAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGt                                                                                                                  >  1:298509/1‑151 (MQ=255)
                              caAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGt                                                                                                                  <  1:310440/150‑1 (MQ=255)
                                  aCGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTccc                                                                                                               >  1:371865/1‑149 (MQ=255)
                                         ggCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCccagc                                                                                                            <  1:801497/145‑1 (MQ=255)
                                         ggCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGa                                                                                                        <  2:277216/149‑1 (MQ=255)
                                          gCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAg                                                                                                     <  2:21811/151‑1 (MQ=255)
                                           cGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAgt                                                                                                    >  2:102890/1‑151 (MQ=255)
                                            gTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAgtg                                                                                                   <  1:643954/151‑1 (MQ=255)
                                               ttATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGCGtt                                                                                                 >  2:838287/1‑150 (MQ=255)
                                                  taACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATg                                                                                             >  1:117848/1‑151 (MQ=255)
                                                    aCAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGc                                                                                            >  1:6143/1‑150 (MQ=255)
                                                      agccagcTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAgg                                                                                         >  2:1107278/1‑151 (MQ=255)
                                                      agccagcTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAgg                                                                                         <  1:314179/151‑1 (MQ=255)
                                                      agccagcTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAgg                                                                                         <  2:1011249/151‑1 (MQ=255)
                                                         cagcTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGccc                                                                                      >  2:894067/1‑151 (MQ=255)
                                                          agcTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCca                                                                                     <  2:628312/151‑1 (MQ=255)
                                                              ggTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCcagca                                                                                  <  1:1019928/150‑1 (MQ=255)
                                                                   cGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGt                                                                             <  1:712080/150‑1 (MQ=255)
                                                                     ttACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGc                                                                          >  1:96503/1‑151 (MQ=255)
                                                                       aCAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTg                                                                        <  1:876230/151‑1 (MQ=255)
                                                                         aGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGa                                                                       >  2:179475/1‑150 (MQ=255)
                                                                              cAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCgg                                                                 <  1:961410/151‑1 (MQ=255)
                                                                               aGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGga                                                                <  1:1000374/151‑1 (MQ=255)
                                                                                gCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGgag                                                               <  2:638836/151‑1 (MQ=255)
                                                                                      aCGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGAc                                                         >  1:302257/1‑151 (MQ=255)
                                                                                       cgcgTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGAcc                                                        >  2:553645/1‑151 (MQ=255)
                                                                                        gcgTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGAcc                                                        >  1:356934/1‑150 (MQ=255)
                                                                                         cgTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACcag                                                      <  1:294951/151‑1 (MQ=255)
                                                                                          gTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACcagc                                                     <  2:596636/151‑1 (MQ=255)
                                                                                              aaaaTGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCg                                                 >  2:769960/1‑151 (MQ=255)
                                                                                              aaaaTGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCg                                                 <  1:73975/151‑1 (MQ=255)
                                                                                               aaaTGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGt                                                >  2:1105162/1‑151 (MQ=255)
                                                                                               aaaTGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGt                                                >  2:1084639/1‑151 (MQ=255)
                                                                                                  tGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTaa                                              >  1:1089463/1‑150 (MQ=255)
                                                                                                   gCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTaaa                                             >  1:1236354/1‑150 (MQ=255)
                                                                                                    cGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTaaa                                             >  1:197869/1‑149 (MQ=255)
                                                                                                    cGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTaaa                                             >  1:1285/1‑149 (MQ=255)
                                                                                                    cGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTaaa                                             >  1:548188/1‑149 (MQ=255)
                                                                                                     gATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAg                                           <  2:926554/150‑1 (MQ=255)
                                                                                                        gTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCt                                        <  2:534359/150‑1 (MQ=255)
                                                                                                         tataACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTgg                                      >  2:1199671/1‑151 (MQ=255)
                                                                                                         tataACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTgg                                      >  2:630068/1‑151 (MQ=255)
                                                                                                          ataACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTgg                                      >  1:1044068/1‑150 (MQ=255)
                                                                                                            aaCGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGaa                                    <  1:357239/150‑1 (MQ=255)
                                                                                                            aaCGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGaa                                    <  2:45594/150‑1 (MQ=255)
                                                                                                            aaCGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAg                                   <  1:1169507/151‑1 (MQ=255)
                                                                                                                cTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGa                                >  2:605949/1‑150 (MQ=255)
                                                                                                                  tCGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGaa                               <  1:735316/149‑1 (MQ=255)
                                                                                                                   cGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAacac                            >  2:745938/1‑151 (MQ=255)
                                                                                                                   cGCAGCAGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAacac                            >  1:1206292/1‑151 (MQ=255)
                                                                                                                    gcagcaGCTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACAcc                           >  1:414740/1‑151 (MQ=255)
                                                                                                                      agcagcTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCg                          <  1:902340/150‑1 (MQ=255)
                                                                                                                        cagcTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCt                       >  2:1201257/1‑151 (MQ=255)
                                                                                                                         agcTGCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCtt                      >  2:852205/1‑151 (MQ=255)
                                                                                                                             gCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCa                   <  1:1122346/150‑1 (MQ=255)
                                                                                                                             gCAACAAATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAg                  <  1:762242/151‑1 (MQ=255)
                                                                                                                               aacaaATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAgg                 >  1:587278/1‑150 (MQ=255)
                                                                                                                               aacaaATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGa                <  1:681133/151‑1 (MQ=255)
                                                                                                                               aacaaATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGa                <  2:552075/151‑1 (MQ=255)
                                                                                                                                acaaATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGAt               <  2:340032/151‑1 (MQ=255)
                                                                                                                                 caaATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGAt               >  1:532210/1‑150 (MQ=255)
                                                                                                                                 caaATTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATa              <  1:1182536/151‑1 (MQ=255)
                                                                                                                                    aTTCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCt           <  2:50548/151‑1 (MQ=255)
                                                                                                                                      tCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCtt          <  1:1110932/150‑1 (MQ=255)
                                                                                                                                      tCGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTg         >  1:140516/1‑151 (MQ=255)
                                                                                                                                       cGCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGc        <  1:90863/151‑1 (MQ=255)
                                                                                                                                        gCAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCa       <  2:414974/151‑1 (MQ=255)
                                                                                                                                         cAGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCa       >  1:38929/1‑150 (MQ=255)
                                                                                                                                          aGTCGTCCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCaa      >  2:189597/1‑150 (MQ=255)
                                                                                                                                           gtcgtcCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGACCCAGCAGGCGTTGCTGAATGTGGAGAGTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCTTCTTGCGGGATTCCTTGCaaaa    <  2:720383/151‑1 (MQ=255)
                                                                                                                                             cgtcCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAg   <  2:371852/150‑1 (MQ=255)
                                                                                                                                             cgtcCGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAGc  <  2:96504/151‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCCCAGGCGCTGGACGATTTGCAAAACGCACAAAACGATCTGGCGTCTTATAACAGCCAGCTGGTTTCGTTACAGACGCAGCCCGAACGCGTGCAAAATGCGATGTATAACGCTTCGCAGCAGCTGCAACAAATTCGCAGTCGTCTGGATGGGACTGATGTCGGCGAGACAGCCTTACGTCCCAGCCAGAAAGTGTTAATGCAGGCCCAGCAGGCGTTGCTGAATGCGGAGATTGACCAGCAGCGTAAAAGCCTGGAAGGGAACACCGTCTTGCAGGATACCTTGCAAAAGC  >  NC_000913/486947‑487238

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