Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,836,906 C→T 5.4% S278S (AGC→AGT) ‡ ynjB → putative ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,836,9060CT 5.4% 245.1 / 3.3 74S278S (AGC→AGT) ‡ynjBputative ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (38/32);  new base T (2/2);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAA  >  NC_000913/1836760‑1837050
                                                                                                                                                  |                                                                                                                                                
tatCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAgcg                                                                                                                                                 >  1:370885/1‑148 (MQ=255)
        tgtgCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAtttt                                                                                                                                         <  1:606996/148‑2 (MQ=255)
          tgctgcATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCt                                                                                                                                      <  4:261227/149‑1 (MQ=255)
               cATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCaa                                                                                                                                  >  1:172516/1‑148 (MQ=255)
                aTCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGc                                                                                                                                >  3:725551/1‑149 (MQ=255)
                  ccGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCa                                                                                                                               >  2:250702/1‑148 (MQ=255)
                   cGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAg                                                                                                                              <  1:99686/148‑1 (MQ=255)
                        ttGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAg                                                                                                                        <  3:511765/149‑1 (MQ=255)
                             gcgcGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtttt                                                                                                                    >  1:280322/1‑148 (MQ=255)
                                  aaggaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCtt                                                                                                               >  2:49069/1‑148 (MQ=255)
                                     gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                           >  3:397440/1‑149 (MQ=255)
                                     gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                           >  3:250122/1‑149 (MQ=255)
                                     gaaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                           >  4:230155/1‑149 (MQ=255)
                                      aaaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgc                                                                                                           >  2:401164/1‑148 (MQ=255)
                                       aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTcgcg                                                                                                          >  2:460787/1‑148 (MQ=255)
                                       aaggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                         <  4:725551/149‑1 (MQ=255)
                                        aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCa                                                                                                                                   <  1:261892/122‑1 (MQ=255)
                                        aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCa                                                                                                                                   >  2:261892/1‑122 (MQ=255)
                                        aggaTTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGa                                                                                                         >  2:231395/1‑148 (MQ=255)
                                               cccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtttt                                                                                                                    <  2:817254/130‑1 (MQ=255)
                                               cccGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGtttt                                                                                                                    >  1:817254/1‑130 (MQ=255)
                                                 cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttt                                                                                                                                          <  1:617447/106‑1 (MQ=255)
                                                 cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCGCTGGCGAttt                                                                                                                                          >  2:617447/1‑106 (MQ=255)
                                                 cGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGgatgat                                                                                                <  2:172516/148‑1 (MQ=255)
                                                  gCCTTCACCCGTGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCGTGCGTAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTAACAGTTTTGGCTCTCGCGAGGGGATGATTg                                                                                              <  4:757810/149‑1 (MQ=255)
                                                  gCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGAtt                                                                                               <  1:430944/148‑1 (MQ=255)
                                                   ccTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTg                                                                                              >  1:715973/1‑148 (MQ=255)
                                                    cTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGc                                                                                            >  3:353520/1‑149 (MQ=255)
                                                       cACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCaa                                                                                          <  1:250702/148‑1 (MQ=255)
                                                         cccGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGt                                                                                       <  3:118896/149‑1 (MQ=255)
                                                           cgcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGc                                                                                     >  3:157165/1‑149 (MQ=255)
                                                            gcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGc                                                                                     <  1:97472/148‑1 (MQ=255)
                                                            gcgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCa                                                                                    >  4:319469/1‑149 (MQ=255)
                                                             cgcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAt                                                                                   <  4:742691/149‑1 (MQ=255)
                                                              gcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAt                                                                                   <  2:370885/148‑1 (MQ=255)
                                                              gcgGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCAt                                                                                   <  1:788546/148‑1 (MQ=255)
                                                                   tggatgCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATGGGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGAGGATTGGCAACGTGCATTTCGTc                                                                             <  4:771546/149‑1 (MQ=255)
                                                                         ctctGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCAtt                                                                       >  3:194076/1‑149 (MQ=255)
                                                                           ctgctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTc                                                                      >  1:729204/1‑148 (MQ=255)
                                                                             gctgAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCt                                                                    >  2:863675/1‑148 (MQ=255)
                                                                                  aaGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCaa                                                               >  2:740040/1‑148 (MQ=255)
                                                                                     ccGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTg                                                                                                                   >  4:696887/1‑93 (MQ=255)
                                                                                     ccGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTg                                                                                                                   <  3:696887/93‑1 (MQ=255)
                                                                                        gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGaa                                                         >  2:334353/1‑148 (MQ=255)
                                                                                        gCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGaa                                                         <  2:483422/148‑1 (MQ=255)
                                                                                         cacaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATg                                                       <  4:250122/149‑1 (MQ=255)
                                                                                          acaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTc                                                                                >  4:400114/1‑123 (MQ=255)
                                                                                          acaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTc                                                                                <  3:400114/123‑1 (MQ=255)
                                                                                          acaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATg                                                       <  1:166744/148‑1 (MQ=255)
                                                                                          acaTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGc                                                      <  4:323655/149‑1 (MQ=255)
                                                                                            aTTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGcc                                                     <  1:739508/148‑1 (MQ=255)
                                                                                              tGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAg                                                   >  1:901355/1‑148 (MQ=255)
                                                                                               gcgcCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAAt                                                        <  3:826824/142‑1 (MQ=255)
                                                                                               gcgcCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAAt                                                        >  4:826824/1‑142 (MQ=255)
                                                                                                     gTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTgcgg                                           >  3:187889/1‑149 (MQ=255)
                                                                                                       cGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTgcggcg                                         >  3:812650/1‑149 (MQ=255)
                                                                                                        gCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGa                                        >  3:776764/1‑149 (MQ=255)
                                                                                                         cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGa                                        <  1:231395/148‑1 (MQ=255)
                                                                                                         cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGa                                        <  1:740040/148‑1 (MQ=255)
                                                                                                         cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGa                                        <  2:715973/148‑1 (MQ=255)
                                                                                                         cTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGa                                        <  1:859320/148‑1 (MQ=255)
                                                                                                            aCCTTTAACCCCGCGCATGCGCAGCGAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGGTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGt                                    >  4:760687/1‑149 (MQ=255)
                                                                                                               tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAgg                                                                                                       <  4:205816/79‑1 (MQ=255)
                                                                                                               tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAgg                                                                                                       >  3:205816/1‑79 (MQ=255)
                                                                                                               tttAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAg                                  >  1:372346/1‑148 (MQ=255)
                                                                                                                     cccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCaa                           >  3:375355/1‑149 (MQ=255)
                                                                                                                      ccGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAAt                          >  4:495021/1‑149 (MQ=255)
                                                                                                                       cgcgCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTg                                                                                              >  2:396828/1‑80 (MQ=255)
                                                                                                                       cgcgCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTg                                                                                              <  1:396828/80‑1 (MQ=255)
                                                                                                                            aTGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCt                     <  1:460787/148‑1 (MQ=255)
                                                                                                                                 cagcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGctct                >  2:236889/1‑148 (MQ=255)
                                                                                                                                 cagcaAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGctct                >  1:544765/1‑148 (MQ=255)
                                                                                                                                        aaaTCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt        <  4:54447/149‑1 (MQ=255)
                                                                                                                                         aaTCGCCGCTGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGAt        >  2:603815/1‑148 (MQ=255)
                                                                                                                                              ccAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCaa  >  3:236680/1‑149 (MQ=255)
                                                                                                                                                  cggcgATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCAtt                                                                       >  2:304759/1‑76 (MQ=255)
                                                                                                                                                  cggcgATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCAtt                                                                       <  1:304759/76‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
TATCTTGATGTGCTGCATCCGTATTTGTGGCGCGAAGGAAAGGATTTCCCGCCTTCACCCGCGCGGATGGATGCTCTGCTGAAAGCCGGCACATTGCGCCTGTCGCTGACCTTTAACCCCGCGCATGCGCAGCAAAAAATCGCCAGCGGCGATTTGCCTGCAAGCAGTTACAGTTTTGGCTTTCGCGAGGGGATGATTGGCAACGTGCATTTCGTCACCATTCCTGCCAACGCGAATGCCAGTGCTGCGGCGAAGGTAGTTGCCAATTTCCTGCTCTCACCCGATGCGCAA  >  NC_000913/1836760‑1837050

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