Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,564,981 A→G 7.4% G130G (GGT→GGC pitB ← phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,564,9810AG 7.4% 265.3 / 5.9 90G130G (GGT→GGCpitBphosphate transporter
Reads supporting (aligned to +/- strand):  ref base A (38/45);  new base G (3/4);  total (41/49)
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

CACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTAC  >  NC_002947/1564832‑1565127
                                                                                                                                                     |                                                                                                                                                  
caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGGGAGCACCAGCGACATGGCGATATCGATCGCTTTCTGACAGTTGACGCCCTCGCCCAGCGGAATGTCGTTGATCATGGCGTTGGCCAGGCCAaca                                                                                                                                                    >  2:2304580/1‑150 (MQ=255)
caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAaca                                                                                                                                                    >  1:1085814/1‑150 (MQ=255)
caccaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAaca                                                                                                                                                    >  2:1237489/1‑150 (MQ=255)
 accaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAacac                                                                                                                                                   >  1:2242016/1‑150 (MQ=255)
 accaTTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAacac                                                                                                                                                   >  2:1155609/1‑150 (MQ=255)
     tttCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                               <  2:1327067/150‑1 (MQ=255)
       tCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGATCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGa                                                                                                                                             <  2:1430454/150‑1 (MQ=255)
       tCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGa                                                                                                                                             <  1:1183004/150‑1 (MQ=255)
         aGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATg                                                                                                                                           >  2:2091278/1‑150 (MQ=255)
          ggCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGGGACCACCAGCGACATGGCGATATCGATCGCTTTATGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCTAGGCCAACATCGAGGATgg                                                                                                                                          >  1:514067/1‑150 (MQ=255)
          ggCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATgg                                                                                                                                          >  2:599180/1‑150 (MQ=255)
           gCCGATCAGCACCAGGGCAGCAAAGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGa                                                                                                                                         >  1:642707/1‑150 (MQ=255)
              gATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGACGGTATCGATCGCTTTCTGCGAGTTGACGCCGTCGCCCAGCGGAATATCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGCAAcc                                                                                                                                      >  2:444289/1‑150 (MQ=255)
               aTCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGTGTTGGCCTGGCCAACGcc                                                                                                                                                  <  1:2254998/137‑1 (MQ=255)
               aTCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGTGTTGGCCTGGCCAACGcc                                                                                                                                                  >  2:2254998/1‑137 (MQ=255)
                tCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGa                                                                                                                                    <  1:1155609/150‑1 (MQ=255)
                  aGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATc                                                                                                                                  >  2:1325291/1‑150 (MQ=255)
                   gCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCa                                                                                                                                 <  2:2028765/150‑1 (MQ=255)
                     aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGc                                                                                                                               <  1:679118/150‑1 (MQ=255)
                     aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGc                                                                                                                               >  1:2081709/1‑150 (MQ=255)
                     aCCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAATTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGc                                                                                                                               >  1:883038/1‑150 (MQ=255)
                        aGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtg                                                                                                                            >  2:244498/1‑150 (MQ=255)
                         gggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtgt                                                                                                                           <  1:2202070/150‑1 (MQ=255)
                         gggCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAATTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCgtgt                                                                                                                           <  2:883038/150‑1 (MQ=255)
                               caACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTTGCCAGGCCAACACCGAGGATGGAAc                                                                                                                                       >  2:1251973/1‑132 (MQ=255)
                               caACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTTGCCAGGCCAACACCGAGGATGGAAc                                                                                                                                       <  1:1251973/132‑1 (MQ=255)
                                   ggCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGa                                                                                                                 >  2:2025827/1‑150 (MQ=255)
                                        aGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCgg                                                                                                            <  2:2333794/150‑1 (MQ=255)
                                                 ttGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGAAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaa                                                                                                   >  1:1654609/1‑150 (MQ=255)
                                                 ttGGCGAGACCACCAGCGACATGGCCCGAGCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaa                                                                                                   >  1:1230039/1‑150 (MQ=255)
                                                  tGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCaaa                                                                                                  <  2:642707/150‑1 (MQ=255)
                                                       agaCCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAcc                                                                                             <  2:2151967/150‑1 (MQ=255)
                                                         accaccAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAg                                                                                           >  1:1669218/1‑150 (MQ=255)
                                                          ccaccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTAGCCCGGCTGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAgg                                                                                          <  2:1230039/150‑1 (MQ=255)
                                                           caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                               >  1:864870/1‑96 (MQ=255)
                                                           caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                               <  2:864870/96‑1 (MQ=255)
                                                           caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAAc                                                                                                                                       >  1:749650/1‑104 (MQ=255)
                                                           caccaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAAc                                                                                                                                       <  2:749650/104‑1 (MQ=255)
                                                            accaGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                                        >  1:903762/1‑150 (MQ=255)
                                                                  gACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAgg                                                                                  <  1:1325291/150‑1 (MQ=255)
                                                                   aCATGGCGATATCGATCGCTTTCTGCCAGTTCACGTCGTCGCCCGGCGGAATGGCGTTGATCGGGGCGTTGGCCAGGCAAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGt                                                                                 <  1:2304580/150‑1 (MQ=255)
                                                                     aTGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGCTGCCACCGTACTCGCTACCCTGGTTc                                                                               >  2:501015/1‑150 (MQ=255)
                                                                      tGGCGATATCGCCCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGCCGTTCTCCAGGCTCACAAAGAGGATGGAACCGATCAGCGTGTGAGAACTGTAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTcc                                                                              <  1:1333055/150‑1 (MQ=255)
                                                                      tGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTcc                                                                              >  1:897277/1‑150 (MQ=255)
                                                                        gCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAATTACCAGGTGCCCAGGTTCCAg                                                                            <  1:1237489/150‑1 (MQ=255)
                                                                          gATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGt                                                                          <  1:2091278/150‑1 (MQ=255)
                                                                          atcttcgatcgCCCTCTGAGAGGTGACCCCCCCGCCCAGCTGATTGTCTTTGAGCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGt                                                                          <  2:1668121/146‑1 (MQ=255)
                                                                              tcgatcgCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATg                                                                      <  1:501015/150‑1 (MQ=255)
                                                                               cgatcgCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGagat                                                                     >  1:1668121/1‑147 (MQ=255)
                                                                                gatcgCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGc                                                                    >  2:1442298/1‑150 (MQ=255)
                                                                                      tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCa                                                                             <  1:2147808/135‑1 (MQ=255)
                                                                                      tttCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCa                                                                             >  2:2147808/1‑135 (MQ=255)
                                                                                        tCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcagcca                                                            <  2:897277/150‑1 (MQ=255)
                                                                                          tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                               <  1:481753/65‑1 (MQ=255)
                                                                                          tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAg                                                                                                                                               >  2:481753/1‑65 (MQ=255)
                                                                                          tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAcc                                                                                             >  1:1050851/1‑115 (MQ=255)
                                                                                          tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTAcc                                                                                             <  2:1050851/115‑1 (MQ=255)
                                                                                          tGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGcagccagc                                                          >  2:1934528/1‑150 (MQ=255)
                                                                                               gTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACga                                                     <  1:244498/150‑1 (MQ=255)
                                                                                                  gACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                                        <  2:977697/112‑1 (MQ=255)
                                                                                                  gACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTg                                                                                        >  1:977697/1‑112 (MQ=255)
                                                                                                  gACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaa                                                  <  2:546334/150‑1 (MQ=255)
                                                                                                    cGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaaga                                                <  1:2025827/150‑1 (MQ=255)
                                                                                                     gCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGAGCAGCGTGTGCGAACTGGAGGCGGGCATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGAc                                               <  1:315114/150‑1 (MQ=255)
                                                                                                       cGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGc                                                                                                      >  1:668709/1‑93 (MQ=255)
                                                                                                       cGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGc                                                                                                      <  2:668709/93‑1 (MQ=255)
                                                                                                        gTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCAt                                            >  1:1643889/1‑150 (MQ=255)
                                                                                                         tCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                                           <  2:2081709/150‑1 (MQ=255)
                                                                                                         tCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGCATGGAACCGCTTACCGCGCACGATCTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                                           <  1:2156175/150‑1 (MQ=255)
                                                                                                           gCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGc                                         <  1:599180/150‑1 (MQ=255)
                                                                                                           gCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGc                                         >  2:1000245/1‑150 (MQ=255)
                                                                                                                    aaTGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTg                                >  2:2164167/1‑150 (MQ=255)
                                                                                                                    aaTGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCACGGCCTGGCCATg                                <  1:1749556/150‑1 (MQ=255)
                                                                                                                        tCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTACGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                                           >  1:1232949/1‑135 (MQ=255)
                                                                                                                        tCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTACGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATg                                           <  2:1232949/135‑1 (MQ=255)
                                                                                                                        ctgTTGATCAGATCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCg                            <  2:722235/148‑1 (MQ=255)
                                                                                                                          gTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACTATGGCCAGGCCGTGGccatg                          <  1:968341/150‑4 (MQ=255)
                                                                                                                          gTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgt                          >  2:1862083/1‑150 (MQ=255)
                                                                                                                           ttGGTCAGGGCGGTGGCAAGGCCAACACCGAGGATGGAACGGATCAGCGTGTGGGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCTAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGgtg                         <  2:2242016/150‑1 (MQ=255)
                                                                                                                              aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGAt                                                                       >  1:90365/1‑101 (MQ=255)
                                                                                                                              aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTc                      <  1:1000245/150‑1 (MQ=255)
                                                                                                                              aTCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGACCAGGTTCCAGGTGAt                                                                       <  2:90365/101‑1 (MQ=255)
                                                                                                                               tCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGATCCAAGGTACCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTca                     >  1:1795643/1‑150 (MQ=255)
                                                                                                                                ctcggtgTTGGCCTGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCACGTGCCCAGGTTCGAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTcac                    >  2:665309/7‑150 (MQ=255)
                                                                                                                                 tcggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACga                                                     >  1:2081266/6‑116 (MQ=255)
                                                                                                                                 tcggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACga                                                     <  2:2081266/111‑1 (MQ=255)
                                                                                                                                   ggtgTTGGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAg                                                                            >  2:736946/4‑91 (MQ=255)
                                                                                                                                   ggtgTTCGCCTGGCCAACGCCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAg                                                                            <  1:736946/88‑1 (MQ=255)
                                                                                                                                   ggCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACAtt                 <  1:1934528/150‑1 (MQ=255)
                                                                                                                                          gccaggccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaa                                                  >  2:999903/1‑110 (MQ=255)
                                                                                                                                          gccaggccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGagaa                                                  <  1:999903/110‑1 (MQ=255)
                                                                                                                                               gccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGc                                    <  2:1205101/119‑1 (MQ=255)
                                                                                                                                               gccaACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGc                                    >  1:1205101/1‑119 (MQ=255)
                                                                                                                                                 caACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTa   >  2:477151/1‑150 (MQ=255)
                                                                                                                                                  aaCACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTTTAc  <  2:514067/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CACCATTTCAGGCCGATCAGCACCAGGGCAGCAACGGCGAAGCCGGCCATTGGCGAGACCACCAGCGACATGGCGATATCGATCGCTTTCTGCCAGTTGACGCCGTCGCCCAGCGGAATGTCGTTGATCAGGGCGTTGGCCAGGCCAACACCGAGGATGGAACCGATCAGCGTGTGCGAACTGGAGGCGGGGATGCCAAAGTACCAGGTGCCCAGGTTCCAGGTGATGGCCGCAGCCAGCAACGAGAAGACCATGGCCAGGCCGTGGCCGGTGTTCACATTGATCAGCAGCTCTAC  >  NC_002947/1564832‑1565127

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