Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 860,657 A→G R650R (CGT→CGC ybiW ← putative pyruvate formate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913860,6570AG98.7% 265.3 / ‑4.1 75R650R (CGT→CGCybiWputative pyruvate formate lyase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (42/32);  total (43/32)
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.39e-01

GATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCA  >  NC_000913/860509‑860802
                                                                                                                                                    |                                                                                                                                                 
gATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACcgcg                                                                                                                                                 <  1:616232/151‑1 (MQ=255)
      cAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGt                                                                                                                                            <  1:1346104/150‑1 (MQ=255)
             ggAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCgg                                                                                                                                     >  1:687133/1‑150 (MQ=255)
             ggAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGa                                                                                                                                    <  1:236258/151‑1 (MQ=255)
              gAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGAt                                                                                                                                   <  2:822015/151‑1 (MQ=255)
                  ccGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATg                                                                                                                               <  2:1369723/151‑1 (MQ=255)
                  ccGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATg                                                                                                                               >  2:1173085/1‑151 (MQ=255)
                   cGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATgg                                                                                                                              >  2:1000908/1‑151 (MQ=255)
                     ggCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTa                                                                                                                            <  1:778726/151‑1 (MQ=255)
                      gCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTAc                                                                                                                           >  2:468038/1‑151 (MQ=255)
                          gcgcCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGt                                                                                                                        <  1:699012/150‑1 (MQ=255)
                          gcgcCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGt                                                                                                                        <  1:766863/150‑1 (MQ=255)
                            gcCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTc                                                                                                                     <  2:245144/151‑1 (MQ=255)
                                tCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAg                                                                                                                   <  1:163517/149‑1 (MQ=255)
                                tCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGtt                                                                                                                 >  2:944257/1‑151 (MQ=255)
                                 cTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTtc                                                                                                                >  1:28773/1‑151 (MQ=255)
                                  tGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcg                                                                                                               >  1:1403176/1‑151 (MQ=255)
                                   gCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgt                                                                                                              <  1:1189734/151‑1 (MQ=255)
                                    ccAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                                                                                             >  2:346732/1‑151 (MQ=255)
                                    ccAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTtcgtc                                                                                                             >  2:1452498/1‑151 (MQ=255)
                                       gCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGat                                                                                                          <  2:1302945/151‑1 (MQ=255)
                                         ggggTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGatat                                                                                                        <  2:785384/151‑1 (MQ=255)
                                              gtgGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                                                                    <  2:1383218/150‑1 (MQ=255)
                                              gtgGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                                                                    <  2:175861/150‑1 (MQ=255)
                                               tgGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCt                                                                                                  >  2:1013240/1‑151 (MQ=255)
                                               tgGGCTTTACGCCCGTCCGGAGTTGCCATAGTCTGCGCGCCAGACGGTACGTTAGCGGAGATTGATGACGTACCCGCGGAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGt                                                                                                    >  2:584162/1‑149 (MQ=255)
                                                gggCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTg                                                                                                 <  2:197505/151‑1 (MQ=255)
                                                  gCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                                                               >  1:21203/1‑151 (MQ=255)
                                                   cTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGAt                                                                                               >  1:99741/1‑150 (MQ=255)
                                                     ttACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATa                                                                                              <  2:353863/149‑1 (MQ=255)
                                                      tACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgc                                                                                           >  1:1086623/1‑151 (MQ=255)
                                                        cGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                                                                         <  2:140211/151‑1 (MQ=255)
                                                         gCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGcgcg                                                                                        >  2:136944/1‑151 (MQ=255)
                                                         gCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGcgcg                                                                                        <  1:1444616/151‑1 (MQ=255)
                                                          cccGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAgcgc                                                                                         <  1:605374/149‑1 (MQ=255)
                                                            cGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCTGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGc                                                                                     >  1:71138/1‑151 (MQ=255)
                                                                 ggTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                                                                 >  2:535784/1‑150 (MQ=255)
                                                                 ggTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCcagc                                                                                 >  2:1187376/1‑150 (MQ=255)
                                                                       gCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTAt                                                                          <  2:131453/151‑1 (MQ=255)
                                                                         cATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCg                                                                        <  1:75179/151‑1 (MQ=255)
                                                                                tGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTat                                                                 >  1:767630/1‑151 (MQ=255)
                                                                                 gcgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTat                                                                 <  1:1000818/150‑1 (MQ=255)
                                                                                 gcgcgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTa                                                                  <  1:1173181/149‑1 (MQ=255)
                                                                                    cgcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTatcatc                                                             <  1:590877/151‑1 (MQ=255)
                                                                                     gcCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcg                                                            >  2:1358528/1‑151 (MQ=255)
                                                                                      ccAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcgt                                                           >  2:674727/1‑151 (MQ=255)
                                                                                      ccAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCAtcg                                                            >  2:347136/1‑150 (MQ=255)
                                                                                           cGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTg                                                      >  2:986039/1‑151 (MQ=255)
                                                                                                  ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACt                                               >  2:1181699/1‑151 (MQ=255)
                                                                                                  ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACt                                               >  1:755280/1‑151 (MQ=255)
                                                                                                  ttAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACt                                               >  2:1468227/1‑151 (MQ=255)
                                                                                                     gCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCg                                            >  2:747887/1‑151 (MQ=255)
                                                                                                       ggAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGgc                                          <  1:1454168/151‑1 (MQ=255)
                                                                                                             tgatgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCg                                    >  1:858176/1‑151 (MQ=255)
                                                                                                              gatgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGt                                   <  1:184122/151‑1 (MQ=255)
                                                                                                                tgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGtt                                  >  2:1384225/1‑150 (MQ=255)
                                                                                                                tgaCGTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGtt                                  >  2:1193817/1‑150 (MQ=255)
                                                                                                                    gTACCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCa                             <  1:781873/151‑1 (MQ=255)
                                                                                                                      aCCCGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAg                            <  2:650738/150‑1 (MQ=255)
                                                                                                                        ccGCGTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCg                         <  1:1136723/151‑1 (MQ=255)
                                                                                                                         cgcgTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGc                        >  2:1042699/1‑151 (MQ=255)
                                                                                                                           cgTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTg                      >  1:822693/1‑151 (MQ=255)
                                                                                                                            gTAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCt                       >  1:875021/1‑149 (MQ=255)
                                                                                                                             tAATAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGa                     <  2:964487/150‑1 (MQ=255)
                                                                                                                              aaTAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGAc                    >  2:1467256/1‑150 (MQ=255)
                                                                                                                              aaTAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACg                   >  2:1390820/1‑151 (MQ=255)
                                                                                                                                tAGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGc                  >  2:217801/1‑150 (MQ=255)
                                                                                                                                 aGTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGc                  >  1:354164/1‑149 (MQ=255)
                                                                                                                                  gTTGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAgc               >  1:1384032/1‑151 (MQ=255)
                                                                                                                                   ttGCCGCCAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAgct              >  1:591382/1‑151 (MQ=255)
                                                                                                                                          cAACCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTg       <  1:134736/151‑1 (MQ=255)
                                                                                                                                           aaCCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTg       >  1:710226/1‑150 (MQ=255)
                                                                                                                                           aaCCGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTg       >  2:1125895/1‑150 (MQ=255)
                                                                                                                                              cGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTc   >  2:1387946/1‑151 (MQ=255)
                                                                                                                                              cGGACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTCGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTc   >  1:721610/1‑151 (MQ=255)
                                                                                                                                               ggACCGCGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCa  <  2:655854/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GATGGTCAGTACCGGAGGCCGGGCTTGCGCCTTCTGCCAGCGGGGTGTGGGCTTTACGCCCGTCCGGTGTTGCCATAGTCTGCGCGCCAAACGGTACGTTAGCGGAGATTGATGACGTACCCGCGTAATAGTTGCCGCCAACCGGACCACGACCGTAGCGCGGATTATGGTACTGTTTCAGTTCGTCGATATAGGTCTGATAAGCGCGAGCCAGCAGCGTATCGACAGTATCATCGTCGTTGCCGTACTTCGGCGCACCGTTAATCAGCCGCTGACGCAGCTGCTCGTGAGTCA  >  NC_000913/860509‑860802

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