Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 858,664 G→T 5.0% S133* (TCA→TAA)  ybiU ← DUF1479 domain‑containing protein YbiU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913858,6640GT 5.0% 473.8 / 7.7 139S133* (TCA→TAA) ybiUDUF1479 domain‑containing protein YbiU
Reads supporting (aligned to +/- strand):  ref base G (67/65);  new base T (2/5);  total (69/70)
Fisher's exact test for biased strand distribution p-value = 4.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 144 total aligned reads displayed.

TGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGT  >  NC_000913/858515‑858813
                                                                                                                                                     |                                                                                                                                                     
  aCGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGaa                                                                                                                                                     >  2:2046216/1‑150 (MQ=255)
       cacaTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCt                                                                                                                                                >  2:741652/1‑150 (MQ=255)
         caTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGa                                                                                                                                              >  2:72466/1‑150 (MQ=255)
         caTCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGa                                                                                                                                              >  1:784311/1‑150 (MQ=255)
           tCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGc                                                                                                                                            <  2:1102132/150‑1 (MQ=255)
              gggTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTc                                                                                                                                         >  2:244184/1‑150 (MQ=255)
               ggTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTcc                                                                                                                                        <  1:993719/150‑1 (MQ=255)
                   aaaCCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGCTCAGCGCTTCACGTCCCgag                                                                                                                                    <  1:32784/150‑2 (MQ=255)
                      ccATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCgaagaa                                                                                                                                 >  2:526118/1‑150 (MQ=255)
                        aTTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAt                                                                                                                               <  2:557502/150‑1 (MQ=255)
                         ttGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAAtt                                                                                                                              <  1:2046216/150‑1 (MQ=255)
                          tGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGAACGTGAAGCGCTGAGCGTCCCGAAGAAATTa                                                                                                                             >  2:1355182/1‑150 (MQ=255)
                            cTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTAt                                                                                                                            >  2:429484/1‑149 (MQ=255)
                            cTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATc                                                                                                                           >  2:2041677/1‑150 (MQ=255)
                             tttCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCg                                                                                                                          >  1:1867410/1‑150 (MQ=255)
                                ccATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCg                                                                                                                       >  2:1404652/1‑150 (MQ=255)
                                 cATCACTTTCAAATGTCCACAGACGATTGAGAAACGACAGCGCATTCGCCATTTCTACACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgg                                                                                                                      >  1:1223261/1‑150 (MQ=255)
                                  aTCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCggg                                                                                                                     <  1:1963514/150‑1 (MQ=255)
                                   tCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgggg                                                                                                                    >  2:1255679/1‑150 (MQ=255)
                                   tCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCgggg                                                                                                                    >  1:379385/1‑150 (MQ=255)
                                     aCTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTGCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGcc                                                                                                                  <  1:907245/150‑1 (MQ=255)
                                      cTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCt                                                                                                                 >  1:2147675/1‑150 (MQ=255)
                                         tCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTg                                                                                                              >  2:2158698/1‑150 (MQ=255)
                                           aaaTGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTa                                                                                                            <  2:39155/150‑1 (MQ=255)
                                                 ccACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTg                                                                                                              >  2:1101340/1‑142 (MQ=255)
                                                 ccACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTc                                                                                                      >  1:2113964/1‑150 (MQ=255)
                                                  cacaGACGATTGAGAAACGACTGCGTATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCg                                                                                                     >  1:1110444/1‑150 (MQ=255)
                                                  cacaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCg                                                                                                     <  1:2184416/150‑1 (MQ=255)
                                                   acaGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGt                                                                                                    <  2:457395/150‑1 (MQ=255)
                                                    caGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTc                                                                                                   <  1:2395064/150‑1 (MQ=255)
                                                     aGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCa                                                                                                  <  2:2394692/150‑1 (MQ=255)
                                                         gATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGAACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGc                                                                                              <  1:1355182/150‑1 (MQ=255)
                                                           ttGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAATTATCGGGACGTGAAGCGCTGCGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGCCCTCGTCAAGGCgg                                                                                            >  1:1462270/1‑150 (MQ=255)
                                                               gaAACGACTGCGCATTCGCCATTTCTTCCCTCTGGCGGGCCTGCATTTGCGCCTGCTACCTGTAGATGGGGTAAATCTCGGGACGTGAAGCTCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCGTTGTAGACCTCGTCAAAGCGGTTGc                                                                                        <  1:1458941/150‑1 (MQ=255)
                                                               gaAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGc                                                                                        <  2:648455/150‑1 (MQ=255)
                                                                aaaCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCg                                                                                       <  2:379385/150‑1 (MQ=255)
                                                                  aCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTATATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGt                                                                                     <  2:586133/150‑1 (MQ=255)
                                                                  aCGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGt                                                                                     <  1:2377345/150‑1 (MQ=255)
                                                                   cgccTGCGCATTCGCCATTTCTGCACTGTGGCGGGCCTGCTTTTGCGCGTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                                    <  1:610757/147‑1 (MQ=255)
                                                                   cGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTc                                                                                    <  1:1068938/150‑1 (MQ=255)
                                                                     aCTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCa                                                                                  >  1:1497986/1‑150 (MQ=255)
                                                                      cTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTATATCTCGGTACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAg                                                                                 <  2:2128520/150‑1 (MQ=255)
                                                                      cTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAg                                                                                 >  2:467612/1‑150 (MQ=255)
                                                                       ttcgcATTGGCCATTTCTTCACTCTGGCGGGCCTGGATTTGCGCGTGCGACCAGTAGATGGGTTAAATCTCGGGAGGTGAAGCGCTTAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGTTTGCGGTCCAGa                                                                                <  2:1145097/148‑1 (MQ=255)
                                                                        gcgcATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGAt                                                                               >  1:387575/1‑150 (MQ=255)
                                                                           cATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGTGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCTGTCCAGATAGt                                                                            <  2:1637273/150‑1 (MQ=255)
                                                                              tCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGATGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCa                                                                         <  1:2356819/150‑1 (MQ=255)
                                                                               cGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGACAGTCCAg                                                                        >  1:2058156/1‑150 (MQ=255)
                                                                                 ccATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCa                                                                      >  1:1260578/1‑150 (MQ=255)
                                                                                   aTTTCTTCACTCTGGGGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATc                                                                    <  2:2275619/150‑1 (MQ=255)
                                                                                    tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCGCTAAGAAATTATCTCGGGGGCCTTTGTAGTCCTCTTCAAAGCTGTTGCTGTCCATATAGTACTGGATCg                                                                   <  2:335539/150‑1 (MQ=255)
                                                                                    tttCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCg                                                                   <  1:716224/150‑1 (MQ=255)
                                                                                       cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                >  1:2225770/1‑150 (MQ=255)
                                                                                       cttcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACt                                                                >  2:1974565/1‑150 (MQ=255)
                                                                                         tcACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGa                                                              >  2:1813411/1‑150 (MQ=255)
                                                                                            ctctGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATcc                                                           >  1:879213/1‑150 (MQ=255)
                                                                                                ggCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGc                                                       >  2:1675402/1‑150 (MQ=255)
                                                                                                  cGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCt                                                     <  1:1404652/150‑1 (MQ=255)
                                                                                                     gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgc                                                  >  2:934621/1‑150 (MQ=255)
                                                                                                     gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgc                                                  >  2:456362/1‑150 (MQ=255)
                                                                                                     gCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgc                                                  <  1:1415301/150‑1 (MQ=255)
                                                                                                      ccTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                 <  1:154748/150‑1 (MQ=255)
                                                                                                      ccTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcg                                                 >  2:2048513/1‑150 (MQ=255)
                                                                                                       cTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGAGCTGAGCGTCCCGAAGAAAGTATCGCCGGGGCCTTTGTAGACCTCGTCAACGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAgcgg                                                >  2:1257542/1‑149 (MQ=255)
                                                                                                           aTTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGt                                            <  1:1557088/150‑1 (MQ=255)
                                                                                                            tttGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGtt                                           <  1:1623465/150‑1 (MQ=255)
                                                                                                            tttGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGtt                                           <  2:2147675/150‑1 (MQ=255)
                                                                                                            attGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATAATCGCCGGGGCCTTGGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGAACCCAGCCTAGCGCTTggt                                           >  2:1891321/2‑148 (MQ=255)
                                                                                                               gcgcCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTcgc                                        >  1:1074570/1‑150 (MQ=255)
                                                                                                                  ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                     <  2:241932/150‑1 (MQ=255)
                                                                                                                  ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                     <  1:2998/150‑1 (MQ=255)
                                                                                                                  ccTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggg                                     >  2:2066587/1‑150 (MQ=255)
                                                                                                                    tGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGCCCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCggggg                                   >  2:776237/1‑150 (MQ=255)
                                                                                                                      cGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGaa                                 >  1:1368544/1‑150 (MQ=255)
                                                                                                                        aCCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAAt                               <  2:1497986/150‑1 (MQ=255)
                                                                                                                         ccAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATg                              <  1:1740870/150‑1 (MQ=255)
                                                                                                                          cAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATgg                             >  2:584671/1‑150 (MQ=255)
                                                                                                                           aGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGc                            >  1:1158722/1‑150 (MQ=255)
                                                                                                                            gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGcc                           <  2:176620/150‑1 (MQ=255)
                                                                                                                            gTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTACAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGcc                           <  1:495984/150‑1 (MQ=255)
                                                                                                                             tAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCt                          <  2:387575/150‑1 (MQ=255)
                                                                                                                                 tGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGATTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTa                      >  2:1972941/1‑150 (MQ=255)
                                                                                                                                 tGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTa                      <  1:2245301/150‑1 (MQ=255)
                                                                                                                                    ggTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCa                   <  2:1260578/150‑1 (MQ=255)
                                                                                                                                     gTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAc                  <  1:1256315/150‑1 (MQ=255)
                                                                                                                                      tAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCAcc                 <  1:1675402/150‑1 (MQ=255)
                                                                                                                                        aaTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCATTTATCACCGc               >  1:541495/1‑150 (MQ=255)
                                                                                                                                        aaTCTCGGGACGTGAAGCGCGGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGc               <  2:2029693/150‑1 (MQ=255)
                                                                                                                                         aTCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGca              >  1:530090/1‑150 (MQ=255)
                                                                                                                                          tctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcac             >  1:603645/1‑150 (MQ=255)
                                                                                                                                          tctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcac             <  2:784311/150‑1 (MQ=255)
                                                                                                                                          tctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcac             <  2:1477539/150‑1 (MQ=255)
                                                                                                                                           ttcGGGACGCTCAGCGCTTCACGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca            >  1:1702502/2‑150 (MQ=255)
                                                                                                                                           ctcGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGcaca            >  1:1204002/1‑150 (MQ=255)
                                                                                                                                              gggACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGAGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAAcc         >  2:909987/1‑150 (MQ=255)
                                                                                                                                                    tGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACg   >  1:55395/1‑150 (MQ=255)
                                                                                                                                                    tGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACg   <  1:1255679/150‑1 (MQ=255)
                                                                                                                                                    tGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACg   >  1:1078373/1‑150 (MQ=255)
                                                                                                                                                     gAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGt  >  1:47400/1‑150 (MQ=255)
                                                                                                                                                     gAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGt  >  2:1357957/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGACGCTCACATCCGGGTTAAACCATTGCTTTCCATCACTTTCAAATGTCCACAGACGATTGAGAAACGACTGCGCATTCGCCATTTCTTCACTCTGGCGGGCCTGCATTTGCGCCTGCGACCAGTAGATGGGGTAAATCTCGGGACGTGAAGCGCTGAGCGTCCCGAAGAAATTATCGCCGGGGCCTTTGTAGACCTCGTCAAAGCGGTTGCGGTCCAGATAGTCCAGCATCGACTGATCCCAGCCTAGCGCTTGTTCGCGGGGGAAATGGCCTTTTATCACCGCACAACCGCGACGT  >  NC_000913/858515‑858813

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