Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,054,242 Δ1 bp 8.6% coding (772/1320 nt) PP_3573 → monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,054,2420A. 8.6% 263.3 / 20.5 93coding (772/1320 nt)PP_3573monooxygenase
Reads supporting (aligned to +/- strand):  ref base A (41/44);  new base . (3/5);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.82e-01

CAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCG  >  NC_002947/4054094‑4054391
                                                                                                                                                     |                                                                                                                                                     
cAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGAc                                                                                                                                                      >  2:809273/1‑150 (MQ=255)
cAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGCGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGAc                                                                                                                                                      >  1:1282056/1‑150 (MQ=255)
   ccGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCt                                                                                                                                                   >  2:1151028/1‑150 (MQ=255)
       gACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCACTCGCCCACCAGGGTCATGCCCTGGtg                                                                                                                                               <  1:731226/150‑2 (MQ=255)
       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                                                                                                                                       <  1:196611/147‑1 (MQ=255)
                  gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                       <  1:196604/147‑1 (MQ=255)
                  gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                       >  2:196604/1‑147 (MQ=255)
                  gtggCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGAc                                                                                                                                       >  2:196611/1‑147 (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                                                                                                                              <  2:1605842/99‑1 (MQ=255)
                                                                           tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGtt                                                                                                                              >  1:1605842/1‑99 (MQ=255)
                                                                           tgtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATcg                                                                           >  2:896088/1‑150 (MQ=255)
                                                                            gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCt                                                                                                                                                   >  1:940556/1‑77 (MQ=255)
                                                                            gtgACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCt                                                                                                                                                   <  2:940556/77‑1 (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                                                           >  2:940758/1‑142 (MQ=255)
                                                                                                   ccGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTAc                                                           <  1:940758/142‑1 (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                                                         <  2:474047/114‑1 (MQ=255)
                                                                                                                                 ccTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCt                                                         >  1:474047/1‑114 (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                   >  2:556343/1‑150 (MQ=255)
                                                                                                                                   tCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCa                   >  1:1588886/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                                                                       >  1:1122248/1‑92 (MQ=255)
                                                                                                                                         accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCgcggc                                                                       <  2:1122248/92‑1 (MQ=255)
                                                                                                                                         accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                               >  1:1691695/1‑132 (MQ=255)
                                                                                                                                         accaGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCt                               <  2:1691695/132‑1 (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)
                                                                                                                                                     |                                                                                                                                                     
CAACCGCGACTTCTGCTGGTGGCTGGGCGTGCTGGGCGAGTGGGACGCGGAAATTGCCAAGCCCGGTCGCGAGCATGTGACCATCGCGGTCAGTGGCGCCCGTGGTGGCCACACCGTGGACTTCCGTGCCCTCGCCCACCAGGGCATGACCCTGGTCGGCCTGACCCAGTCGTTCGAAAACGGCGTCGCGCGCTTCCAGGACAACCTGGTCGAGAACATCAATCGCGGCGACGACAACTACCTGGCCCTGCTGGATGCGGCCGATGCCTACATCGAACGCAACGGCCTGGACCTGCCG  >  NC_002947/4054094‑4054391

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: