Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,897 C→T 18.2% M13I (ATG→ATA) ‡ yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8970CT18.2% 347.9 / 59.7 153M13I (ATG→ATAyaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (69/56);  new base T (20/8);  total (89/64)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 156 total aligned reads displayed.

CAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383801‑383996
                                                                                                |                                                                                                   
 aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACTGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTg                                                                                                 <  1:3688369/100‑1 (MQ=255)
  gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                <  2:476301/100‑1 (MQ=255)
  gACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                >  1:2051141/1‑100 (MQ=255)
     acGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTAt                                                                                             >  2:2084227/1‑100 (MQ=255)
      cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTa                                                                                              <  2:2264618/98‑1 (MQ=255)
      cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            >  1:3353907/1‑100 (MQ=255)
      cGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                            >  2:2817984/1‑100 (MQ=255)
        ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTGTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggt                                                                                          <  2:2040618/100‑2 (MQ=255)
        ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                          <  1:479810/100‑1 (MQ=255)
        ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                          >  1:2945243/1‑100 (MQ=255)
        ggCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGa                                                                                          >  2:3580856/1‑100 (MQ=255)
         gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtt                                                                                         <  1:3127982/100‑3 (MQ=255)
         gCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                         <  2:1091980/100‑1 (MQ=255)
          ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        <  2:2910090/100‑1 (MQ=255)
          ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        >  1:1806178/1‑100 (MQ=255)
          ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACa                                                                                        <  1:2972229/100‑1 (MQ=255)
           cTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttt                                                                                       <  2:2842917/100‑5 (MQ=255)
             tCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttta                                                                                     <  1:2782628/100‑7 (MQ=39)
             tCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggttttt                                                                                      <  1:2358598/99‑6 (MQ=39)
             tCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACTTAc                                                                                     >  1:2596938/1‑100 (MQ=255)
               gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAg                                                                                   >  1:1071399/1‑100 (MQ=255)
                ttATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttttaagt                                                                                  <  1:781213/100‑10 (MQ=38)
                  ataCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                >  1:739954/1‑100 (MQ=255)
                  ataCGCAGGTATAATGGCGTCAATACAGCCTTTCGTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATa                                                                                <  2:3155328/100‑1 (MQ=255)
                   taCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaa                                                                               >  2:3029255/1‑100 (MQ=255)
                   taCGCAGGTATAATAGAGTCAATACACCCTTTCTTTATACTGGGTTGACGTCTTAGCGCACTTAAAAACCATAGGAACATTGCTATGGACATTCAGATaa                                                                               >  2:2222526/1‑100 (MQ=255)
                     cGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAat                                                                             >  2:1094148/1‑100 (MQ=255)
                      gCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAata                                                                            <  1:3376330/100‑1 (MQ=255)
                         ggTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATcc                                                                         >  1:101231/1‑100 (MQ=255)
                                 ggCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTc                                                                 >  2:609308/1‑100 (MQ=255)
                                   cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               <  2:336101/100‑1 (MQ=255)
                                   cGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCAt                                                               >  1:3524068/1‑100 (MQ=255)
                                    gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATg                                                              <  1:3062952/100‑1 (MQ=255)
                                    gTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATg                                                              >  1:2210703/1‑100 (MQ=255)
                                         tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAtt                                                         >  1:1309417/1‑100 (MQ=255)
                                         naCAGCCTTTCTTTTGAATGAGTTTACGGCTTAGCGTCCTTAACTACCTTAGCAACCTTGCTATGCGCCTACAGCTCACTATTCTGACTTTCATGCagt                                                          >  2:2445547/2‑97 (MQ=255)
                                          aCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAAttt                                                        >  1:2618868/1‑100 (MQ=255)
                                           cAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTg                                                       >  1:2991332/1‑100 (MQ=255)
                                             gCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTg                                                       >  1:2612853/1‑98 (MQ=255)
                                             gCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                     >  2:2254164/1‑100 (MQ=255)
                                              ccTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTc                                                    <  1:1718462/100‑1 (MQ=255)
                                                tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaa                                                  >  2:315431/1‑100 (MQ=255)
                                                tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaa                                                  >  1:793330/1‑100 (MQ=255)
                                                 ttctttTGACTGGGTATACGTCTTAGCGTACTTAAAAACCATAGCACCATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                 <  1:1744622/100‑1 (MQ=255)
                                                  tctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                <  2:3353881/100‑1 (MQ=255)
                                                  tctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAg                                                <  1:2543685/100‑1 (MQ=255)
                                                    ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGaa                                              >  1:2339533/1‑100 (MQ=255)
                                                    ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGaa                                              >  2:1298452/1‑100 (MQ=255)
                                                     tttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGa                                               <  2:2319750/98‑1 (MQ=255)
                                                     tttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAAt                                             >  1:617846/1‑100 (MQ=255)
                                                        gACTGGGTGTACGGCGTAGCGTACTTAAAATCCATAGCAGCATTGCTACGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGa                                               <  2:3022149/95‑1 (MQ=255)
                                                            gggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTa                                      >  2:3284207/1‑100 (MQ=255)
                                                            gggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTa                                      <  1:1795831/100‑1 (MQ=255)
                                                            gggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTa                                      <  1:2451730/100‑1 (MQ=255)
                                                             ggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATCTTa                                      >  2:821076/1‑99 (MQ=255)
                                                              gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  2:965754/100‑1 (MQ=255)
                                                              gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  1:1238826/100‑1 (MQ=255)
                                                              gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  2:467851/100‑1 (MQ=255)
                                                              gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATg                                    <  1:583344/100‑1 (MQ=255)
                                                               tttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGa                                   >  1:57448/1‑100 (MQ=255)
                                                               tttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGa                                   <  1:745842/100‑1 (MQ=255)
                                                                 tACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                 <  2:252819/100‑1 (MQ=255)
                                                                 tACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAt                                 >  2:2085171/1‑100 (MQ=255)
                                                                  aCGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAAtt                                <  1:802804/100‑1 (MQ=255)
                                                                     tCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACa                             >  1:901452/1‑100 (MQ=255)
                                                                             tACTTAAACACCATAGCAACATTGCTATGGACATCCCGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAAAATCATTTTATGAATTACAAAGCCTTc                     >  2:160129/1‑100 (MQ=255)
                                                                             tACTTAAAAACCATAGCATCATTGCTATGTACATACAGATAAATATCCAGGTTTTCATGAAATTTGCTCAAAGAATCATTTTATGAATTACAATGCCttt                     <  2:311072/100‑2 (MQ=255)
                                                                              aCTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCa                    >  2:2452939/1‑100 (MQ=255)
                                                                               ctgtatgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  1:1447762/9‑100 (MQ=38)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCCATTTGCTCAAAGAATCATTTTATGAACTACAAAGCCTTCAc                   >  2:747831/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  2:1562369/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  2:1679876/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  1:1244509/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  2:3742018/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  1:1473915/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAc                   >  1:1963834/1‑100 (MQ=255)
                                                                               cTTAAAAACCATAGCAACATTGCTAGGTACATCCAGAAAAATATACAAATTTTCTTGAAATTTGCTCAAAGAATCATTTTATGAATTACAAATCATTCa                    >  2:1894269/1‑99 (MQ=255)
                                                                                tgtatgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATACAATTTGACCAAAGCAACATTTTATGCATTACAAAGCCTTCAcc                  >  2:2297860/8‑100 (MQ=25)
                                                                                 tAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAccc                 >  1:845602/1‑100 (MQ=255)
                                                                                  aaaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCa                <  1:3452797/100‑1 (MQ=255)
                                                                                   aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAg               >  1:2543283/1‑100 (MQ=255)
                                                                                   aaaaCCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAg               >  1:2537124/1‑100 (MQ=255)
                                                                                    tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  2:8710/4‑100 (MQ=255)
                                                                                    tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  1:3509233/4‑100 (MQ=255)
                                                                                    tgtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGa              >  2:2756512/4‑100 (MQ=255)
                                                                                    tgtacATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGATCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAg               >  2:3631150/5‑99 (MQ=255)
                                                                                     gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  1:1678680/3‑100 (MQ=255)
                                                                                     gtccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGAt             >  2:2713010/3‑100 (MQ=255)
                                                                                      tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  1:3598071/2‑100 (MQ=255)
                                                                                      tccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGATCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc            >  2:1319271/2‑100 (MQ=255)
                                                                                       ccATAGCAATGTTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCg           >  1:1327463/1‑100 (MQ=255)
                                                                                        cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc          <  1:598641/100‑1 (MQ=255)
                                                                                        cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc          <  2:1692981/100‑1 (MQ=255)
                                                                                        cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc          <  1:2074540/100‑1 (MQ=255)
                                                                                        cATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGc          >  1:3465891/1‑100 (MQ=255)
                                                                                         aTAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCt         >  2:328090/1‑100 (MQ=255)
                                                                                          tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTa        <  2:1748765/100‑1 (MQ=255)
                                                                                          tAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTa        <  1:1990263/100‑1 (MQ=255)
                                                                                                cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc  <  1:2832948/100‑1 (MQ=255)
                                                                                                cATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc  >  2:2889235/1‑100 (MQ=255)
                                                                                                |                                                                                                   
CAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACC  >  NC_000913/383801‑383996

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