Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,054,230 C→G 6.1% H254D (CAC→GAC)  PP_3573 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,054,2300CG 6.1% 287.0 / 3.9 91H254D (CAC→GAC) PP_3573monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (42/43);  new base G (3/3);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGTGCCTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGG  >  NC_002947/4054081‑4054378
                                                                                                                                                     |                                                                                                                                                     
cgcgTGCCTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGccc                                                                                                                                                       <  1:1736757/150‑1 (MQ=255)
    tGCCTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCccacca                                                                                                                                                   <  2:571914/150‑1 (MQ=255)
      ccgaCCGCAGCCGCGACTTCTGCTGTTGGCTGGGCGTGCTGGTCTTGTGGTACGCGGAAATTTCCAAGCCCTGTCGCGAGCATGTTACCATCGCGGTCAGTGGCGCCCTTGGTGGCCACACCGTGGTCTTCCTTGCTCTCGCCCTCCAgg                                                                                                                                                 <  2:390841/147‑1 (MQ=255)
       cTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAggg                                                                                                                                                >  1:1694873/1‑150 (MQ=255)
       cTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAggg                                                                                                                                                >  2:624857/1‑150 (MQ=255)
       cTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAg                                                                                                                                                  >  2:1521870/1‑148 (MQ=255)
       cTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAg                                                                                                                                                  <  1:1521870/148‑1 (MQ=255)
        tACNGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑c                                                                                                                                              >  2:1198544/1‑150 (MQ=255)
        tACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑c                                                                                                                                              >  1:770981/1‑150 (MQ=255)
          ccGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CAt                                                                                                                                            <  1:87471/150‑1 (MQ=255)
          ccGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CAt                                                                                                                                            >  1:1458662/1‑150 (MQ=255)
             cAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGAc                                                                                                                                         >  2:809273/1‑150 (MQ=255)
             cAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGCGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGAc                                                                                                                                         >  1:1282056/1‑150 (MQ=255)
                ccGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCt                                                                                                                                      >  2:1151028/1‑150 (MQ=255)
                    gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCACTCGCCCACCAGGGTCATG‑CCCTGGtg                                                                                                                                  <  1:731226/150‑2 (MQ=255)
                    gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTCCCCTCGCCCACCAGGG‑CCTGACCCTGGTc                                                                                                                                  >  2:586291/1‑150 (MQ=255)
                     aCTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCg                                                                                                                                 <  1:1822474/150‑1 (MQ=255)
                       ttCTGCTGGTGGCTGGGCGTGCGGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGc                                                                                                                               <  2:1541137/150‑1 (MQ=255)
                       ttCTGCTGGGGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCGAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATG‑CCCTGGTGGGc                                                                                                                               <  2:149833/150‑1 (MQ=255)
                        tCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGcc                                                                                                                              <  2:1116431/150‑1 (MQ=255)
                         ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCt                                                                                                                             <  2:902315/150‑1 (MQ=255)
                              ggtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATG‑CCCTGGTGGGcgagggca                                                                                                                        <  1:907232/150‑8 (MQ=255)
                               gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                          <  1:196604/147‑1 (MQ=255)
                               gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                          <  1:196611/147‑1 (MQ=255)
                               gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                          >  2:196611/1‑147 (MQ=255)
                               gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                          >  2:196604/1‑147 (MQ=255)
                                  gCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTc                                                                                                                    <  2:1833111/150‑1 (MQ=255)
                                   cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATG‑CCCTGGTGGGcgagggc                                                                                                                         <  1:787013/144‑7 (MQ=255)
                                   cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATG‑CCCTGGTGGGcgagggc                                                                                                                         >  2:787013/1‑138 (MQ=255)
                                      ggCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                <  1:1906727/150‑1 (MQ=255)
                                       gCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCCGACCCAGTAGTTCg                                                                                                               <  2:1268375/150‑1 (MQ=255)
                                         gTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGaa                                                                                                             <  2:1723966/150‑1 (MQ=255)
                                             tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                         <  2:38207/150‑1 (MQ=255)
                                             tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                         <  2:1030084/150‑1 (MQ=255)
                                                  gAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTc                                                                                                    >  1:1580811/1‑150 (MQ=255)
                                                     tGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgcg                                                                                                 >  1:278090/1‑150 (MQ=255)
                                                           gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                           >  1:1824289/1‑150 (MQ=255)
                                                             ggAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                         <  1:1149507/150‑1 (MQ=255)
                                                              gAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAg                                                                                        <  1:722760/150‑1 (MQ=255)
                                                                      cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCt                                                                                >  1:290086/1‑150 (MQ=255)
                                                                      cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCt                                                                                >  1:1731122/1‑150 (MQ=255)
                                                                       aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                               <  1:698925/150‑1 (MQ=255)
                                                                               tCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGaa                                                                       <  1:749558/150‑1 (MQ=255)
                                                                                 gcgAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACa                                                                     <  1:809273/150‑1 (MQ=255)
                                                                                   gAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATc                                                                   >  1:935108/1‑150 (MQ=255)
                                                                                    aGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCa                                                                  >  2:1110050/1‑150 (MQ=255)
                                                                                        tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGtt                                                                                                                 >  1:1605842/1‑99 (MQ=255)
                                                                                        tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGtt                                                                                                                 <  2:1605842/99‑1 (MQ=255)
                                                                                        tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                              >  2:896088/1‑150 (MQ=255)
                                                                                         gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCt                                                                                                                                      <  2:940556/77‑1 (MQ=255)
                                                                                         gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCt                                                                                                                                      >  1:940556/1‑77 (MQ=255)
                                                                                          tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCCCACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                          <  2:1163094/88‑1 (MQ=255)
                                                                                          tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGAc                                                                                                                          >  1:1163094/1‑88 (MQ=255)
                                                                                           gACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                >  1:1412414/1‑97 (MQ=255)
                                                                                           gACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                <  2:1412414/97‑1 (MQ=255)
                                                                                            aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                          <  1:1582077/150‑1 (MQ=255)
                                                                                              cATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcga                                                        >  1:686751/1‑150 (MQ=255)
                                                                                               aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgac                                                       >  1:1275384/1‑150 (MQ=255)
                                                                                               aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGAGGcgac                                                       <  2:935108/150‑1 (MQ=255)
                                                                                                tCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                          >  1:540133/1‑146 (MQ=255)
                                                                                                tCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                          <  2:540133/146‑1 (MQ=255)
                                                                                                 cgcgGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                     <  1:1151028/150‑1 (MQ=255)
                                                                                                    ggTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACaa                                                  <  2:1580811/150‑1 (MQ=255)
                                                                                                               cccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGGCGGTCGAAAACCGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGccc                                       >  1:1570429/1‑150 (MQ=255)
                                                                                                                ccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                              >  2:940758/1‑142 (MQ=255)
                                                                                                                ccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                              <  1:940758/142‑1 (MQ=255)
                                                                                                                  gtggtgGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgc                                    <  1:624857/150‑1 (MQ=255)
                                                                                                                   tggtggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                   <  1:1746972/150‑1 (MQ=255)
                                                                                                                       ggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGAt                               <  2:1458662/150‑1 (MQ=255)
                                                                                                                        gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGcaca                                                                     <  2:506102/111‑1 (MQ=255)
                                                                                                                        gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACa                                                                     >  1:506102/1‑111 (MQ=255)
                                                                                                                        gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATg                              <  1:1198544/150‑1 (MQ=255)
                                                                                                                          cacaCCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCg                            <  2:1694873/150‑1 (MQ=255)
                                                                                                                           acacCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                          <  2:249512/135‑1 (MQ=255)
                                                                                                                           acacCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                          >  1:249512/1‑135 (MQ=255)
                                                                                                                               cGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGa                       >  1:1793425/1‑150 (MQ=255)
                                                                                                                                gTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTAGAGCATATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGAt                      >  1:1057679/1‑150 (MQ=255)
                                                                                                                                      gagggggCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATAAATCGAGGCGACGAAAACTAACTGGCCCTGCTGGATGCGGCCGATGCCTAc                >  1:544188/7‑150 (MQ=255)
                                                                                                                                         cGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc             <  2:1814489/150‑1 (MQ=255)
                                                                                                                                            gCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGaa          <  2:1687825/150‑1 (MQ=255)
                                                                                                                                             cccTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAAc         <  2:325896/150‑1 (MQ=255)
                                                                                                                                              ccTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCt                                            >  1:474047/1‑114 (MQ=255)
                                                                                                                                              ccTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCt                                            <  2:474047/114‑1 (MQ=255)
                                                                                                                                              ccTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTAGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACg        >  1:724831/1‑150 (MQ=255)
                                                                                                                                               cTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTGGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc       >  1:1095333/1‑150 (MQ=255)
                                                                                                                                               cTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc       >  1:5586/1‑150 (MQ=255)
                                                                                                                                               cTCGCCCACCACAG‑CACGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGGCAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc       >  2:1262170/1‑150 (MQ=255)
                                                                                                                                                tCGCCCACCAGGG‑CATGACCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGAGAACCTGGTCGACAACATCACTCGAGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa      >  1:516914/1‑150 (MQ=255)
                                                                                                                                                tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa      >  2:556343/1‑150 (MQ=255)
                                                                                                                                                tCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa      >  1:1588886/1‑150 (MQ=255)
                                                                                                                                                  gCCGACCAGGGTCATG‑CCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                              >  2:1833088/1‑108 (MQ=38)
                                                                                                                                                  gCCGACCAGGGTCATG‑CCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc    >  1:1252107/1‑150 (MQ=255)
                                                                                                                                                  gCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTAGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc    <  2:724831/150‑1 (MQ=255)
                                                                                                                                                    ccaccaGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACgg  >  2:1302620/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGTGCCTACCGCAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGG‑CATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGG  >  NC_002947/4054081‑4054378

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: ‑