Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,942,649 A→G 7.7% F81L (TTT→CTT)  nudB ← dihydroneopterin triphosphate pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,942,6490AG 7.7% 813.9 / 45.3 272F81L (TTT→CTT) nudBdihydroneopterin triphosphate pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (126/122);  new base G (11/10);  total (139/133)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 282 total aligned reads displayed.

CGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCGGGATCGT  >  CP009273/1942500‑1942798
                                                                                                                                                     |                                                                                                                                                     
  ccgcAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGct                                                                                                                                                     >  1:1021048/1‑150 (MQ=255)
   cgcAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAActc                                                                                                                                                    >  1:115400/1‑150 (MQ=255)
     cagcagGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTa                                                                                                                                                  <  1:722589/150‑1 (MQ=255)
     cagcagGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTa                                                                                                                                                  <  2:630014/150‑1 (MQ=255)
         agGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGt                                                                                                                                              <  1:2269732/150‑1 (MQ=255)
            cgcATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGCTCTACCGTgcg                                                                                                                                           >  1:866924/1‑150 (MQ=255)
               aTCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTg                                                                                                                                        >  2:3920303/1‑150 (MQ=255)
                tCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGa                                                                                                                                       <  1:1691445/150‑1 (MQ=255)
                tCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGa                                                                                                                                       >  1:3675624/1‑150 (MQ=255)
                  aaGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACa                                                                                                                                     >  1:4166121/1‑150 (MQ=255)
                  aaGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACa                                                                                                                                     <  2:4181287/150‑1 (MQ=255)
                  aaGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACa                                                                                                                                     >  2:3211416/1‑150 (MQ=255)
                     ccACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTc                                                                                                                                  >  1:553684/1‑150 (MQ=255)
                     ccACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTc                                                                                                                                  >  1:1683456/1‑150 (MQ=255)
                        cTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAAt                                                                                                                               >  2:5036140/1‑150 (MQ=255)
                         ttGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAAtt                                                                                                                              <  1:2478368/150‑1 (MQ=255)
                           gTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTaa                                                                                                                            >  1:1248237/1‑150 (MQ=255)
                             aaGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAgg                                                                                                                          >  1:1486654/1‑150 (MQ=255)
                             aaGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAgg                                                                                                                          >  2:3376167/1‑150 (MQ=255)
                             aaGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAgg                                                                                                                          >  2:3377357/1‑150 (MQ=255)
                              aGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGt                                                                                                                         >  2:832076/1‑150 (MQ=255)
                               gCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTc                                                                                                                        >  1:2842779/1‑150 (MQ=255)
                               gCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTc                                                                                                                        <  1:4721414/150‑1 (MQ=255)
                                 cAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAg                                                                                                                      <  2:414580/150‑1 (MQ=255)
                                   gATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAgtt                                                                                                                    <  1:4650485/150‑1 (MQ=255)
                                    aTGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAgttg                                                                                                                   <  1:3218453/150‑1 (MQ=255)
                                      gTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAgttgtt                                                                                                                 <  2:196022/150‑1 (MQ=255)
                                       ttCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTc                                                                                                                <  2:740497/150‑1 (MQ=255)
                                        tCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCa                                                                                                               >  1:82866/1‑150 (MQ=255)
                                         cAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAg                                                                                                              >  1:2185459/1‑150 (MQ=255)
                                          aGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGc                                                                                                             >  2:1787265/1‑150 (MQ=255)
                                            tGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGa                                                                                                           <  1:349041/150‑1 (MQ=255)
                                             gAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGac                                                                                                          >  2:2347898/1‑150 (MQ=255)
                                              aaaaCGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGaca                                                                                                         <  2:932910/150‑1 (MQ=255)
                                                aaCGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGacaac                                                                                                       <  1:812342/150‑1 (MQ=255)
                                                   gATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATc                                                                                                    >  2:1332441/1‑150 (MQ=255)
                                                       tGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATg                                                                                                <  2:671408/150‑1 (MQ=255)
                                                          cGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTg                                                                                             <  2:2569311/150‑1 (MQ=255)
                                                          cGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTg                                                                                             >  2:1133317/1‑150 (MQ=255)
                                                          cGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTg                                                                                             >  1:95911/1‑150 (MQ=255)
                                                           gCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGa                                                                                            >  1:44589/1‑150 (MQ=255)
                                                           gCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGa                                                                                            >  1:1793632/1‑150 (MQ=255)
                                                            cTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGAc                                                                                           <  1:65408/150‑1 (MQ=255)
                                                             tCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGAcc                                                                                          <  1:1858780/150‑1 (MQ=255)
                                                               gTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACctc                                                                                        >  2:2752590/1‑150 (MQ=255)
                                                                tGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACctct                                                                                       >  1:4741830/1‑150 (MQ=255)
                                                                tGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACctct                                                                                       <  1:4022998/150‑1 (MQ=255)
                                                                tGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACctct                                                                                       <  1:3744896/150‑1 (MQ=255)
                                                                tGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACctct                                                                                       <  2:3959745/150‑1 (MQ=255)
                                                                   ggAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTcc                                                                                    >  2:5110662/1‑150 (MQ=255)
                                                                      aGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCttt                                                                                 <  1:1062213/150‑1 (MQ=255)
                                                                        cgcAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTAc                                                                               <  1:2346629/150‑1 (MQ=255)
                                                                        cgcAAGACAGAACCATGAGTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTAc                                                                               <  2:443474/150‑1 (MQ=255)
                                                                           aaGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTc                                                                            <  1:1955866/150‑1 (MQ=255)
                                                                           aaGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTc                                                                            >  2:1488418/1‑150 (MQ=255)
                                                                               cAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCgcgc                                                                        >  2:3056302/1‑150 (MQ=255)
                                                                               cAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCgcgc                                                                        <  2:1291326/150‑1 (MQ=255)
                                                                                aGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCa                                                                       >  2:3648179/1‑150 (MQ=255)
                                                                                  aaCCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATg                                                                     <  2:2662894/150‑1 (MQ=255)
                                                                                    ccATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGc                                                                   >  1:5121627/1‑150 (MQ=255)
                                                                                     cATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCa                                                                  <  2:4972729/150‑1 (MQ=255)
                                                                                        gATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCt                                                               <  1:4229414/150‑1 (MQ=255)
                                                                                            ccGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTgcg                                                           <  1:481129/150‑1 (MQ=255)
                                                                                              gTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTgcggc                                                         <  1:596624/150‑1 (MQ=255)
                                                                                              gTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTgcggc                                                         >  2:1433126/1‑150 (MQ=255)
                                                                                               tATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTgcggcg                                                        <  2:1184736/150‑1 (MQ=255)
                                                                                                aTTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGgcgc                                                       >  2:2889753/1‑150 (MQ=255)
                                                                                                   aCGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAGCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGt                                                    >  2:4738639/1‑150 (MQ=255)
                                                                                                    cgcgTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGtt                                                   <  1:1577036/150‑1 (MQ=255)
                                                                                                     gcgTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAACCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGttt                                                  <  2:4840089/150‑1 (MQ=255)
                                                                                                     gcgTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGttt                                                  <  2:3876583/150‑1 (MQ=255)
                                                                                                     gcgTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGttt                                                  <  1:174862/150‑1 (MQ=255)
                                                                                                     gcgTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGttt                                                  >  1:1450887/1‑150 (MQ=255)
                                                                                                      cgTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTc                                                 <  2:1683456/150‑1 (MQ=255)
                                                                                                           cGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCct                                            <  2:573249/150‑1 (MQ=255)
                                                                                                               cGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTc                                        <  1:2672716/150‑1 (MQ=255)
                                                                                                               cGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTc                                        <  1:3932139/150‑1 (MQ=255)
                                                                                                               cGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAAGTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTc                                        <  2:2984216/150‑1 (MQ=255)
                                                                                                                ggCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTcc                                       >  1:1351514/1‑150 (MQ=255)
                                                                                                                 gcgcATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCa                                      <  1:5011646/150‑1 (MQ=255)
                                                                                                                  cgcATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCAc                                     >  1:153797/1‑150 (MQ=255)
                                                                                                                  cgcATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCAc                                     <  2:4725246/150‑1 (MQ=255)
                                                                                                                  cgcATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCAc                                     >  2:4631330/1‑150 (MQ=255)
                                                                                                                    cATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGc                                   >  2:3183953/1‑150 (MQ=255)
                                                                                                                     aTAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCt                                  <  1:2713363/150‑1 (MQ=255)
                                                                                                                     aTAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCt                                  <  2:2163055/150‑1 (MQ=255)
                                                                                                                         cGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCg                              >  2:1901598/1‑150 (MQ=255)
                                                                                                                               cGTAAATGTGAAAAAATTTCAAACTCTCCCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTAcc                        <  2:357323/150‑1 (MQ=255)
                                                                                                                               cGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTAcc                        <  2:4305991/150‑1 (MQ=255)
                                                                                                                                 tAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGa                      >  1:486556/1‑150 (MQ=255)
                                                                                                                                   aaTGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACt                    >  2:2371177/1‑150 (MQ=255)
                                                                                                                                     tgtgAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGc                  >  2:3485617/1‑150 (MQ=255)
                                                                                                                                          aaaaaTTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGaa             >  1:4786672/1‑150 (MQ=255)
                                                                                                                                           aaaaTTTCAAGCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGaaa            >  1:4234299/1‑150 (MQ=255)
                                                                                                                                            aaaTTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAAt           <  2:2079386/150‑1 (MQ=255)
                                                                                                                                              aTTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCg         <  1:2119753/150‑1 (MQ=255)
                                                                                                                                                  cAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCGGGAt     >  1:2556672/1‑150 (MQ=255)
                                                                                                                                                   aaaCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCGGGATc    <  1:2647919/150‑1 (MQ=255)
                                                                                                                                                    aaCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCGGGATCg   >  2:2185109/1‑150 (MQ=255)
                                                                                                                                                    aaCTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCGGGATCg   >  1:3017436/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACCATGATTCCGTATTACGCGTCACGCCCGGCGCATAGCGATGACGTAAATGTGAAAAAATTTCAAACTCTACCGTGCGCTGACAGTCAATTAAGGTCAGTTGTTCAGCGACAACATCAATGGTGACCTCTTCCTTTACTTCGCGCATGGCAGCTTGCGGCGCGGTTTCACCCTCTTCCACGCTGCCGGTTACCGACTGCCAGAAATCGGGATCGT  >  CP009273/1942500‑1942798

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: