Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,321,045 A→C L68* (TTA→TGA)  phnH ← ribophosphonate triphosphate synthase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,321,0450AC98.6% 273.3 / ‑4.2 71L68* (TTA→TGA) phnHribophosphonate triphosphate synthase subunit
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (34/35);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

CGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGAC  >  NC_000913/4320897‑4321192
                                                                                                                                                    |                                                                                                                                                   
cGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCAt                                                                                                                                                   <  1:725307/151‑1 (MQ=255)
 ggCGCAACGGCGGTGCCGGGGGAAAGGGCGTTGAGCTGTTTGCTGGAAAACGCCTCATCCGTCACCGCGAAGGGCGCCTGTTCCGGCTGGCTGAACAGCGGCGCGGTGGTATGAAAAAGCCGGCTCTGGTTGGCGGTATCGGTCTTTTAt                                                                                                                                                   >  1:1197494/1‑150 (MQ=255)
  gcgcAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATgg                                                                                                                                                 <  2:653050/151‑1 (MQ=255)
       aCGGCGGGGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGt                                                                                                                                              <  1:499537/149‑1 (MQ=255)
           cggTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGaa                                                                                                                                          >  1:159531/1‑149 (MQ=255)
            ggTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGc                                                                                                                                       >  1:56982/1‑151 (MQ=255)
                ccGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCcaca                                                                                                                                   <  1:538579/151‑1 (MQ=255)
                ccGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCcac                                                                                                                                    >  1:1125316/1‑150 (MQ=255)
                   gTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCg                                                                                                                                >  2:1294306/1‑151 (MQ=255)
                   gTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCg                                                                                                                                >  1:484729/1‑151 (MQ=255)
                     ggAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCgg                                                                                                                               <  2:750226/150‑1 (MQ=255)
                     ggAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGc                                                                                                                              <  2:580490/151‑1 (MQ=255)
                       aaaGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCgt                                                                                                                            >  2:1168086/1‑151 (MQ=255)
                        aaGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCgtg                                                                                                                           <  2:965846/151‑1 (MQ=255)
                          gggCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTc                                                                                                                         <  1:364285/151‑1 (MQ=255)
                             cGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGtt                                                                                                                      >  1:1004142/1‑151 (MQ=255)
                               ttGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTAt                                                                                                                    >  1:294981/1‑151 (MQ=255)
                                tGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTAt                                                                                                                    <  1:331208/150‑1 (MQ=255)
                                  aGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCgg                                                                                                                 >  1:120982/1‑151 (MQ=255)
                                     tGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGcc                                                                                                               >  1:272135/1‑150 (MQ=255)
                                     tGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGcc                                                                                                               >  1:500189/1‑150 (MQ=255)
                                     tGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGcc                                                                                                               >  2:843098/1‑150 (MQ=255)
                                      gTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCa                                                                                                              <  1:421507/150‑1 (MQ=255)
                                          gCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTc                                                                                                         >  2:187708/1‑151 (MQ=255)
                                          gCTGGAAATCGCCTCATCCGTCACCGCGAAGGCCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAAAGCAGGCTCTGGTTGACGAAATCGTTATTTCATGGGGGAGAAAGCCACACCGGCGGGTCGTTATCGGCCAGCGt                                                                                                          >  2:554773/1‑150 (MQ=255)
                                                  tCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCAc                                                                                                 >  2:406126/1‑151 (MQ=255)
                                                   cGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACg                                                                                                <  2:290127/151‑1 (MQ=255)
                                                    gCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGc                                                                                                                <  2:692766/134‑1 (MQ=255)
                                                    gCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGc                                                                                                                >  1:692743/1‑134 (MQ=255)
                                                    gCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCAc                                                                                                 >  1:1091278/1‑149 (MQ=255)
                                                     ccTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCCCTGGTTGACGTTATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGc                                                                                               >  2:693125/1‑150 (MQ=255)
                                                        cATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtgg                                                                                            <  1:191881/150‑1 (MQ=255)
                                                        cATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtg                                                                                             <  1:1276048/149‑1 (MQ=255)
                                                           ccGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtggt                                                                                        <  2:680705/151‑1 (MQ=255)
                                                           ccGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtgg                                                                                         <  2:572223/150‑1 (MQ=255)
                                                           ccGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtgg                                                                                         <  2:270448/150‑1 (MQ=255)
                                                            cGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCtggtggtg                                                                                       >  2:476572/1‑151 (MQ=255)
                                                                 ccGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGat                                                                                  <  1:1206593/151‑1 (MQ=255)
                                                                           cGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGt                                                                        <  2:470303/151‑1 (MQ=255)
                                                                                  tCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGGGGCGATATTCAGCGGGTGCCAGc                                                                 >  2:557364/1‑151 (MQ=255)
                                                                                    cGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCcg                                                               >  2:539176/1‑151 (MQ=255)
                                                                                       ctggctgACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCcgcg                                                             >  1:987401/1‑150 (MQ=255)
                                                                                             gACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGa                                                        <  1:351358/149‑1 (MQ=255)
                                                                                             gACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGa                                                        <  1:637159/149‑1 (MQ=255)
                                                                                              aCCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCt                                                     <  1:1272524/151‑1 (MQ=255)
                                                                                                  gcggcgCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATg                                                 >  2:149293/1‑151 (MQ=255)
                                                                                                     gcgcGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAg                                              <  2:806843/151‑1 (MQ=255)
                                                                                                     gcgcGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAg                                              <  2:848919/151‑1 (MQ=255)
                                                                                                     gcgcGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAg                                              <  2:190282/151‑1 (MQ=255)
                                                                                                       gcgTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGGGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCCGAtaacggg                                            >  2:3306/1‑145 (MQ=255)
                                                                                                        cgTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGACGTTATCGGCCAGCGTCAGCAGCACGCTGGGTGAGGCGATATTCAGCGGCTGCCAGCCGCGTTTGAGCTGATGCAg                                              >  1:887442/1‑148 (MQ=255)
                                                                                                         gTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCg                                          <  1:257993/151‑1 (MQ=255)
                                                                                                         gTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGGTTGAGCTGATGCAGGGCg                                          <  2:632482/151‑1 (MQ=255)
                                                                                                              tATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAAt                                     >  2:219412/1‑151 (MQ=255)
                                                                                                               aTGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATc                                    >  2:400679/1‑151 (MQ=255)
                                                                                                                  aaaaCGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACg                                 >  2:491483/1‑151 (MQ=255)
                                                                                                                    aaCGCAGGCTCTGGTTGACGATATCGTTATGTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGcc                               <  2:703567/151‑1 (MQ=255)
                                                                                                                     aCGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGccc                              <  2:120059/151‑1 (MQ=255)
                                                                                                                      cGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCg                             >  2:628984/1‑151 (MQ=255)
                                                                                                                      cGCAGGCTCTGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCg                             <  1:798523/151‑1 (MQ=255)
                                                                                                                              ctGGTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcgc                       <  2:456309/149‑1 (MQ=255)
                                                                                                                                ggTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGgctcgct                      <  1:129254/148‑1 (MQ=255)
                                                                                                                                ggTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCa                    <  1:204583/150‑1 (MQ=255)
                                                                                                                                ggTTGACGATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCAt                   <  1:936647/151‑1 (MQ=255)
                                                                                                                                       gATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCttt            <  2:816555/151‑1 (MQ=255)
                                                                                                                                       gATATCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCttt            >  2:1249433/1‑151 (MQ=255)
                                                                                                                                        atatCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTa           <  2:751287/151‑1 (MQ=255)
                                                                                                                                          atCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAAc         >  1:571229/1‑151 (MQ=255)
                                                                                                                                           tCGTTATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAAc         >  2:40472/1‑150 (MQ=255)
                                                                                                                                               tATTTCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGcg    >  2:1170589/1‑151 (MQ=255)
                                                                                                                                                 tttCATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGcgac  >  1:1112160/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGGCGCAACGGCGGTGCCGGTGGAAAGGGCGTTGAGCTGTTCGCTGGAAATCGCCTCATCCGTCACCGCGAAGGTCGCCTGTTCCGGCTGGCTGACCAGCGGCGCGTTGGTATGAAAACGCAGGCTCTGGTTGACGATATCGTTATTTAATGGGGTAGAAAGCCACACCGGCGTGTCGTTATCGGCCAGCGTCAGCAGCACGCTGGTGGTGGCGATATTCAGCGGTTGCCAGCCGCGTTTGAGCTGATGCAGGGCGACAATCACGCCCGGCTCGCTCATGGCCTTTAACAGGCGAC  >  NC_000913/4320897‑4321192

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