Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,054,250 C→G 6.8% V260V (GTC→GTG PP_3573 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,054,2500CG 6.8% 305.3 / 8.4 96V260V (GTC→GTGPP_3573monooxygenase
Reads supporting (aligned to +/- strand):  ref base C (43/46);  new base G (3/4);  total (46/50)
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

GACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCGGAAGA  >  NC_002947/4054101‑4054396
                                                                                                                                                      |                                                                                                                                                  
gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTCCCCTCGCCCACCAGGGCCTGACCCTGGTc                                                                                                                                                    >  2:586291/1‑150 (MQ=255)
 aCTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCg                                                                                                                                                   <  1:1822474/150‑1 (MQ=255)
   ttCTGCTGGTGGCTGGGCGTGCGGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGc                                                                                                                                                 <  2:1541137/150‑1 (MQ=255)
   ttCTGCTGGGGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCGAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATGCCCTGGTGGGc                                                                                                                                                 <  2:149833/150‑1 (MQ=255)
    tCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGcc                                                                                                                                                <  2:1116431/150‑1 (MQ=255)
     ctgctgGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCt                                                                                                                                               <  2:902315/150‑1 (MQ=255)
          ggtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATGCCCTGGTGGGcgagggca                                                                                                                                          <  1:907232/150‑8 (MQ=255)
           gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                            >  2:196604/1‑147 (MQ=255)
           gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                            >  2:196611/1‑147 (MQ=255)
           gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                            <  1:196604/147‑1 (MQ=255)
           gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                            <  1:196611/147‑1 (MQ=255)
              gCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTc                                                                                                                                      <  2:1833111/150‑1 (MQ=255)
               cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATGCCCTGGTGGGcgagggc                                                                                                                                           <  1:787013/144‑7 (MQ=255)
               cTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCGACCAGGGTCATGCCCTGGTGGGcgagggc                                                                                                                                           >  2:787013/1‑138 (MQ=255)
                  ggCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                                  <  1:1906727/150‑1 (MQ=255)
                   gCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCCGACCCAGTAGTTCg                                                                                                                                 <  2:1268375/150‑1 (MQ=255)
                     gTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGaa                                                                                                                               <  2:1723966/150‑1 (MQ=255)
                         tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                                           <  2:38207/150‑1 (MQ=255)
                         tGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACg                                                                                                                           <  2:1030084/150‑1 (MQ=255)
                              gAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTc                                                                                                                      >  1:1580811/1‑150 (MQ=255)
                                 tGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTcgcg                                                                                                                   >  1:278090/1‑150 (MQ=255)
                                       gcgGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTc                                                                                                             >  1:1824289/1‑150 (MQ=255)
                                         ggAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCa                                                                                                           <  1:1149507/150‑1 (MQ=255)
                                          gAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAg                                                                                                          <  1:722760/150‑1 (MQ=255)
                                                  cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCt                                                                                                  >  1:1731122/1‑150 (MQ=255)
                                                  cAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCt                                                                                                  >  1:290086/1‑150 (MQ=255)
                                                   aaGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTg                                                                                                 <  1:698925/150‑1 (MQ=255)
                                                           tCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGaa                                                                                         <  1:749558/150‑1 (MQ=255)
                                                             gcgAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACa                                                                                       <  1:809273/150‑1 (MQ=255)
                                                               gAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATc                                                                                     >  1:935108/1‑150 (MQ=255)
                                                                aGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCa                                                                                    >  2:1110050/1‑150 (MQ=255)
                                                                    tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGtt                                                                                                                                   >  1:1605842/1‑99 (MQ=255)
                                                                    tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGtt                                                                                                                                   <  2:1605842/99‑1 (MQ=255)
                                                                    tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                                >  2:896088/1‑150 (MQ=255)
                                                                      tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCCCACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                            <  2:1163094/88‑1 (MQ=255)
                                                                      tgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                            >  1:1163094/1‑88 (MQ=255)
                                                                       gACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                                  <  2:1412414/97‑1 (MQ=255)
                                                                       gACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTc                                                                                                                                  >  1:1412414/1‑97 (MQ=255)
                                                                        aCCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                                            <  1:1582077/150‑1 (MQ=255)
                                                                          cATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcga                                                                          >  1:686751/1‑150 (MQ=255)
                                                                           aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgac                                                                         >  1:1275384/1‑150 (MQ=255)
                                                                           aTCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGAGGcgac                                                                         <  2:935108/150‑1 (MQ=255)
                                                                            tCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                                            >  1:540133/1‑146 (MQ=255)
                                                                            tCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                                            <  2:540133/146‑1 (MQ=255)
                                                                             cgcgGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGcgacga                                                                       <  1:1151028/150‑1 (MQ=255)
                                                                                ggTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACaa                                                                    <  2:1580811/150‑1 (MQ=255)
                                                                                           cccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGGCGGTCGAAAACCGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGccc                                                         >  1:1570429/1‑150 (MQ=255)
                                                                                            ccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                                <  1:940758/142‑1 (MQ=255)
                                                                                            ccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                                >  2:940758/1‑142 (MQ=255)
                                                                                              gtggtgGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgc                                                      <  1:624857/150‑1 (MQ=255)
                                                                                               tggtggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                                     <  1:1746972/150‑1 (MQ=255)
                                                                                                   ggCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGAt                                                 <  2:1458662/150‑1 (MQ=255)
                                                                                                    gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGcaca                                                                                       <  2:506102/111‑1 (MQ=255)
                                                                                                    gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACa                                                                                       >  1:506102/1‑111 (MQ=255)
                                                                                                    gCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATg                                                <  1:1198544/150‑1 (MQ=255)
                                                                                                      cacaCCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCg                                              <  2:1694873/150‑1 (MQ=255)
                                                                                                       acacCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                                            <  2:249512/135‑1 (MQ=255)
                                                                                                       acacCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTgg                                                            >  1:249512/1‑135 (MQ=255)
                                                                                                           cGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGa                                         >  1:1793425/1‑150 (MQ=255)
                                                                                                            gTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTAGAGCATATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGAt                                        >  1:1057679/1‑150 (MQ=255)
                                                                                                                  gagggggCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATAAATCGAGGCGACGAAAACTAACTGGCCCTGCTGGATGCGGCCGATGCCTAc                                  >  1:544188/7‑150 (MQ=255)
                                                                                                                     cGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc                               <  2:1814489/150‑1 (MQ=255)
                                                                                                                        gCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGaa                            <  2:1687825/150‑1 (MQ=255)
                                                                                                                         cccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAAc                           <  2:325896/150‑1 (MQ=255)
                                                                                                                          ccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCt                                                              >  1:474047/1‑114 (MQ=255)
                                                                                                                          ccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCt                                                              <  2:474047/114‑1 (MQ=255)
                                                                                                                          ccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTAGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACg                          >  1:724831/1‑150 (MQ=255)
                                                                                                                           cTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTGGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc                         >  1:1095333/1‑150 (MQ=255)
                                                                                                                           cTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc                         >  1:5586/1‑150 (MQ=255)
                                                                                                                           cTCGCCCACCACAGCACGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGGCAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGc                         >  2:1262170/1‑150 (MQ=255)
                                                                                                                            tCGCCCACCAGGGCATGACCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGAGAACCTGGTCGACAACATCACTCGAGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa                        >  1:516914/1‑150 (MQ=255)
                                                                                                                            tCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa                        >  1:1588886/1‑150 (MQ=255)
                                                                                                                            tCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa                        >  2:556343/1‑150 (MQ=255)
                                                                                                                              gccgaccaGGGTCATGCCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                                <  1:1833088/104‑1 (MQ=38)
                                                                                                                              gCCGACCAGGGTCATGCCCTGGTGGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                                >  2:1833088/1‑108 (MQ=38)
                                                                                                                              gCCGACCAGGGTCATGCCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc                      >  1:1252107/1‑150 (MQ=255)
                                                                                                                              gCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTAGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAAc                      <  2:724831/150‑1 (MQ=255)
                                                                                                                                ccaccaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACgg                    >  2:1302620/1‑150 (MQ=255)
                                                                                                                                  accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGAc                                                                                                       <  2:477491/65‑1 (MQ=255)
                                                                                                                                  accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGAc                                                                                                       >  1:477491/1‑65 (MQ=255)
                                                                                                                                  accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                                            <  2:1122248/92‑1 (MQ=255)
                                                                                                                                  accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                                            >  1:1122248/1‑92 (MQ=255)
                                                                                                                                  accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                                    <  2:1691695/132‑1 (MQ=255)
                                                                                                                                  accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                                    >  1:1691695/1‑132 (MQ=255)
                                                                                                                                       ggCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                                     <  1:1460534/110‑1 (MQ=255)
                                                                                                                                       ggCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                                     >  2:1460534/1‑110 (MQ=255)
                                                                                                                                       ggCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGAc             <  2:770981/150‑1 (MQ=255)
                                                                                                                                              aCCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCg       >  1:1495157/1‑150 (MQ=255)
                                                                                                                                               cccTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                                     <  1:1746435/103‑1 (MQ=255)
                                                                                                                                               cccTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCctgct                                                     >  2:1746435/1‑103 (MQ=255)
                                                                                                                                               cccTGGTCGGCCTGACCCAGTCGTTCGAAAACGGAGTCGCGCCCGTCCAGGACAACCTGGTCGAGAACATCAATGGCGGCGACGACAACTACGTGGACCTACTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCGGGACCTGcccg      >  1:275991/1‑148 (MQ=255)
                                                                                                                                                ccTGGTCGGCCTGACCCATTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCCACGGCCTGGACCTGCCggg     >  2:556811/1‑149 (MQ=255)
                                                                                                                                                ccTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc                               <  2:235528/124‑1 (MQ=255)
                                                                                                                                                ccTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATc                               >  1:235528/1‑124 (MQ=255)
                                                                                                                                                ccTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCGGa     >  1:1540485/1‑150 (MQ=255)
                                                                                                                                                  tGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCGGAag   <  2:1588886/150‑1 (MQ=255)
                                                                                                                                                   ggTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCATACATCGAACGCAAGGGCCTGGACCTGCCCGAaga  >  1:587905/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCGGAAGA  >  NC_002947/4054101‑4054396

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: