Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,491,144 G→C 9.6% S267T (AGT→ACT)  cysNC → sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,491,1440GC 9.6% 349.5 / 24.2 114S267T (AGT→ACT) cysNCsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (43/60);  new base C (5/6);  total (48/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 114 total aligned reads displayed.

GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAT  >  NC_002947/1490998‑1491291
                                                                                                                                                   |                                                                                                                                                    
gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                                   >  2:536493/1‑150 (MQ=255)
gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                                   >  1:1937477/1‑150 (MQ=255)
gCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGAGTGTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTgg                                                                                                                                                   >  1:1733061/1‑150 (MQ=255)
   gTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgt                                                                                                                                               <  1:297143/150‑1 (MQ=255)
    tGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtg                                                                                                                                              <  2:2309969/150‑1 (MQ=255)
    tGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCACGCCACTGGCgatg                                                                                                                                              <  2:1979851/150‑4 (MQ=255)
     ggTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtgg                                                                                                                                             <  1:536493/150‑1 (MQ=255)
      gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCgtggt                                                                                                                                            <  2:1473048/150‑1 (MQ=255)
      gTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCACGCCACTGGCgatggt                                                                                                                                            <  1:4034/150‑6 (MQ=255)
        aCGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGc                                                                                                                                          <  2:434172/150‑1 (MQ=255)
        aCGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTCCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCACGCCACTGGCGATGGTGc                                                                                                                                          <  1:907194/150‑1 (MQ=255)
          gCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGcac                                                                                                                                        >  1:776506/1‑150 (MQ=255)
            cGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACaa                                                                                                                                      >  2:407042/1‑150 (MQ=255)
             ggCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACaaa                                                                                                                                     <  1:571172/150‑1 (MQ=255)
              gCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                                    <  2:1157785/150‑1 (MQ=255)
              gCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                                    <  1:1856544/150‑1 (MQ=255)
              gCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                                    >  1:1809163/1‑150 (MQ=255)
               cccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAgg                                                                                                                                   <  2:1384002/150‑1 (MQ=255)
                ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGt                                                                                                                                  <  1:1466814/150‑1 (MQ=255)
                  gACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGa                                                                                                                                <  2:1073191/150‑1 (MQ=255)
                  gACGCCGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACGCCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCCACCTGAACTGCCGCGGCTTCGCGGGCATCATCGCCAGTGGCGTGGTGCACAAAGGTGa                                                                                                                                <  1:858684/150‑1 (MQ=255)
                   aCGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGAc                                                                                                                               <  1:2189770/150‑1 (MQ=255)
                     gCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                                    >  1:11004/1‑143 (MQ=255)
                     gCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACga                                                                                                                             <  1:2153192/150‑1 (MQ=255)
                            gAAATCCTCGAACCCGGCGAGGTTTCGGCTGACCGCAACTGCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCCGTGGCGTGGTGCACAAAGGTGACGAGATTGTc                                                                                                                      <  2:276160/150‑1 (MQ=255)
                            gAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTc                                                                                                                      <  2:1333654/150‑1 (MQ=255)
                                 ccTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgct                                                                                                                 <  1:1109321/150‑1 (MQ=255)
                                  cTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACTTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctg                                                                                                                >  2:1596834/1‑150 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGCCGGAATGGTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGGCCCGACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGAGGAGATTGTCGtgcggg                                                                                                               >  1:63402/1‑147 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGa                                                                                                                                <  1:441774/133‑1 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGa                                                                                                                                >  2:441774/1‑133 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                               >  2:5536/1‑150 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                               >  1:2222421/1‑150 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCCGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                               >  1:1695950/1‑150 (MQ=255)
                                   ccGAAACCTCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                                    >  2:612034/2‑129 (MQ=255)
                                         ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                         <  2:2062069/150‑1 (MQ=255)
                                                 ggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAagag                                                                                                  <  2:701142/150‑1 (MQ=255)
                                                        gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                           <  2:2222443/150‑1 (MQ=255)
                                                        gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                           >  2:458530/1‑150 (MQ=255)
                                                              cGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTTGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAg                                                                                     >  1:6548/1‑150 (MQ=255)
                                                              cGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAg                                                                                     >  2:1863305/1‑150 (MQ=255)
                                                               gCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACgag                                                                                                                            <  2:2047358/110‑1 (MQ=255)
                                                               gCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACgag                                                                                                                            >  1:2047358/1‑110 (MQ=255)
                                                               gCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACATCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGt                                                                                    >  2:326271/1‑150 (MQ=255)
                                                                    ttCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                               >  2:2437433/1‑150 (MQ=255)
                                                                          gACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCAcc                                                                         >  1:861559/1‑150 (MQ=255)
                                                                              tGCGTTTCGCGGTGCAGTCCGTCAACCGTCCCAACCTGAACTTCCGCGGCTCCGCGGGCACCATCCCCAGTGGCGTGGTGCGCAAGGGTGACGACATTGTCGTGCGGCCGTCGGGCAGGAGCAGCCGGGCCAAGTCCATCGTCAcccgcg                                                                     >  2:218637/1‑146 (MQ=255)
                                                                                 gTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATGGTCCTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                                  <  2:1796914/150‑1 (MQ=255)
                                                                                  tttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAgg                                                                 >  1:1582325/1‑150 (MQ=255)
                                                                                   ttCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGt                                                                <  2:1695950/150‑1 (MQ=255)
                                                                                    tCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTg                                                               <  2:2062726/150‑1 (MQ=255)
                                                                                           gCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGa                                                        <  2:1315467/150‑1 (MQ=255)
                                                                                               tACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAc                                                    >  2:361400/1‑150 (MQ=255)
                                                                                               tACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAc                                                    >  2:52036/1‑150 (MQ=255)
                                                                                                  gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCt                                                 <  2:1622967/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCCACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                                <  1:188731/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                                >  2:373187/1‑150 (MQ=255)
                                                                                                   tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                                <  1:361400/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                                <  2:1874832/150‑1 (MQ=255)
                                                                                                    cAACCGTCCCAACCTGAACTTCCGCAGCTCCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                               >  1:1245018/1‑150 (MQ=255)
                                                                                                     aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCCCCAGTGGCGTGGTGCCCAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAACCGGGTCAAGTCCATCGTCACCTACGAAAGGGAACTGGAAAACGCTGGc                                              >  2:1419941/1‑150 (MQ=255)
                                                                                                           cccAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggc                                        <  2:1163682/150‑1 (MQ=255)
                                                                                                            ccAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggcc                                       >  1:2058247/1‑150 (MQ=255)
                                                                                                             cAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggcca                                      >  2:26618/1‑150 (MQ=255)
                                                                                                              aaCCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggccag                                     <  1:52036/150‑1 (MQ=255)
                                                                                                               aCCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggccagg                                    <  1:1596834/150‑1 (MQ=255)
                                                                                                                ccTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCcaggccaggc                                   >  2:1749943/1‑150 (MQ=255)
                                                                                                                 cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCaggccaggcc                                  >  1:2361706/1‑150 (MQ=255)
                                                                                                                  tGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                                  <  1:1470992/117‑1 (MQ=255)
                                                                                                                  tGAACTTCCGCGGCTTCGCGGGACATATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCg                                 >  1:1519910/1‑150 (MQ=255)
                                                                                                                   gAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                                <  1:2434878/150‑1 (MQ=255)
                                                                                                                   gAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                                >  1:2189220/1‑150 (MQ=255)
                                                                                                                     aCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGCCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGa                              <  2:2141800/150‑1 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                               <  1:551977/95‑1 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                               >  2:551977/1‑95 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaa                                                      >  2:1508910/1‑120 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaa                                                      <  1:1508910/120‑1 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCctgggccagcg                                   <  1:332451/139‑3 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGa                        <  1:326271/150‑1 (MQ=255)
                                                                                                                            cgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGAc                       <  2:729895/150‑1 (MQ=255)
                                                                                                                               cTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCAt                    <  2:1582325/150‑1 (MQ=255)
                                                                                                                                 ctttgtgCACCACGCCACTGGCGATGGTGCCCAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                                >  2:1664627/7‑136 (MQ=255)
                                                                                                                                 ctttgtgCACCACGCCACTGGCGATGGTGCCCAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                                <  1:1664627/130‑1 (MQ=255)
                                                                                                                                   gcgGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggccag                                     <  1:1318264/129‑1 (MQ=255)
                                                                                                                                   gcgGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCccaggccag                                     >  2:1318264/1‑129 (MQ=255)
                                                                                                                                      tgCACCACGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                            >  1:2148862/2‑103 (MQ=255)
                                                                                                                                      tgCACCACGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                            <  2:2148862/102‑1 (MQ=255)
                                                                                                                                      tgCACCACGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc             >  2:620925/2‑150 (MQ=255)
                                                                                                                                        cACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACga           <  2:1733061/150‑1 (MQ=255)
                                                                                                                                        cACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACga           <  1:1736504/150‑1 (MQ=255)
                                                                                                                                         aCCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACgag          >  1:592244/1‑150 (MQ=255)
                                                                                                                                          accacGCCACTGGCGATGGTGCCCAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCAt                    >  2:1538436/5‑140 (MQ=255)
                                                                                                                                          accacGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACgag          >  2:2388304/5‑150 (MQ=255)
                                                                                                                                           cATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCaa                                                                                      >  1:1756268/1‑72 (MQ=255)
                                                                                                                                           cATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCaa                                                                                      <  2:1756268/72‑1 (MQ=255)
                                                                                                                                             tCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAc                                                    >  1:529760/1‑104 (MQ=255)
                                                                                                                                              cGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGa     <  1:5536/150‑1 (MQ=255)
                                                                                                                                               gCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGAc    <  1:407042/150‑1 (MQ=255)
                                                                                                                                                ccAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCCCATACGAAGTTGATATTGTGATCGCGTTGTCAATCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACa   <  1:1863305/150‑1 (MQ=255)
                                                                                                                                                 cAGGGGCGTGGTGCACAAAGGTGACGAGATTGTCGAGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAt  >  2:1276817/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAGATCGACAT  >  NC_002947/1490998‑1491291

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: