Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 797,512 G→A 17.3% intergenic (+33/‑44) PP_0682 → / → PP_0683 inner membrane protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947797,5120GA17.3% 159.7 / 27.5 75intergenic (+33/‑44)PP_0682/PP_0683inner membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/27);  new base A (7/6);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

TTCATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAA  >  NC_002947/797363‑797659
                                                                                                                                                     |                                                                                                                                                   
ttCATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTgggg                                                                                                                                                     <  2:780124/150‑1 (MQ=255)
   aTGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcg                                                                                                                                                  >  2:1339056/1‑150 (MQ=255)
         gcgcAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATg                                                                                                                                            <  1:98423/150‑1 (MQ=255)
           gcAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGcc                                                                                                                                          >  1:945707/1‑150 (MQ=255)
            cAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGccc                                                                                                                                         >  1:1337276/1‑150 (MQ=255)
                 cATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCAGTTGGGGCGGCCATGccaggtcg                                                                                                                                    >  1:384882/1‑144 (MQ=255)
                             ttACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                        <  1:759860/150‑1 (MQ=255)
                                   gATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAggag                                                                                                                            <  1:1032062/140‑4 (MQ=255)
                                   gATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAggag                                                                                                                            >  2:1032062/1‑137 (MQ=255)
                                         gTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             >  2:583056/1‑131 (MQ=255)
                                         gTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             <  1:583056/133‑3 (MQ=255)
                                             gTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGt                                                                                                        <  2:945707/150‑1 (MQ=255)
                                             gTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGt                                                                                                        <  1:1339056/150‑1 (MQ=255)
                                               gATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTgg                                                                                                      >  1:344141/1‑150 (MQ=255)
                                                  catcCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgt                                                                                                   <  2:434910/150‑1 (MQ=255)
                                                       acaTGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtcagtc                                                                                              >  2:1249162/1‑150 (MQ=255)
                                                       acaTGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCgtcagtc                                                                                              <  2:1337276/150‑1 (MQ=255)
                                                           gACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcgg                                                                                                                                                 <  2:1037801/95‑1 (MQ=255)
                                                           gACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcgg                                                                                                                                                 >  1:1037801/1‑95 (MQ=255)
                                                             ccTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGTCAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCt                                                                                        >  2:1639518/1‑150 (MQ=255)
                                                              cTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGc                                                                                                     >  2:1181054/1‑136 (MQ=255)
                                                              cTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGc                                                                                                     <  1:1181054/136‑1 (MQ=255)
                                                              cTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCGGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             <  2:676543/112‑3 (MQ=39)
                                                              cTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCGGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             >  1:676543/1‑110 (MQ=39)
                                                                   tGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCa                                                                                  >  2:391948/1‑150 (MQ=255)
                                                                      tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGc                                                                                                                                       >  2:83074/1‑94 (MQ=255)
                                                                      tctcGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGc                                                                                                                                       <  1:83074/94‑1 (MQ=255)
                                                                          ggCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTccc                                                                           >  2:1507533/1‑150 (MQ=255)
                                                                               tCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCt                                                                      <  2:344141/150‑1 (MQ=255)
                                                                                     tGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTc                                                                        >  2:754835/1‑142 (MQ=255)
                                                                                     tGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTc                                                                        <  1:754835/142‑1 (MQ=255)
                                                                                        gCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCcaca                                                                                                                                >  2:1473111/1‑83 (MQ=255)
                                                                                        gCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCcaca                                                                                                                                <  1:1473111/83‑1 (MQ=255)
                                                                                          tACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtc                                                                               <  2:647300/130‑1 (MQ=255)
                                                                                          tACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATtc                                                                               >  1:647300/1‑130 (MQ=255)
                                                                                           aCCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGa                                                                                                          <  1:1426350/102‑1 (MQ=255)
                                                                                           aCCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGa                                                                                                          >  2:1426350/1‑102 (MQ=255)
                                                                                            ccTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcg                                                                                                                                                  <  1:260714/61‑1 (MQ=255)
                                                                                            ccTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGgggcg                                                                                                                                                  >  2:260714/1‑61 (MQ=255)
                                                                                             cTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                        <  1:462955/86‑1 (MQ=255)
                                                                                             cTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGa                                                                                                                        >  2:462955/1‑86 (MQ=255)
                                                                                             cTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             <  2:1516262/81‑3 (MQ=37)
                                                                                             cTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGAGCGGGCATGCCCGCCCCCACAgga                                                                                                                             >  1:1516262/1‑79 (MQ=37)
                                                                                              tGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGc                                                       >  1:397258/1‑150 (MQ=255)
                                                                                                        gACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGcc                                                                                         >  1:666843/1‑106 (MQ=255)
                                                                                                        gACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGcc                                                                                         <  2:666843/106‑1 (MQ=255)
                                                                                                         aCCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTg                                            <  1:1639518/150‑1 (MQ=255)
                                                                                                           cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTgcgc                                                                                   >  2:111837/1‑109 (MQ=255)
                                                                                                           cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTgcgc                                                                                   <  1:111837/109‑1 (MQ=255)
                                                                                                           cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGt                                          >  2:1394972/1‑150 (MQ=255)
                                                                                                           cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCCCTCGCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGt                                          >  1:836204/1‑150 (MQ=255)
                                                                                                           cAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTACGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGt                                          >  2:1400618/1‑150 (MQ=255)
                                                                                                                aCTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCa                                     >  1:1135430/1‑150 (MQ=255)
                                                                                                                aCTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCa                                     >  1:124473/1‑150 (MQ=255)
                                                                                                                     gTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTg                                            <  2:1603935/138‑1 (MQ=255)
                                                                                                                     gTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTg                                            >  1:1603935/1‑138 (MQ=255)
                                                                                                                        cccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACaaaaaa                                                            <  1:1453423/119‑1 (MQ=255)
                                                                                                                        cccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACaaaaaa                                                            >  2:1453423/1‑119 (MQ=255)
                                                                                                                        cccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCCCCcaca                                                                                                                                >  2:95389/1‑51 (MQ=255)
                                                                                                                        cccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCCCCcaca                                                                                                                                <  1:95389/51‑1 (MQ=255)
                                                                                                                         ccccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGt                            >  1:869222/1‑150 (MQ=255)
                                                                                                                          cccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTg                           >  2:2732/1‑150 (MQ=255)
                                                                                                                          cccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTg                           >  1:702575/1‑150 (MQ=255)
                                                                                                                           ccAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGc                          >  2:87400/1‑150 (MQ=255)
                                                                                                                             aGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGctc                        <  2:702575/150‑1 (MQ=255)
                                                                                                                               gCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCtctc                      <  1:1507533/150‑1 (MQ=255)
                                                                                                                                ccTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTg                                            >  1:1152099/1‑127 (MQ=255)
                                                                                                                                ccTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTg                                            <  2:1152099/127‑1 (MQ=255)
                                                                                                                                 cTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCt                                                                          >  1:560230/1‑96 (MQ=255)
                                                                                                                                 cTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCt                                                                          <  2:560230/96‑1 (MQ=255)
                                                                                                                                  tCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATg                   >  2:328083/1‑150 (MQ=255)
                                                                                                                                    aCTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGaa                 >  1:264463/1‑150 (MQ=255)
                                                                                                                                     cTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAAc                <  1:391948/150‑1 (MQ=255)
                                                                                                                                      tGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTgc                                                                                     >  1:1108279/1‑80 (MQ=255)
                                                                                                                                      tGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTgc                                                                                     <  2:1108279/80‑1 (MQ=255)
                                                                                                                                           cTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCAt          <  1:1249162/150‑1 (MQ=255)
                                                                                                                                           cTCCTGCGGTGGCGCGAATGCCCGCTCCCACCGATATCGGGGTCCGACGGACGAGTGGCGGCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCAt          <  2:836204/150‑1 (MQ=255)
                                                                                                                                            tCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATc         >  1:1083562/1‑150 (MQ=255)
                                                                                                                                               tgtgGGAGCGGGCATGCCCGCCCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGc      >  1:1414728/1‑150 (MQ=255)
                                                                                                                                                 tgGGAGCGGGCATGCCCGCCCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTc                                                                        <  1:139850/82‑1 (MQ=255)
                                                                                                                                                 tgGGAGCGGGCATGCCCGCCCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTc                                                                        >  2:139850/1‑82 (MQ=255)
                                                                                                                                                   ggggCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAAt                                                           >  1:68839/1‑93 (MQ=255)
                                                                                                                                                   ggggCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAAt                                                           <  2:68839/93‑1 (MQ=255)
                                                                                                                                                   ggggCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCaa  >  2:1424394/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCATGGATGCGCAGAACATCTCCACCGATTACCTGATGTGGTACGTGATCATCCACATGACCTTCGTGGTCTCGGCGTTCTGCATGGGCTACCTGGACAAGCTGACCAAGCACTGAGTCCCCCCAGGCCTCACTGCAACTCCTGTGGGGGCGGGCATGCCCGCTCCCACAGATATCGAGGCTCAACTGACGAGTGGCGTCAGTCATGCCTTGCGCATTCTCCCTTCCTGTACAAAAAATGAGCTCTCATGCGTGGTTCCCATCGCGAGGTGCTCTCATGAACCTGCATCAGCTCAA  >  NC_002947/797363‑797659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑