Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 535,821 G→T 7.4% E85* (GAA→TAA) ‡ rplA → 50S ribosomal protein L1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947535,8210GT 7.4% 205.3 / 5.5 68E85* (GAA→TAA) ‡rplA50S ribosomal protein L1
Reads supporting (aligned to +/- strand):  ref base G (26/36);  new base T (2/3);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.24e-01

CAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGG  >  NC_002947/535676‑535962
                                                                                                                                                  |                                                                                                                                             
cAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAgc                                                                                                                                           >  1:684740/1‑150 (MQ=255)
     tcgtcgAGTCCTCCGACATCGCCGTTCCCCTCGGTGTTGACCCGCGTAAATCCGCCCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGctc                                                                                                                                      <  2:654945/150‑1 (MQ=255)
     tcgtcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACTGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCATCCAGGGTCCAGCAGCTGAAGCTGctc                                                                                                                                      <  1:332883/150‑1 (MQ=255)
      cgtcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGctct                                                                                                                                     <  1:561564/150‑1 (MQ=255)
       gtcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTTCAgctgctgga                                                                                                                                    <  2:735720/150‑4 (MQ=255)
       gtcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTg                                                                                                                                    >  2:857573/1‑150 (MQ=255)
        tcgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCCGGTCCTGAACA‑AGCTGCTGGTGCTGCCACACGCCATTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTgg                                                                                                                                   <  1:678655/150‑1 (MQ=255)
         cgAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGc                                                                                                                                  <  2:684740/150‑1 (MQ=255)
            gTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGAGTCCAGCAGCTGAAGCTGCTCTGGCTGc                                                                                                                               >  2:307989/1‑150 (MQ=255)
                       tCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATc                                                                                                                    >  2:1288545/1‑150 (MQ=255)
                            gTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTa                                                                                                               <  1:307989/150‑1 (MQ=255)
                             ttAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAg                                                                                                              <  1:1290853/150‑1 (MQ=255)
                             ttAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAg                                                                                                              <  2:419113/150‑1 (MQ=255)
                             ttAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAg                                                                                                              <  2:959463/150‑1 (MQ=255)
                                   tCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc                                                                                                                         <  2:539784/133‑1 (MQ=255)
                                   tCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc                                                                                                                         >  1:539784/1‑133 (MQ=255)
                                    cGGTGTAGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGa                                                                                                       <  2:868727/150‑1 (MQ=255)
                                        gtTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGGACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATgg                                                                                                        >  2:919175/1‑145 (MQ=255)
                                        gtTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATgg                                                                                                        <  1:919175/145‑1 (MQ=255)
                                           gACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTg                                                                                                >  2:254152/1‑150 (MQ=255)
                                              ccGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACAAAGGGTCAAACAGATAATGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCt                                                                                             >  1:1387104/1‑150 (MQ=255)
                                                      aTCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTAt                                                                                                          >  1:1364188/1‑129 (MQ=255)
                                                      aTCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTAt                                                                                                          <  2:1364188/129‑1 (MQ=255)
                                                      aTCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGGATGGACGATCTGGCTGCCGAAAt                                                                                     <  2:590451/150‑1 (MQ=255)
                                                       tCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATg                                                                                    >  2:1273721/1‑150 (MQ=255)
                                                            ccAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAgg                                                                               >  2:561191/1‑150 (MQ=255)
                                                             cAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAggc                                                                              <  1:1273721/150‑1 (MQ=255)
                                                                    tacgtaGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGAcc                                                                       >  1:500752/1‑150 (MQ=255)
                                                                      cgtaGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTg                                                                     <  1:254152/150‑1 (MQ=255)
                                                                       gtaGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGa                                                                    >  1:179355/1‑150 (MQ=255)
                                                                        taGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa                                                                   >  2:814901/1‑150 (MQ=255)
                                                                            gcTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt                                                               >  2:567911/1‑150 (MQ=255)
                                                                            gcTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt                                                               <  1:1288545/150‑1 (MQ=255)
                                                                            gcTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt                                                               <  2:1310256/150‑1 (MQ=255)
                                                                            gcTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTAt                                                               >  2:983453/1‑150 (MQ=255)
                                                                             cTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATg                                                              >  1:1141627/1‑150 (MQ=255)
                                                                              tAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGaa                                                                                                                                             <  1:338888/70‑1 (MQ=255)
                                                                              tAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGaa                                                                                                                                             >  2:338888/1‑70 (MQ=255)
                                                                                        gcCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTat                                                    >  1:282968/1‑150 (MQ=255)
                                                                                         ccACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatc                                                   <  1:561191/150‑1 (MQ=255)
                                                                                         ccACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTatc                                                   <  2:179355/150‑1 (MQ=255)
                                                                                             acGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc                                                                                                                         >  2:1321938/1‑76 (MQ=255)
                                                                                             acGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGgcgc                                                                                                                         <  1:1321938/76‑1 (MQ=255)
                                                                                                                 cgcgTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTCATCGCATCGCCTGATGCCATGCGCGtt                           <  2:1141627/150‑1 (MQ=255)
                                                                                                                  gcgTTGCTGTCTTCACCCAGGGTACAGCAGCTGAAGCTGCCCTGGCTGCCGGCGCTGATCGCGTAGGTATGGCCGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATAACGTCGTTCTCGCATCGCCTGATGCCATGCGCGTTg                          >  2:551675/1‑150 (MQ=255)
                                                                                                                   cgTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGt                         <  1:1066720/150‑1 (MQ=255)
                                                                                                                    gTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTa                        <  1:857573/150‑1 (MQ=255)
                                                                                                                     ttGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAg                       <  2:975215/150‑1 (MQ=255)
                                                                                                                       gCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGTCCTATACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTATGt                     <  2:1267276/150‑1 (MQ=255)
                                                                                                                             ttCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTg               >  1:563253/1‑150 (MQ=255)
                                                                                                                                aCCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa                                                                   <  2:321395/95‑1 (MQ=255)
                                                                                                                                aCCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGaa                                                                   >  1:321395/1‑95 (MQ=255)
                                                                                                                                aCCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt            >  2:2979/1‑150 (MQ=255)
                                                                                                                                aCCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGt            <  2:1408365/150‑1 (MQ=255)
                                                                                                                                 cccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTat                                                    >  2:846914/1‑109 (MQ=255)
                                                                                                                                 cccAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGAGCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTat                                                    <  1:846914/109‑1 (MQ=255)
                                                                                                                                        tCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATNGCATCGNCTGATGCCATGCNCGTTGTAGGTCAGCTGGGTCAGGTTCt    <  1:983453/150‑1 (MQ=255)
                                                                                                                                        cagagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATgc                                <  2:4839/119‑1 (MQ=255)
                                                                                                                                        cagagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATgc                                >  1:4839/4‑122 (MQ=255)
                                                                                                                                         agagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTg   <  1:1265937/148‑1 (MQ=255)
                                                                                                                                         agagcagcTTCAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTg   >  2:1265937/3‑150 (MQ=255)
                                                                                                                                          cagcagCTGAAGCTGCTTTGGCTGCCGGCGCCGATCACGTAGGTTTGCACGATCTTGCTGCCGAAATGAAAGGCGGCGACCTGAGCTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg  <  2:447733/150‑1 (MQ=255)
                                                                                                                                          cagcagCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAg         <  1:1205686/143‑1 (MQ=255)
                                                                                                                                          cagcagCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAg         >  2:1205686/1‑143 (MQ=255)
                                                                                                                                          cagcagCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg  <  1:814901/150‑1 (MQ=255)
                                                                                                                                          aagcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTgg  <  1:567911/149‑1 (MQ=255)
                                                                                                                                           agcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATgcgc                              <  1:285965/121‑1 (MQ=255)
                                                                                                                                           agcagcTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATgcgc                              >  2:285965/1‑121 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CAAGTTCGTCGAGTCCTACGACATCGCCGTTAACCTCGGTGTTGACCCGCGTAAATCCGACCAGGTCGTACGTAGCGCTAC‑CGTGCTGCCACACGGCACTGGCAAGACCGTTCGCGTTGCTGTCTTCACCCAGGGTCCAGCAGCTGAAGCTGCTCTGGCTGCCGGCGCTGATCGCGTAGGTATGGACGATCTGGCTGCCGAAATGAAAGGCGGCGACCTGAACTATGACGTCGTTATCGCATCGCCTGATGCCATGCGCGTTGTAGGTCAGCTGGGTCAGGTTCTGG  >  NC_002947/535676‑535962

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