Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,491,138 Δ1 bp 7.7% coding (794/1902 nt) 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,1380T. 7.7% 296.4 / 20.3 104coding (794/1902 nt)cysNCsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase
Reads supporting (aligned to +/- strand):  ref base T (41/55);  new base . (3/5);  total (44/60)
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.89e-02

Only 100 of 106 total aligned reads displayed.

GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAG  >  NC_002947/1490998‑1491283
                                                                                                                                             |                                                                                                                                                  
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)
        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:1809163/1‑150 (MQ=255)
               cccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAgg                                                                                                                           <  2:1384002/150‑1 (MQ=255)
                ccGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGt                                                                                                                          <  2:1670069/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:1756510/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                                                                                                                        >  2:441774/1‑133 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGa                                                                                                                        <  1:441774/133‑1 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                       >  1:2222421/1‑150 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                       >  2:5536/1‑150 (MQ=255)
                                   tCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCCGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGtgctgc                                                                                                       >  1:1695950/1‑150 (MQ=255)
                                   ccGAAACCTCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                            <  1:612034/128‑1 (MQ=255)
                                   ccGAAACCTCGACGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAg                                                                                                                            >  2:612034/2‑129 (MQ=255)
                                         ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 <  1:2095176/150‑1 (MQ=255)
                                         ccGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCgg                                                                                                 <  2:2062069/150‑1 (MQ=255)
                                                 ggTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAagag                                                                                          <  2:701142/150‑1 (MQ=255)
                                                        gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                   >  2:458530/1‑150 (MQ=255)
                                                        gCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCgg                                                                                   <  2:2222443/150‑1 (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)
                                                                    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:52036/1‑150 (MQ=255)
                                                                                               tACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAc                                            >  2:361400/1‑150 (MQ=255)
                                                                                                  gTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCt                                         <  2:1622967/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCCACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                        <  1:188731/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                        <  2:1874832/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                        <  1:361400/150‑1 (MQ=255)
                                                                                                   tCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTg                                        >  2:373187/1‑150 (MQ=255)
                                                                                                    cAACCGTCCCAACCTGAACTTCCGCAGCTCCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTgg                                       >  1:1245018/1‑150 (MQ=255)
                                                                                                     aaCCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCCCCAGTGGCGTGGTGCCCAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAACCGGGTCAAGTCCATCGTCACCTACGAAAGGGAACTGGAAAACGCTGGc                                      >  2:1419941/1‑150 (MQ=255)
                                                                                                       ccGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGGGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGccc                                    <  1:471011/150‑1 (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)
                                                                                                                 cTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCaggccaggcc                          >  1:2361706/1‑150 (MQ=255)
                                                                                                                  tGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                          <  1:1470992/117‑1 (MQ=255)
                                                                                                                  tGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAg                                                          >  2:1470992/1‑117 (MQ=255)
                                                                                                                  tGAACTTCCGCGGCTTCGCGGGACATATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCg                         >  1:1519910/1‑150 (MQ=255)
                                                                                                                   gAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                        >  1:2189220/1‑150 (MQ=255)
                                                                                                                   gAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGt                        <  1:2434878/150‑1 (MQ=255)
                                                                                                                     aCTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGCCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGa                      <  2:2141800/150‑1 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                       >  2:551977/1‑95 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATc                                                                       <  1:551977/95‑1 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaa                                              <  1:1508910/120‑1 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGaaa                                              >  2:1508910/1‑120 (MQ=255)
                                                                                                                           gcgGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCctgggccagcg                           >  2:332451/1‑137 (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)
                                                                                                                                     gggCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGa      <  2:1937477/150‑1 (MQ=255)
                                                                                                                                      tgCACCACGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                    >  1:2148862/2‑103 (MQ=255)
                                                                                                                                      tgCACCACGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAAc                                                    <  2:2148862/102‑1 (MQ=255)
                                                                                                                                      tgCACCACGCCACTGGCGATGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGAc     >  2:620925/2‑150 (MQ=255)
                                                                                                                                        cACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACga   <  1:1736504/150‑1 (MQ=255)
                                                                                                                                        cACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACga   <  2:1733061/150‑1 (MQ=255)
                                                                                                                                         aCCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACgag  >  1:592244/1‑150 (MQ=255)
                                                                                                                                         aCCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACgag  <  1:2437433/150‑1 (MQ=255)
                                                                                                                                           cATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCaa                                                                              >  1:1756268/1‑72 (MQ=255)
                                                                                                                                           cATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCaa                                                                              <  2:1756268/72‑1 (MQ=255)
                                                                                                                                             tCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAc                                            <  2:529760/104‑1 (MQ=255)
                                                                                                                                             tCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAAc                                            >  1:529760/1‑104 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
GCCGTGGTACGCCGGCCCGACGCTGATGGAAATCCTCGAAACCGTCGAGGTTTCGGCTGACCGCAACTTCACCGACCTGCGTTTCCCGGTGCAGTACGTCAACCGTCCCAACCTGAACTTCCGCGGCTTCGCGGGCACCATCGCCAGTGGCGTGGTGCACAAAGGTGACGAGATTGTCGTGCTGCCGTCGGGCAAGAGCAGCCGGGTCAAGTCCATCGTCACCTACGAAGGTGAACTGGAAAACGCTGGCCCAGGCCAGGCCGTGACCCTGACCATGGAAGACGAG  >  NC_002947/1490998‑1491283

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: