Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,718,959 C→A 11.8% L38F (TTG→TTT)  ysaB ← uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,718,9590CA11.8% 667.8 / 44.2 210L38F (TTG→TTT) ysaBuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base C (94/90);  new base A (17/8);  total (111/99)
Fisher's exact test for biased strand distribution p-value = 1.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 214 total aligned reads displayed.

AAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATT  >  CP009273/3718812‑3719106
                                                                                                                                                   |                                                                                                                                                   
aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTc                                                                                                                                                     >  1:389049/1‑148 (MQ=255)
aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTc                                                                                                                                                     <  1:2393402/148‑1 (MQ=255)
aaaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTc                                                                                                                                                     <  1:2159640/148‑1 (MQ=255)
 aaaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCa                                                                                                                                                    >  1:1268502/1‑148 (MQ=255)
  aaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCaa                                                                                                                                                   <  2:538204/148‑1 (MQ=255)
  aaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCaa                                                                                                                                                   <  2:2524045/148‑1 (MQ=255)
  aaGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCaa                                                                                                                                                   >  2:2678587/1‑148 (MQ=255)
       cAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCaacgaac                                                                                                                                              <  1:2325582/148‑1 (MQ=255)
          aCAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCa                                                                                                                                           >  2:1801531/1‑148 (MQ=255)
           cAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCag                                                                                                                                          <  2:2572371/148‑1 (MQ=255)
               cTGCTCGTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATACTTGCATAACGCTTCCATATTCAACGAAAGCAGAGGa                                                                                                                                      >  2:1471725/1‑148 (MQ=255)
                tGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGAt                                                                                                                                     <  1:2598829/148‑1 (MQ=255)
                tGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGAt                                                                                                                                     >  2:2184336/1‑148 (MQ=255)
                  ctctTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCa                                                                                                                                   <  1:34101/148‑1 (MQ=255)
                  ctctTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCa                                                                                                                                   >  1:405209/1‑148 (MQ=255)
                     tttACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAAcc                                                                                                                                <  2:2051353/148‑1 (MQ=255)
                         cGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTa                                                                                                                            >  1:2906402/1‑148 (MQ=255)
                         cGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTa                                                                                                                            <  2:2778241/148‑1 (MQ=255)
                         cGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCAGAGGATCAACCTTTa                                                                                                                            >  2:3265712/1‑148 (MQ=255)
                            aaGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTAccc                                                                                                                         >  2:2854401/1‑148 (MQ=255)
                             aGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCg                                                                                                                        <  2:1129569/148‑1 (MQ=255)
                             aGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCg                                                                                                                        >  2:269474/1‑148 (MQ=255)
                              ggTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAAAGTTTTGCCGCCAGATCCTTGCATAACGCTTCCATATTCAACGAAAGCAGAGGATCAACCTTTACCCGt                                                                                                                       >  1:838938/1‑148 (MQ=255)
                               gTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCAGAGGATCAACCATTACCCggt                                                                                                                      >  1:1431/1‑146 (MQ=255)
                                tataACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCAGAGGATCAACCTTTACCCGTtg                                                                                                                     >  1:1970436/1‑148 (MQ=255)
                                 ataACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTtgt                                                                                                                    >  2:1812965/1‑148 (MQ=255)
                                  taACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTtgtg                                                                                                                   >  1:2674978/1‑148 (MQ=255)
                                   aaCCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGc                                                                                                                  <  1:2500814/148‑1 (MQ=255)
                                    aCCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCt                                                                                                                 >  1:273441/1‑148 (MQ=255)
                                     ccGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCtt                                                                                                                >  1:238245/1‑148 (MQ=255)
                                      cgcgCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCttt                                                                                                               <  2:2365790/148‑1 (MQ=255)
                                          cATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGc                                                                                                           >  2:2321937/1‑148 (MQ=255)
                                          cATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGc                                                                                                           >  1:691355/1‑148 (MQ=255)
                                              tCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACgg                                                                                                       >  2:1176015/1‑148 (MQ=255)
                                               cATAGCTTCCCTGGAACTGTTCGGAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACggg                                                                                                      <  1:1316601/148‑1 (MQ=255)
                                               cATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACggg                                                                                                      >  2:2701502/1‑148 (MQ=255)
                                                aTAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGt                                                                                                     <  1:2199441/148‑1 (MQ=255)
                                                aTAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGt                                                                                                     >  1:193226/1‑148 (MQ=255)
                                                 tAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTg                                                                                                    <  1:749638/148‑1 (MQ=255)
                                                 tAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTg                                                                                                    >  1:1458023/1‑148 (MQ=255)
                                                  aGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGa                                                                                                   >  1:3042602/1‑148 (MQ=255)
                                                   gCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGAc                                                                                                  <  2:755162/148‑1 (MQ=255)
                                                      tCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGt                                                                                               >  1:2282578/1‑148 (MQ=255)
                                                      tCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGt                                                                                               <  1:939106/148‑1 (MQ=255)
                                                        ccTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGtat                                                                                             <  2:432028/148‑1 (MQ=255)
                                                        ccTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGtat                                                                                             <  1:1989910/148‑1 (MQ=255)
                                                         cTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGtata                                                                                            <  1:3068687/148‑1 (MQ=255)
                                                           ggAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTTCTTTCTGCACGGGTGACGGTATAGa                                                                                          <  1:2834877/148‑1 (MQ=255)
                                                             aaCTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGGGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGaac                                                                                        >  2:593350/1‑148 (MQ=255)
                                                              aCTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGaaca                                                                                       <  2:359783/148‑1 (MQ=255)
                                                               cTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGaacaa                                                                                      <  2:1010429/148‑1 (MQ=255)
                                                                  ttCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAAccc                                                                                   <  1:1799880/148‑1 (MQ=255)
                                                                  ttCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAAAGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAAccc                                                                                   >  1:2251985/1‑148 (MQ=255)
                                                                        ggCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTTCATAACGCTTCCATATTCAACGAACGCAGAGGATCAAACTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGc                                                                             >  1:3067206/1‑148 (MQ=255)
                                                                            gTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCa                                                                         >  1:2858728/1‑148 (MQ=255)
                                                                              gACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATa                                                                       <  1:2076737/148‑1 (MQ=255)
                                                                               aCGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATaa                                                                      <  2:2458936/148‑1 (MQ=255)
                                                                                cGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAg                                                                     <  1:187636/148‑1 (MQ=255)
                                                                                  tCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAgcg                                                                   >  1:2639150/1‑148 (MQ=255)
                                                                                   cGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAgcgc                                                                  >  2:2038544/1‑148 (MQ=255)
                                                                                     aTTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCa                                                                <  1:457006/148‑1 (MQ=255)
                                                                                      ttttttGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCAt                                                               <  1:3029819/148‑1 (MQ=255)
                                                                                       tttttGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCAtt                                                              >  2:5114/1‑148 (MQ=255)
                                                                                       tttttGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCAtt                                                              >  1:298827/1‑148 (MQ=255)
                                                                                         tttGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCGTCCATATTAAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGc                                                            <  2:1257749/148‑1 (MQ=255)
                                                                                          ttGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGcc                                                           <  2:30552/148‑1 (MQ=255)
                                                                                           tGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCg                                                          >  1:2023252/1‑148 (MQ=255)
                                                                                           tGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCg                                                          <  1:2562325/148‑1 (MQ=255)
                                                                                            gCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCgg                                                         >  2:387053/1‑148 (MQ=255)
                                                                                              tCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGgaa                                                       <  1:2394516/148‑1 (MQ=255)
                                                                                              tCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTAAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGgaa                                                       >  2:706120/1‑148 (MQ=255)
                                                                                               cGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGgaaa                                                      >  2:2217861/1‑148 (MQ=255)
                                                                                               cGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGgaaa                                                      >  2:2218136/1‑148 (MQ=255)
                                                                                                   ggTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGgaaagaaa                                                  >  1:1560992/1‑148 (MQ=255)
                                                                                                         ataACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTtc                                            <  2:3086269/148‑1 (MQ=255)
                                                                                                          taACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTtca                                           <  1:538078/148‑1 (MQ=255)
                                                                                                          taACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTtca                                           <  2:2657640/148‑1 (MQ=255)
                                                                                                              gTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTtcatcat                                       <  1:3176258/148‑1 (MQ=255)
                                                                                                                  tGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCAtt                                   <  2:1362013/148‑1 (MQ=255)
                                                                                                                    ccgccTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTaa                                 <  2:1995869/148‑1 (MQ=255)
                                                                                                                      gccTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAAt                               <  2:2121462/148‑1 (MQ=255)
                                                                                                                          gATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTAcc                           >  1:1708359/1‑148 (MQ=255)
                                                                                                                           aTCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCg                          >  1:1195532/1‑148 (MQ=255)
                                                                                                                           aTCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCg                          >  1:1036149/1‑148 (MQ=255)
                                                                                                                            tCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGt                         <  2:59666/148‑1 (MQ=255)
                                                                                                                            tCCTTGCATAACGCTTCCATAGTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGt                         <  2:1733730/148‑1 (MQ=255)
                                                                                                                             ccTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTc                        >  2:2666614/1‑148 (MQ=255)
                                                                                                                             ccTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTc                        >  2:3128836/1‑148 (MQ=255)
                                                                                                                                tGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTgg                     >  1:168038/1‑148 (MQ=255)
                                                                                                                                tGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTgg                     <  1:2935163/148‑1 (MQ=255)
                                                                                                                                 gCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGt                    <  1:891724/148‑1 (MQ=255)
                                                                                                                                     aaCGCTTCCATATTAAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGc                <  1:1146089/148‑1 (MQ=255)
                                                                                                                                      aCGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCa               >  2:2389397/1‑148 (MQ=255)
                                                                                                                                       cGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATATAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATAACCGTCTGGTTTGCAg              >  2:376218/1‑148 (MQ=255)
                                                                                                                                        gCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAgg             >  2:2614635/1‑148 (MQ=255)
                                                                                                                                            ccATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAg         >  1:3242700/1‑148 (MQ=255)
                                                                                                                                             cATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGt        >  2:677935/1‑148 (MQ=255)
                                                                                                                                               tatTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCt      >  2:1659323/1‑148 (MQ=255)
                                                                                                                                                atTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCtt     <  2:2906402/148‑1 (MQ=255)
                                                                                                                                                   cAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCttatt  <  1:3103658/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAAGAACAGACAAACTGCTCTTTACGGAAGGTATAACCGCGCATTTCATAGCTTCCCTGGAACTGTTCGAAGGCGGTGACGTCGATTTTTTGCTCGCCGGTGTTATAACGTTTTGCCGCCTGATCCTTGCATAACGCTTCCATATTCAACGAACGCAGAGGATCAACCTTTACCCGTTGTGCTTTCTGCACGGGTGACGGTATAGAACAACCCACTAGCACCATAAGCGCCATTGCCGGAAAGAAAGCGTTCATCATCATTAAATTACCGTCTGGTTTGCAGGTGAGTCTTATT  >  CP009273/3718812‑3719106

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: ‑