Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,892,358 C→T 6.1% G293D (GGT→GAT)  yoaA ← putative ATP‑dependent helicase, DinG family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,892,3580CT 6.1% 237.3 / 4.7 82G293D (GGT→GAT) yoaAputative ATP‑dependent helicase, DinG family
Reads supporting (aligned to +/- strand):  ref base C (30/47);  new base T (2/3);  total (32/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

ACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACT  >  NC_000913/1892212‑1892506
                                                                                                                                                  |                                                                                                                                                    
aCAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGCACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcg                                                                                                                                                    >  1:1287931/1‑149 (MQ=255)
 cAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgc                                                                                                                                                   <  1:259995/149‑1 (MQ=255)
  aaCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgcg                                                                                                                                                  >  1:1551713/1‑149 (MQ=255)
  aaCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACcgcg                                                                                                                                                  >  1:1271022/1‑149 (MQ=255)
  aaCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCAAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACctcg                                                                                                                                                  >  1:722216/1‑149 (MQ=255)
       ggCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAAc                                                                                                                                             <  1:1173519/149‑1 (MQ=255)
         cgcgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGGGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCt                                                                                                                                           <  1:816079/149‑1 (MQ=255)
          gcgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGAGTAGCTAACAGCGCACGCAGGTTACCGCGATAACCTg                                                                                                                                          >  2:411627/1‑149 (MQ=255)
           cgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTgg                                                                                                                                         <  1:655568/149‑1 (MQ=255)
           cgcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTgg                                                                                                                                         <  1:448334/149‑1 (MQ=255)
            gcTCAAATGCCGCATCCAGCAAGGCGGACCGCCCCAGTGACAGTTTAGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAAGCTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGc                                                                                                                                        <  2:1568793/149‑1 (MQ=255)
            gcTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcg                                                                                                                                        <  1:1040008/149‑3 (MQ=255)
             cTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGACTTTGCGGATTAGCTAACAGCTCACGCAGGTTATCGCGGTAACCTgcgt                                                                                                                                       <  1:453776/149‑4 (MQ=255)
                aaaTGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAc                                                                                                                                    >  1:600777/1‑149 (MQ=255)
                  aTGCCGCATCCAGCACGGCGGAACGCCCCAGTGCCCGTTTCGCGACGACATAAGAACGGTCCAGGGTGTCATCGAGCAGGAAATATGCCGGCAGATTGTGGGGATTAGCATACACATCACGCAGGTATCGGCGATAACCGTGCTcccgg                                                                                                                                  >  2:883102/1‑145 (MQ=255)
                     ccGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGt                                                                                                                               <  1:523089/149‑1 (MQ=255)
                        cATCCAGCAATGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGc                                                                                                                            <  1:1405176/149‑1 (MQ=255)
                        cATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGc                                                                                                                            <  1:1423556/149‑1 (MQ=255)
                         aTCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGCTCTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCTGATTAGCTAACAGCTCACGCAGGTTACCGCTATAACCTGGCTCAGCGAGTTGCa                                                                                                                           >  1:846126/1‑149 (MQ=255)
                             aGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACTTCATAACAAAGTTCCAGGGTGTCATCGAGCATTAAAAATTCACGCTGAATTTGCGGATTAGCTAACAGCTCAAGCAGGTTACCGAGATAACCTGGCTCCACGAGATCAAGACg                                                                                                                       >  2:1015395/1‑149 (MQ=255)
                             aGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACg                                                                                                                       >  1:835307/1‑149 (MQ=255)
                              gCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGCATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGa                                                                                                                      >  2:255556/1‑149 (MQ=255)
                                aaGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGaaa                                                                                                                    <  2:1206330/149‑1 (MQ=255)
                                    cGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATcc                                                                                                                <  1:413056/149‑1 (MQ=255)
                                       aaCGCCCCAGTGACAGTATCGCCACGTCATAACAAAGTTCCACGGTGTCATCGAGCAGTAAAAATGCCCGGTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTgc                                                                                                             <  2:832578/149‑1 (MQ=255)
                                             ccAGTGCCAGTTTCGCCACGTCATAAGATAGTTCCAGGGTGTCATCGAGCAGTAATAATGACCGGTGAATTTGCGGATTAGATAACAGCTCCCGCAGGTTACCGCGATAACCTGGCTCACAGAGTTGCAGCCGAAATTCCTGCGCAgtc                                                                                                       >  2:1052378/1‑149 (MQ=255)
                                              cAGTGACAGTTTCGCCACGTCATATCAAAGTTGCAGGGTGTCATCGAGCAGGAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCAGCGAGTTGCAGACGAAAATCCTGCGCActct                                                                                                      <  1:300000/149‑1 (MQ=255)
                                               aGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTg                                                                                                     <  1:1281707/149‑1 (MQ=255)
                                                gTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAGAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTgg                                                                                                    <  2:216530/149‑1 (MQ=255)
                                                gTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTgg                                                                                                    <  1:963264/149‑1 (MQ=255)
                                                 tGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTggg                                                                                                   >  1:429260/1‑149 (MQ=255)
                                                 tGACAGTATCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTggg                                                                                                   <  2:1309732/149‑1 (MQ=255)
                                                      gTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGa                                                                                              <  2:1287931/149‑1 (MQ=255)
                                                          cGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt                                                                                          <  1:461018/149‑1 (MQ=255)
                                                          cGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGATGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGAt                                                                                          <  1:941146/149‑1 (MQ=255)
                                                            ccACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCa                                                                                        <  1:827845/149‑1 (MQ=255)
                                                                 tCATAACACAGTTCCAGGGTGTCATCGAGCAGTAAAACTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCa                                                                                   <  1:153485/149‑1 (MQ=255)
                                                                 tCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCa                                                                                   <  2:451887/149‑1 (MQ=255)
                                                                  cATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCAc                                                                                  >  1:53118/1‑149 (MQ=255)
                                                                    tAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACtt                                                                                >  2:1102066/1‑149 (MQ=255)
                                                                       cAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTg                                                                             <  2:420772/149‑1 (MQ=255)
                                                                        aaaGTTCCAGGGTGTCATCGAGCAGTAAAAAAGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGt                                                                            <  2:1023408/149‑1 (MQ=255)
                                                                           gTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                                                                         <  1:783992/149‑1 (MQ=255)
                                                                           gTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAAc                                                                         <  2:1138718/149‑1 (MQ=255)
                                                                             tCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAActg                                                                       >  2:403637/1‑149 (MQ=255)
                                                                                aGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAActgctg                                                                    >  1:1469084/1‑149 (MQ=255)
                                                                                 gggTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCAAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTgg                                                                   <  2:745558/149‑1 (MQ=255)
                                                                                     gtCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCTGACGAAAATCCTGCGCACTCTGGGCAAGACGCTCAGCGCACTTTTGTAACTGCTGGgtgt                                                               >  2:146860/1‑149 (MQ=255)
                                                                                       cATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCt                                                             >  2:620141/1‑149 (MQ=255)
                                                                                         tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGCTCAGCGCACTTTTGTAACTGATGGGTGTAttt                                                           >  1:1034759/1‑149 (MQ=255)
                                                                                         tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt                                                           <  2:429260/149‑1 (MQ=255)
                                                                                         tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt                                                           <  2:835307/149‑1 (MQ=255)
                                                                                         tCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCttt                                                           <  2:248607/149‑1 (MQ=255)
                                                                                                 tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGACTTTTAATTCgg                                                   >  2:66241/1‑149 (MQ=255)
                                                                                                 tAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGACTTGCAGACGAAAATCCTGCGCACTATGGTCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCgg                                                   >  2:132096/1‑149 (MQ=255)
                                                                                                  aaaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTGACCGAGTTGCAGACGAAAATCCTGCGCACTCTGTGCAAGTCGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTATTTCGGt                                                  <  2:576959/149‑1 (MQ=255)
                                                                                                  aaaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTACCTGCTGGGTGTCTTTGAATTCGGt                                                  >  2:634444/1‑149 (MQ=255)
                                                                                                   aaaaTGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTg                                                 <  1:1317794/149‑1 (MQ=255)
                                                                                                         cccGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAg                                           <  1:623838/149‑1 (MQ=255)
                                                                                                            gCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCg                                        >  2:649183/1‑149 (MQ=255)
                                                                                                                aaTTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATgg                                    <  2:53118/149‑1 (MQ=255)
                                                                                                                  tttGCGTATTATCTAACAGCTCACGCATGTTACCGCGATAACCTGGCTCACCGAGTTGCAGAAGAATATCCTGCGCACTCTGTGCAAGACTATCAGCTCACTTTTGTTACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTg                                  <  2:722216/149‑1 (MQ=255)
                                                                                                                  tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTg                                  <  1:768465/149‑1 (MQ=255)
                                                                                                                  tttGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGAAAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTg                                  <  2:419177/149‑1 (MQ=255)
                                                                                                                    tGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGAt                                <  2:600777/149‑1 (MQ=255)
                                                                                                                    tGCGGATTAGCTAACACCTAACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGAt                                <  2:365863/149‑1 (MQ=255)
                                                                                                                     gCGGATTAGCTAACAGCTCACGCAGGTTACGGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATg                               <  1:1558737/149‑1 (MQ=255)
                                                                                                                        gATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTCCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCt                            >  1:1533835/1‑149 (MQ=255)
                                                                                                                           tAGCTAACAGATCACGCTGGTTACCGCTATATCCTGGCTGACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGTCAAGACGATCAGCTCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTTATGTCtttt                         <  1:1185758/149‑1 (MQ=255)
                                                                                                                              cTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGcc                      <  1:497379/149‑1 (MQ=255)
                                                                                                                               tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCa                     >  1:1427126/1‑149 (MQ=255)
                                                                                                                               tAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGGTGCTGACGAAAATCCTGCGCACTCTGGGCAAGACGATCATCGCACTTTTGTAACTGCTGGGTGTCTGTTAATTCGGTGCGGTAGGCGATGGTGAGTTCTTTTTCCa                     >  1:1127237/1‑149 (MQ=255)
                                                                                                                                aaCAGCTCACGCAGGTTACCGCGATAACCTGGCTCCCCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAg                    >  1:952887/1‑149 (MQ=255)
                                                                                                                                  cAGCTAACGCAGGTTTCCGCGATAACCTGTCTCACCGAGTTGCAGTCGAAACTCCTGCGCACTCTGGGAAAGATGATCAGCGCACTTTTTTAACTGCTTGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTTATGTCTTTTGACATGt                  <  2:846126/149‑1 (MQ=255)
                                                                                                                                   cgCTCACGCAGGCTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTc                 <  2:1551713/148‑1 (MQ=255)
                                                                                                                                     cTCACGAAGTGTACCGCGATTACCTGGCTCACCGAGCTGCAGCCGAAACTCCTGCGCACTCTGTGCAAGACGTTCAGCGCACTTTTGTACCTGCTGGGTGGCTTTTAATTCTGTGCGGTAGGCGATGGTGATGGCTTTTGCCAGGTCGa               <  2:771191/149‑1 (MQ=255)
                                                                                                                                      tCACGCACTTTACCGCGTTAACCTGGCTCACCTAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGCTCAGCGCACTTTTGTAACTGCTGGGTTTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAg              <  1:1084330/149‑1 (MQ=255)
                                                                                                                                        agccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGGACTTTTGTAACTGCTGGGTGTCTGTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGACAGGTCGAgga            >  1:1262096/4‑147 (MQ=255)
                                                                                                                                        aCGCAGGTTACTGAGATAACCTTGCTCACCGATTTGCAGACTAAAATACTGCGCACTCTGGGCAAGTCGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTTATTCGGTGCGGTAGGGGATGGTGATGTCTTTTGCCAGGTGGTGCa            <  1:1015395/149‑1 (MQ=255)
                                                                                                                                        aCGCAGGTTACCGCGATAACCTGGCTCCCCGAGTTCCAGACGAAAATCCTGCGCACTCTGGTCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGAGGTGGATGTCTTTTGCCAGGTCGAGCa            >  1:818927/1‑149 (MQ=255)
                                                                                                                                         gccAGGTTATCGCGGTAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAg           >  1:176916/3‑149 (MQ=255)
                                                                                                                                             ggTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTATGGGCAAGACGATCAGCGCACTTTTGTTACTGCTGGGTGTCTGATAATTATGTGAGGTAGGCGATGGTGATGTATTTTGCCAGGTCGAGAAGTTGt       >  1:1073459/1‑149 (MQ=255)
                                                                                                                                                  ccGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGTGTGTCTTTTATTTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACt  <  2:1469084/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGGGTGTCATCGAGCAGTAAAAATGCCCGCTGAATTTGCGGATTAGCTAACAGCTCACGCAGGTTACCGCGATAACCTGGCTCACCGAGTTGCAGACGAAAATCCTGCGCACTCTGGGCAAGACGATCAGCGCACTTTTGTAACTGCTGGGTGTCTTTTAATTCGGTGCGGTAGGCGATGGTGATGTCTTTTGCCAGGTCGAGCAGTTGTCGACT  >  NC_000913/1892212‑1892506

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