Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,978,810 G→T 15.6% I186I (ATC→ATA aslA ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,978,8100GT15.6% 829.2 / 89.3 258I186I (ATC→ATAaslAputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (127/90);  new base T (17/24);  total (144/114)
Fisher's exact test for biased strand distribution p-value = 5.85e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

GAACGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGA  >  CP009273/3978661‑3978959
                                                                                                                                                     |                                                                                                                                                     
gAACGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCg                                                                                                                                                       <  2:3892336/150‑1 (MQ=255)
     gTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGc                                                                                                                                                  >  2:3663023/1‑150 (MQ=255)
         ggACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGa                                                                                                                                              >  1:4155706/1‑150 (MQ=255)
         ggACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGa                                                                                                                                              <  1:2550808/150‑1 (MQ=255)
          gACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAg                                                                                                                                             >  1:1238896/1‑150 (MQ=255)
            cTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTg                                                                                                                                           <  2:1908078/150‑1 (MQ=255)
                  gCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAg                                                                                                                                     <  1:3575921/150‑1 (MQ=255)
                    cACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGcc                                                                                                                                   >  1:1972789/1‑150 (MQ=255)
                       ttCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTg                                                                                                                                <  1:1421974/150‑1 (MQ=255)
                          cGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATc                                                                                                                             >  1:1229389/1‑150 (MQ=255)
                          cGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATc                                                                                                                             >  1:4114651/1‑150 (MQ=255)
                           ggATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCg                                                                                                                            <  1:653884/150‑1 (MQ=255)
                           ggATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGGGTGACGTAGCCCTGATCg                                                                                                                            >  2:3338552/1‑150 (MQ=255)
                               tGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTgca                                                                                                                        >  1:455088/1‑150 (MQ=255)
                               tGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTgca                                                                                                                        <  2:394896/150‑1 (MQ=255)
                                gACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTgcag                                                                                                                       <  2:3089949/150‑1 (MQ=255)
                                gACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTgcag                                                                                                                       >  1:3860119/1‑150 (MQ=255)
                                  cGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTgcagca                                                                                                                     <  2:2419641/150‑1 (MQ=255)
                                  cGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTgcagca                                                                                                                     >  2:2903635/1‑150 (MQ=255)
                                   gTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTACTGCGGCTGAGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCaa                                                                                                                    >  2:1094345/1‑150 (MQ=255)
                                     gAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATTCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATAGTGCAGCAACt                                                                                                                  >  1:2184456/1‑150 (MQ=255)
                                     gAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACt                                                                                                                  <  2:1710783/150‑1 (MQ=255)
                                     gAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACt                                                                                                                  >  2:214503/1‑150 (MQ=255)
                                     gAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGACACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACt                                                                                                                  >  2:1187422/1‑150 (MQ=255)
                                        cGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTGTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGca                                                                                                               <  2:2969999/150‑2 (MQ=255)
                                           cgcgCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCa                                                                                                            <  2:535261/150‑1 (MQ=255)
                                           cgcgCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCa                                                                                                            >  2:1211783/1‑150 (MQ=255)
                                              gcCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCg                                                                                                         >  2:1606101/1‑150 (MQ=255)
                                               ccATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGt                                                                                                        >  2:4060586/1‑150 (MQ=255)
                                                cATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAATGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTg                                                                                                       <  1:2430789/150‑1 (MQ=255)
                                                cATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTg                                                                                                       <  2:2702754/150‑1 (MQ=255)
                                                  ttCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGt                                                                                                     <  2:3944912/150‑1 (MQ=255)
                                                   tCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGtt                                                                                                    >  1:2251023/1‑150 (MQ=255)
                                                    cGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTa                                                                                                   >  1:352609/1‑150 (MQ=255)
                                                     ggTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTaa                                                                                                  >  1:1416328/1‑150 (MQ=255)
                                                     ggTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTaa                                                                                                  >  2:1631660/1‑150 (MQ=255)
                                                      gtgtACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAAc                                                                                                 <  2:1077182/150‑1 (MQ=255)
                                                      gtgtACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAAc                                                                                                 >  2:3621384/1‑150 (MQ=255)
                                                      gtgtACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAAc                                                                                                 >  1:1391483/1‑150 (MQ=255)
                                                        gtACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAAccc                                                                                               >  2:1842534/1‑150 (MQ=255)
                                                          aCATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCtt                                                                                             <  2:2204785/150‑1 (MQ=255)
                                                           cATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTg                                                                                            <  1:3957366/150‑1 (MQ=255)
                                                               tCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgc                                                                                        >  2:1288147/1‑150 (MQ=255)
                                                                cAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcc                                                                                       >  2:3619116/1‑150 (MQ=255)
                                                                cAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcc                                                                                       >  2:1976528/1‑150 (MQ=255)
                                                                cAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcc                                                                                       <  1:1066029/150‑1 (MQ=255)
                                                                cAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcc                                                                                       >  1:2584611/1‑150 (MQ=255)
                                                                 aGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgccc                                                                                      <  2:8367/150‑1 (MQ=255)
                                                                   acacCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGTGGCAGCGTGGTTAACCCTTGCAgcccgc                                                                                    >  1:3106558/1‑150 (MQ=255)
                                                                    cacCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAgcccgcc                                                                                   <  2:2994532/150‑1 (MQ=255)
                                                                        gAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGt                                                                               >  1:89315/1‑150 (MQ=255)
                                                                        gAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGt                                                                               >  1:1687859/1‑150 (MQ=255)
                                                                          gTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTg                                                                             >  2:1730470/1‑150 (MQ=255)
                                                                           ttAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGc                                                                            >  2:553213/1‑150 (MQ=255)
                                                                            tAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGcc                                                                           <  1:3340889/150‑1 (MQ=255)
                                                                                  cACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACa                                                                     <  2:422719/150‑1 (MQ=255)
                                                                                  cACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACa                                                                     >  2:659795/1‑150 (MQ=255)
                                                                                    cGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACAtt                                                                   <  2:1912014/150‑1 (MQ=255)
                                                                                     ggAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTg                                                                  <  2:2534281/150‑1 (MQ=255)
                                                                                      gAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTgg                                                                 >  2:486308/1‑150 (MQ=255)
                                                                                        aaTCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTggcg                                                               >  1:3153093/1‑150 (MQ=255)
                                                                                          tcatcaAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTggcggc                                                             >  2:3205004/1‑150 (MQ=255)
                                                                                           catcaAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTggcggta                                                            >  1:3645472/1‑148 (MQ=255)
                                                                                             tcaAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATc                                                          <  2:4125681/150‑1 (MQ=255)
                                                                                              caAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCa                                                         >  2:614933/1‑150 (MQ=255)
                                                                                              caAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCACCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCAGTACATTGGCGGCATCa                                                         >  2:3719367/1‑150 (MQ=255)
                                                                                               aaaGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAg                                                        >  1:544991/1‑150 (MQ=255)
                                                                                                 aGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGaa                                                      >  2:3512905/1‑150 (MQ=255)
                                                                                                   ccAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATg                                                    <  1:450207/150‑1 (MQ=255)
                                                                                                    cAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGc                                                   >  1:812274/1‑150 (MQ=255)
                                                                                                     aaCGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGcc                                                  >  2:593565/1‑150 (MQ=255)
                                                                                                       cGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgt                                                <  1:3053521/150‑1 (MQ=255)
                                                                                                       cGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgt                                                >  1:3268074/1‑150 (MQ=255)
                                                                                                         ttCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtgg                                              <  1:2886193/150‑1 (MQ=255)
                                                                                                         ttCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtgg                                              <  2:1000698/150‑1 (MQ=255)
                                                                                                           cTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtggtg                                            <  1:1740643/150‑1 (MQ=255)
                                                                                                           cTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCgtggtg                                            >  2:1397310/1‑150 (MQ=255)
                                                                                                                 cTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGtggatgga                                      >  1:555766/1‑150 (MQ=255)
                                                                                                                 cTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGtggatgga                                      >  2:3935464/1‑150 (MQ=255)
                                                                                                                    cGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTa                                   <  1:2774455/150‑1 (MQ=255)
                                                                                                                      aCTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTAtt                                 <  2:881909/150‑1 (MQ=255)
                                                                                                                       ctctTTGGTTTCCCCCATATGCCATTTTCCTATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTg                                <  2:2385280/150‑1 (MQ=255)
                                                                                                                           ttGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCg                            >  1:1771800/1‑150 (MQ=255)
                                                                                                                              ttttCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTg                         >  2:531030/1‑150 (MQ=255)
                                                                                                                               tttCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTga                        <  1:880319/150‑1 (MQ=255)
                                                                                                                               tttCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTga                        >  2:3988151/1‑150 (MQ=255)
                                                                                                                               tgtCCCCCCTCTGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTga                        <  2:1951081/148‑1 (MQ=255)
                                                                                                                                ttCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTgag                       >  2:2052817/1‑150 (MQ=255)
                                                                                                                                   ccccATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCACCTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAAt                    >  2:2814528/1‑150 (MQ=255)
                                                                                                                                        tatGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAAtggtgg               <  1:1124673/150‑1 (MQ=255)
                                                                                                                                         atGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGgc              >  2:2877468/1‑150 (MQ=255)
                                                                                                                                         atGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGgc              <  1:706018/150‑1 (MQ=255)
                                                                                                                                         atGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGgc              >  1:1924357/1‑150 (MQ=255)
                                                                                                                                          tGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGgcg             <  1:3696691/150‑1 (MQ=255)
                                                                                                                                              aTTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCggg         >  1:3736051/1‑150 (MQ=255)
                                                                                                                                              aTTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCggg         >  2:3669139/1‑150 (MQ=255)
                                                                                                                                                tttCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGtt       >  1:1626431/1‑150 (MQ=255)
                                                                                                                                                  tCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTgg     >  2:2238543/1‑150 (MQ=255)
                                                                                                                                                    cGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTgggg   <  2:2524641/150‑1 (MQ=255)
                                                                                                                                                    cGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTgggg   >  2:3015205/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAACGGTCCGGACTCAGGGCCACTTCCGGATTGACGTGAACGTCGCGCCATTCGGTGTACATATCAGACACCGAGTTAAAGCCACGGAAATCATCAAAGCCAACGTTCTGCGGCTGCGACTCTTTGTTTTCCCCCATATGCCATTTTCCGATGGCCTGAGTGACGTAGCCCTGATCGTGCAGCAACTGCGGCAGCGTGGTTAACCCTTGCAGCCCGCCCGGTTGCCCGTACATTGGCGGCATCAGAATGCCGTGGTGGATGGAGTATTGTCCGGTGAGAATGGTGGCGCGGGTTGGGGA  >  CP009273/3978661‑3978959

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