Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 337,720 G→T 88.2% E295D (GAG→GAT yahJ → putative metallo‑dependent hydrolase domain deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273337,7200GT88.2% 390.9 / 49.1 166E295D (GAG→GATyahJputative metallo‑dependent hydrolase domain deaminase
Reads supporting (aligned to +/- strand):  ref base G (10/10);  new base T (75/71);  total (85/81)
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.53e-01

Only 100 of 168 total aligned reads displayed.

CGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTG  >  CP009273/337575‑337864
                                                                                                                                                 |                                                                                                                                                
cGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaa                                                                                                                                                >  2:1700653/1‑148 (MQ=255)
cGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATaa                                                                                                                                                >  2:1700691/1‑148 (MQ=255)
   ccAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                             >  1:81642/1‑148 (MQ=255)
   ccAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                             <  2:2807146/148‑1 (MQ=255)
   ccAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAAc                                                                                                                                             <  2:2807149/148‑1 (MQ=255)
       tgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGATACGGTAGATAAAACGCca                                                                                                                                         >  2:1149838/1‑148 (MQ=255)
       tgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCca                                                                                                                                         >  1:407688/1‑148 (MQ=255)
       tgtTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCACCAATTATATGGTTGAAACGGTAGATAAAACGCca                                                                                                                                         >  2:1260824/1‑148 (MQ=255)
         tttttGGCGCTATGTAAAAATCCCTCGACACATTTTTCCTTATTGCGCTGTACTACCACAATGGGGTCTATATTCACCTCCACGAAACCACTCCGGCAGGCGTGGACGCAATAAATTATATGGTTGAAACGGTAGATAAAACGCcaca                                                                                                                                       <  1:1529383/144‑1 (MQ=255)
           gATGGCGCGATGTAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGAAAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGc                                                                                                                                     <  1:1571776/148‑1 (MQ=255)
              ggCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGa                                                                                                                                  <  2:1330064/148‑1 (MQ=255)
              ggCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGa                                                                                                                                  >  2:2958646/1‑148 (MQ=255)
               gcgcGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGaa                                                                                                                                 <  1:1318745/148‑1 (MQ=255)
                  cgATGCAAAAATCCCTCGACACCAGGTTCCAGAGTGCGCTGGACTACGACAATGGCGTCTATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATAAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAggg                                                                                                                              <  1:1149241/148‑1 (MQ=255)
                    aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCa                                                                                                                            >  2:2195229/1‑148 (MQ=255)
                    aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCa                                                                                                                            >  2:1203765/1‑148 (MQ=255)
                    aTGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCa                                                                                                                            >  2:1866547/1‑148 (MQ=255)
                      gtaaaaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                          <  1:369487/146‑1 (MQ=255)
                      ggAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                          <  2:2687511/148‑1 (MQ=255)
                      ggAAAAATCCCTCGACACCATGTGCCAGATTGCGCTGGACTCCGACAAAGGCTTCGATATTCACCTGCACTAAACCACTCCGGCAGGCGTGGCAGCCATAAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAg                                                                                                                          <  1:2276542/148‑1 (MQ=255)
                        aaaaaTCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCCATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCt                                                                                                                        >  1:1644869/1‑148 (MQ=255)
                          aaaTCCCTCGACACCATGTTCCATATTGCGGTGGACTACTAAAAAGGCTTCGATATTCACCTGCACGAAACCACTCCGGCAGGGGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGa                                                                                                                      <  1:2043000/148‑1 (MQ=255)
                              cccTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCAt                                                                                                                  >  2:544580/1‑148 (MQ=255)
                                   gACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTc                                                                                                             <  2:1909977/148‑1 (MQ=255)
                                    acacCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCa                                                                                                            <  2:1641205/148‑1 (MQ=255)
                                     cacCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCACGCTGACCATCAGTCAc                                                                                                           >  1:2342747/1‑148 (MQ=255)
                                       ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGc                                                                                                         <  1:3042518/148‑1 (MQ=255)
                                       ccATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGc                                                                                                         >  2:677111/1‑148 (MQ=255)
                                        cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAACACGCCACAGCCGAAGTGCAAGCTGACCATCAGTCACGcc                                                                                                        <  2:1119053/148‑1 (MQ=255)
                                        cATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGcc                                                                                                        >  2:599154/1‑148 (MQ=255)
                                         aTGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACGCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCt                                                                                                       <  2:2788471/148‑1 (MQ=255)
                                             tCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgc                                                                                                   >  1:2944783/1‑148 (MQ=255)
                                             tCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTgc                                                                                                   >  1:2943758/1‑148 (MQ=255)
                                                  aTTTCGCTGGACTACGACAAAGGCGTCGATATTCACCCGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTgg                                                                                              >  2:889828/1‑148 (MQ=255)
                                                    tGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCa                                                                                            <  2:1655485/148‑1 (MQ=255)
                                                        cTGGACTACGACACAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGc                                                                                        <  2:3309529/148‑1 (MQ=255)
                                                              tACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACg                                                                                  <  1:2823812/148‑1 (MQ=255)
                                                              tACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACg                                                                                  >  1:2975659/1‑148 (MQ=255)
                                                              tACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACg                                                                                  <  2:2700308/148‑1 (MQ=255)
                                                               aCGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGa                                                                                 >  1:349493/1‑148 (MQ=255)
                                                                cGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAg                                                                                >  1:2498547/1‑148 (MQ=255)
                                                                  aCAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCa                                                                              <  1:1027524/148‑1 (MQ=255)
                                                                   cAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCaa                                                                             >  1:1895159/1‑148 (MQ=255)
                                                                     aaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACa                                                                           >  2:637978/1‑148 (MQ=255)
                                                                     aaGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACa                                                                           >  2:638822/1‑148 (MQ=255)
                                                                      aGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                          >  1:472345/1‑148 (MQ=255)
                                                                      aGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAg                                                                          >  1:472351/1‑148 (MQ=255)
                                                                        gCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAATGAGCAACAGGt                                                                        <  2:804811/148‑1 (MQ=255)
                                                                        gCGTCGATATTCACCGGCACGACACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGt                                                                        <  2:2003075/148‑1 (MQ=255)
                                                                         cGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTa                                                                       >  2:605305/1‑148 (MQ=255)
                                                                           tCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGa                                                                     <  1:3020392/148‑1 (MQ=255)
                                                                            cGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGAt                                                                    <  1:2647856/148‑1 (MQ=255)
                                                                                 ttCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACt                                                               >  2:195574/1‑148 (MQ=255)
                                                                                   cACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTgg                                                             >  1:2737990/1‑148 (MQ=255)
                                                                                   cACCGGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTgg                                                             <  1:668314/148‑1 (MQ=255)
                                                                                       tGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGaa                                                         <  2:2458188/148‑1 (MQ=255)
                                                                                       tGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGaa                                                         <  2:2458193/148‑1 (MQ=255)
                                                                                            gAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGa                                                    <  2:3036608/148‑1 (MQ=255)
                                                                                             aaaCCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGATGAAACGTTAGATAAAACTCCAAACCTGAAGGGCAAGCTGACCATCAGTCACGCCTATGAGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGag                                                   >  2:2060142/1‑147 (MQ=255)
                                                                                                ccACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtggt                                                <  1:1269297/148‑1 (MQ=255)
                                                                                                 cACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAtggtg                                               <  2:1241890/148‑1 (MQ=255)
                                                                                                    tCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGgtggtg                                            <  2:2125688/148‑1 (MQ=255)
                                                                                                        gCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCACCAGGTAGATGAACTGGCGAACCGGATGGTGGTGcaac                                        >  2:2217055/1‑148 (MQ=255)
                                                                                                        gCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGAGGGTGGTGccac                                        >  2:1725891/1‑148 (MQ=255)
                                                                                                           ggCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACaaa                                     >  1:628692/1‑148 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    <  1:2194664/148‑1 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    >  1:1449961/1‑148 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    >  1:1450256/1‑148 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    <  1:2195229/148‑1 (MQ=255)
                                                                                                            gCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAt                                    >  1:2472457/1‑148 (MQ=255)
                                                                                                              gTGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAAttt                                  <  1:138292/148‑1 (MQ=255)
                                                                                                               tGGCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTc                                 <  1:2643281/148‑1 (MQ=255)
                                                                                                                ggCAGCAATCAATTATATGGTTGAATCGGTAGATAAACCGCCACAGCTGAAGGGCAAGCTGTCCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCATGAAATTTCt                                <  2:1018631/148‑1 (MQ=255)
                                                                                                                 gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTa                               <  1:1700668/148‑1 (MQ=255)
                                                                                                                 gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTa                               <  1:186867/148‑1 (MQ=255)
                                                                                                                 gCAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTa                               >  2:2614798/1‑148 (MQ=255)
                                                                                                                  cAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                              <  1:3097409/148‑1 (MQ=255)
                                                                                                                  cAGCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTAt                              <  1:291751/148‑1 (MQ=255)
                                                                                                                    gCCATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCg                            <  1:1926396/148‑1 (MQ=255)
                                                                                                                     ccATCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           <  1:473967/148‑1 (MQ=255)
                                                                                                                     ccATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTTACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           <  2:407688/148‑1 (MQ=255)
                                                                                                                     ccATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGc                           <  2:409180/148‑1 (MQ=255)
                                                                                                                       aTCAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCt                         >  2:3084594/1‑148 (MQ=255)
                                                                                                                         cAATTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCg                       >  2:1808399/1‑148 (MQ=255)
                                                                                                                         cAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCg                       >  1:2420981/1‑148 (MQ=255)
                                                                                                                          aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGa                      >  1:1090109/1‑148 (MQ=255)
                                                                                                                          aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGa                      <  2:2442340/148‑1 (MQ=255)
                                                                                                                          aaTTATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGa                      <  2:1776411/148‑1 (MQ=255)
                                                                                                                            ttATATGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACg                    <  1:407442/148‑1 (MQ=255)
                                                                                                                             tataTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                   >  1:265121/1‑148 (MQ=255)
                                                                                                                             tataTGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACgg                   >  1:450408/1‑148 (MQ=255)
                                                                                                                               tatGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTg                 >  1:1957233/1‑148 (MQ=255)
                                                                                                                                 tGGTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGcc               <  1:2027901/148‑1 (MQ=255)
                                                                                                                                 gggTTGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGcc               <  1:362183/147‑1 (MQ=255)
                                                                                                                                    ttGAAACGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGAt            <  1:2017858/148‑1 (MQ=255)
                                                                                                                                       aaaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTgg         >  1:1762076/1‑148 (MQ=255)
                                                                                                                                       aaaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTgg         >  1:86616/1‑148 (MQ=255)
                                                                                                                                       aaaCGGTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACAATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTgg         >  1:1762074/1‑148 (MQ=255)
                                                                                                                                            gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGc    <  1:528031/148‑1 (MQ=255)
                                                                                                                                            gTAGATAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGc    >  1:630843/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGACCAGTGTTGATGGCGCGATGGAAAAATCCCTCGACACCATGTTCCAGATTGCGCTGGACTACGACAAAGGCGTCGATATTCACCTGCACGAAACCACTCCGGCAGGCGTGGCAGCCATCAATTATATGGTTGAAACGGTAGAGAAAACGCCACAGCTGAAGGGCAAGCTGACCATCAGTCACGCCTTTGCGCTGGCAACGCTCAACGAGCAACAGGTAGATGAACTGGCGAACCGGATGGTGGTGCAACAAATTTCTATCGCCTCGACGGTGCCGATTGGCACGCTG  >  CP009273/337575‑337864

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