Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 515,571 A→G I68T (ATC→ACC)  qmcA ← PHB domain membrane‑anchored putative protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913515,5710AG100.0% 439.2 / NA 122I68T (ATC→ACC) qmcAPHB domain membrane‑anchored putative protease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (66/56);  total (66/56)

Only 100 of 124 total aligned reads displayed.

ACGGATGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGATATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACAC  >  NC_000913/515423‑515718
                                                                                                                                                    |                                                                                                                                                   
aCGGATGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGgta                                                                                                                                                   <  4:473125/151‑1 (MQ=255)
  ggATGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATc                                                                                                                                                 <  4:1014697/151‑1 (MQ=255)
   gATGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCg                                                                                                                                                <  4:1313461/151‑1 (MQ=255)
    aTGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGa                                                                                                                                               <  2:62520/151‑1 (MQ=255)
     tGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCg                                                                                                                                                >  4:793603/1‑149 (MQ=255)
      gTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAg                                                                                                                                              <  2:123272/150‑1 (MQ=255)
       ttAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGc                                                                                                                                             >  1:754071/1‑150 (MQ=255)
         aGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGATGACTTCATACGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGCAGGGGTATCGAGCACt                                                                                                                                          <  1:460595/151‑1 (MQ=255)
            cATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGc                                                                                                                                       <  1:1147444/151‑1 (MQ=255)
              tGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTc                                                                                                                                     >  1:1917265/1‑151 (MQ=255)
              tGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTc                                                                                                                                     >  3:450330/1‑151 (MQ=255)
                   aGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCcatca                                                                                                                                <  4:1683162/151‑1 (MQ=255)
                   aGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCcatc                                                                                                                                 >  1:665396/1‑150 (MQ=255)
                        gatgatCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTg                                                                                                                           <  1:121992/151‑1 (MQ=255)
                        gatgatCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTg                                                                                                                           <  4:1399949/151‑1 (MQ=255)
                        gatgatCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTg                                                                                                                           >  4:331614/1‑151 (MQ=255)
                         atgatCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTg                                                                                                                           >  3:342539/1‑150 (MQ=255)
                           gatCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGAt                                                                                                                         <  4:1044943/150‑1 (MQ=255)
                              cGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCtt                                                                                                                      <  1:237162/150‑1 (MQ=255)
                              cGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTg                                                                                                                     <  1:1160740/151‑1 (MQ=255)
                                 cAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCg                                                                                                                   >  4:942177/1‑150 (MQ=255)
                                 cAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCg                                                                                                                   >  1:380855/1‑150 (MQ=255)
                                 cAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCg                                                                                                                   >  4:954043/1‑150 (MQ=255)
                                   gCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGAcc                                                                                                                >  4:683069/1‑151 (MQ=255)
                                   gCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGAcc                                                                                                                >  1:39000/1‑151 (MQ=255)
                                     tCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCa                                                                                                               >  1:1192539/1‑150 (MQ=255)
                                       cAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATg                                                                                                            >  4:1316208/1‑151 (MQ=255)
                                        aGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGc                                                                                                           >  4:806290/1‑151 (MQ=255)
                                             gCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATc                                                                                                       <  1:347817/150‑1 (MQ=255)
                                              cTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATcc                                                                                                      >  1:258882/1‑150 (MQ=255)
                                                gACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCAt                                                                                                    <  1:193592/150‑1 (MQ=255)
                                                    tCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAAc                                                                                                <  3:217177/150‑1 (MQ=255)
                                                     cATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACg                                                                                               <  1:1406060/150‑1 (MQ=255)
                                                      aTAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGc                                                                                             >  1:809108/1‑151 (MQ=255)
                                                        aaGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGcac                                                                                           >  4:963370/1‑151 (MQ=255)
                                                         aGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGcacc                                                                                          >  1:1833969/1‑151 (MQ=255)
                                                         aGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGcacc                                                                                          >  1:1173169/1‑151 (MQ=255)
                                                          gCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGcacca                                                                                         <  4:424950/151‑1 (MQ=255)
                                                             gcgcgcgGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGcaccac                                                                                        >  1:1374747/1‑149 (MQ=255)
                                                             gcgcgcgGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCAccacca                                                                                      <  1:900683/151‑1 (MQ=255)
                                                               gcgcggcgcgTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGa                                                                                    <  1:509070/151‑1 (MQ=255)
                                                                  cggcgcgTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTg                                                                                 <  4:1422519/151‑1 (MQ=255)
                                                                  cggcgcgTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTg                                                                                 >  1:1579655/1‑151 (MQ=255)
                                                                       cgTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGccc                                                                            >  3:500965/1‑151 (MQ=255)
                                                                       cgTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGccc                                                                            >  2:404568/1‑151 (MQ=255)
                                                                         tCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCg                                                                          >  2:56874/1‑151 (MQ=255)
                                                                         tCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCg                                                                          >  4:1483546/1‑151 (MQ=255)
                                                                          cAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCgg                                                                         <  3:117640/151‑1 (MQ=255)
                                                                          cAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCgg                                                                         <  4:1123884/151‑1 (MQ=255)
                                                                          cAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCg                                                                          <  4:312582/150‑1 (MQ=255)
                                                                            aTCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCt                                                                       <  1:249288/151‑1 (MQ=255)
                                                                              cACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTg                                                                      <  1:331502/150‑1 (MQ=255)
                                                                              cACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGt                                                                     >  1:1300379/1‑151 (MQ=255)
                                                                               aCCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTa                                                                    <  1:1454014/151‑1 (MQ=255)
                                                                                ccTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTaa                                                                   <  1:359327/151‑1 (MQ=255)
                                                                                ccTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTaa                                                                   <  4:1363716/151‑1 (MQ=255)
                                                                                 cTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAAc                                                                  >  4:800409/1‑151 (MQ=255)
                                                                                 cTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAAc                                                                  >  4:450320/1‑151 (MQ=255)
                                                                                  tGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACg                                                                 >  1:1738763/1‑151 (MQ=255)
                                                                                     aTAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGttt                                                              >  4:1659030/1‑151 (MQ=255)
                                                                                      tAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGtttt                                                             >  4:1103735/1‑151 (MQ=255)
                                                                                      tAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGtttt                                                             >  1:890913/1‑151 (MQ=255)
                                                                                       aaaGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTg                                                            >  1:550274/1‑151 (MQ=255)
                                                                                       aaaGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTg                                                            >  4:80354/1‑151 (MQ=255)
                                                                                        aaGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTgg                                                           <  1:1821795/151‑1 (MQ=255)
                                                                                        aaGAACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATAGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCACCGGCTGTAACGTTTTgg                                                           >  1:1554803/1‑151 (MQ=255)
                                                                                               cGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCg                                                    >  4:303986/1‑151 (MQ=255)
                                                                                                ggCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGa                                                   <  1:1522220/151‑1 (MQ=255)
                                                                                                  cGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGTTCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGAcc                                                 >  1:902632/1‑151 (MQ=255)
                                                                                                     cGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCaaa                                              <  2:81460/151‑1 (MQ=255)
                                                                                                     cGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCaa                                               >  3:495116/1‑150 (MQ=255)
                                                                                                      gATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCaa                                               >  4:1445967/1‑149 (MQ=255)
                                                                                                         ggTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCg                                           <  4:1764215/150‑1 (MQ=255)
                                                                                                         ggTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGt                                          <  3:191319/151‑1 (MQ=255)
                                                                                                         ggTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGt                                          <  3:491330/151‑1 (MQ=255)
                                                                                                            aaCGTGGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCt                                       <  1:380139/151‑1 (MQ=255)
                                                                                                               gTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACt                                    <  1:793263/151‑1 (MQ=255)
                                                                                                                   gCGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTcc                                >  3:258942/1‑151 (MQ=255)
                                                                                                                    cGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTcc                                >  1:1722376/1‑150 (MQ=255)
                                                                                                                    cGTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCa                               <  1:1110888/151‑1 (MQ=255)
                                                                                                                     gTTATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCAc                              >  4:288842/1‑151 (MQ=255)
                                                                                                                       tATCTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACt                             <  4:525333/150‑1 (MQ=255)
                                                                                                                          cTTTCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGAt                          >  2:150811/1‑150 (MQ=255)
                                                                                                                          cTTTCGAGATAACTTCCTGGGAAGGGGGATCGAGCACTTGCTCTATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCATCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGAt                          >  1:318433/1‑150 (MQ=255)
                                                                                                                            ttCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGa                           <  1:1869279/147‑1 (MQ=255)
                                                                                                                            ttCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGc                       >  4:944683/1‑151 (MQ=255)
                                                                                                                             tCGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGcc                      >  1:1150494/1‑151 (MQ=255)
                                                                                                                              cGAGATAACTTCCTGGGAAGGGGTATTGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGccc                     <  4:1350178/151‑1 (MQ=255)
                                                                                                                              cGAGATAACTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGccc                     >  4:1574343/1‑151 (MQ=255)
                                                                                                                                    aaCTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGt               >  4:566254/1‑151 (MQ=255)
                                                                                                                                    aaCTTCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGt               <  1:1179851/151‑1 (MQ=255)
                                                                                                                                        tCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTAc             <  1:800064/149‑1 (MQ=255)
                                                                                                                                        tCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACg            >  1:856601/1‑150 (MQ=255)
                                                                                                                                        tCCTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACg            >  4:1806557/1‑150 (MQ=255)
                                                                                                                                          cTGGGAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGAtt         >  4:633417/1‑151 (MQ=255)
                                                                                                                                             ggAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTg      <  4:39010/151‑1 (MQ=255)
                                                                                                                                              gAAGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTg      >  2:314936/1‑150 (MQ=255)
                                                                                                                                               aaGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGac    >  2:289773/1‑151 (MQ=255)
                                                                                                                                                aGGGGTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGac    >  1:1131874/1‑150 (MQ=255)
                                                                                                                                                 ggggTATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGacac  >  1:1560340/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACGGATGTTAGTCATGGTCAGGTTGATGATCGCCAGCTCCAGATTGCTGACTTCATAAGCCGCGCGCGGCGCGTCAATCACCTGAATAAAGCACACGGCGTCGATGGTAACGTTGGCGTTATCTTTCGAGATAACTTCCTGGGAAGGGATATCGAGCACTTGCTCCATCATATTGATCTTGCGACCAATGCGATCCATAAACGGCACCACCAGACTGAGCCCCGGCTGTAACGTTTTGGTATAGCGACCAAAGCGTTCTACTGTCCACTGATAGCCCTGCGGTACGATTTTGACAC  >  NC_000913/515423‑515718

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