Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,684,105 C→A 100% intergenic (+97/‑4) manA → / → ydgA mannose‑6‑phosphate isomerase/DUF945 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,684,1050CA97.7% 672.3 / 5.9 214intergenic (+97/‑4)manA/ydgAmannose‑6‑phosphate isomerase/DUF945 family protein
Reads supporting (aligned to +/- strand):  ref base C (4/1);  new base A (145/64);  total (149/65)
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.86e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 227 total aligned reads displayed.

GACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACA  >  CP009273/1683960‑1684249
                                                                                                                                                 |                                                                                                                                                
 aCTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTTTATGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGTATAATAAtt                                                                                                                                               <  1:1860927/148‑1 (MQ=255)
 aCTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAAtt                                                                                                                                               <  2:453631/148‑1 (MQ=255)
   tGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTAt                                                                                                                                             <  2:1909890/148‑1 (MQ=255)
       aaaGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCACTCCGCCTGATTTCATTTTATGGATAATAATTATGAAt                                                                                                                                         >  2:3169014/1‑148 (MQ=255)
         aGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATaa                                                                                                                                       <  2:3007746/148‑1 (MQ=255)
         aGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCATCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATaa                                                                                                                                       >  2:1351244/1‑148 (MQ=255)
           gCCACGGCCGTGTAGCCCGTGTTTAAAACAAGCTGTAAGAGCTTACTGAAAAAATTAACAACTCTTGCTAACCTTATTAAAGGCTTATAACACCATCAGGCGGCCCGTCCGCCGGATTTCATTTTATGGATAATAAATATGTATAAAt                                                                                                                                     >  1:760504/1‑148 (MQ=255)
              aCGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGTCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGc                                                                                                                                  <  2:764554/148‑1 (MQ=255)
               cGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCt                                                                                                                                 <  1:100570/148‑1 (MQ=255)
                  ccGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGt                                                                                                                              <  1:2678950/148‑1 (MQ=255)
                     tttAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCCTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGc                                                                                                                           >  2:2253628/1‑148 (MQ=255)
                          cgcgTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGaa                                                                                                                                          >  2:649733/1‑128 (MQ=255)
                              tgtTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGt                                                                                                                  >  2:2400812/1‑148 (MQ=255)
                                tttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCa                                                                                                                >  2:950973/1‑148 (MQ=255)
                                 ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                                               >  1:877345/1‑148 (MQ=255)
                                 ttACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAt                                                                                                               >  2:1618444/1‑148 (MQ=255)
                                   acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                                                                                             <  1:954856/148‑1 (MQ=255)
                                   acaacaAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCAttg                                                                                                             <  1:955048/148‑1 (MQ=255)
                                       caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                                                                                                         >  1:3561524/1‑148 (MQ=255)
                                       caaGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTgc                                                                                                         >  2:3232667/1‑148 (MQ=255)
                                           cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTa                                                                                                     >  2:186451/1‑148 (MQ=255)
                                           cTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTa                                                                                                     >  2:886231/1‑148 (MQ=255)
                                             gTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAgg                                                                                                   <  1:2308839/148‑1 (MQ=255)
                                              tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                                                                                                  >  1:3558924/1‑148 (MQ=255)
                                              tAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGc                                                                                                  >  2:2884935/1‑148 (MQ=255)
                                               aaGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                                                                                                 >  2:2234827/1‑148 (MQ=255)
                                               aaGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCg                                                                                                 <  2:1770095/148‑1 (MQ=255)
                                                 gagCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTa                                                                                               >  2:3288661/1‑148 (MQ=255)
                                                    cTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc                                                                                            <  1:1489095/148‑1 (MQ=255)
                                                    cTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc                                                                                            <  1:662104/148‑1 (MQ=255)
                                                    cTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTc                                                                                            >  1:43721/1‑148 (MQ=255)
                                                     ttACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCt                                                                                           <  1:1730949/148‑1 (MQ=255)
                                                      tACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTg                                                                                          >  1:3324037/1‑148 (MQ=255)
                                                       aCTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTgg                                                                                         >  2:3502910/1‑148 (MQ=255)
                                                        cTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGa                                                                                        >  1:2051429/1‑148 (MQ=255)
                                                         tGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGAc                                                                                       >  1:159730/1‑148 (MQ=255)
                                                          gAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCAGTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACa                                                                                      >  1:2704564/1‑148 (MQ=255)
                                                                aTTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAggcggc                                                                                <  2:1046685/148‑1 (MQ=255)
                                                                 ttAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGgcggcg                                                                               >  2:1192214/1‑148 (MQ=255)
                                                                  tAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGgcgc                                                                              <  2:301247/148‑1 (MQ=255)
                                                                    aCATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAt                                                                            >  2:2908352/1‑148 (MQ=255)
                                                                      aTCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATgg                                                                          >  1:1296574/1‑148 (MQ=255)
                                                                        ctctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                                        <  1:2400812/148‑1 (MQ=255)
                                                                        ctctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                                        >  1:765682/1‑148 (MQ=255)
                                                                        ctctTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGta                                                                        >  1:1548356/1‑148 (MQ=255)
                                                                              cTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAgg                                                                  <  2:2227653/148‑1 (MQ=255)
                                                                               tAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGc                                                                 >  1:2494824/1‑148 (MQ=255)
                                                                                  gCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaag                                                              >  1:1815764/1‑148 (MQ=255)
                                                                                     tattaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCaagaag                                                           <  2:936564/148‑1 (MQ=255)
                                                                                       ttaAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAt                                                         >  2:180606/1‑148 (MQ=255)
                                                                                        taAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                        >  1:2798017/1‑148 (MQ=255)
                                                                                        taAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGAtt                                                        >  2:2943670/1‑148 (MQ=255)
                                                                                         aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTg                                                       >  1:999592/1‑148 (MQ=255)
                                                                                         aaaGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTg                                                       >  1:999813/1‑148 (MQ=255)
                                                                                              cTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAAcc                                                  <  2:1495051/148‑1 (MQ=255)
                                                                                                tataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCa                                                >  2:2056310/1‑148 (MQ=255)
                                                                                                 ataACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAt                                               <  2:1888865/148‑1 (MQ=255)
                                                                                                  taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                                              >  2:534324/1‑148 (MQ=255)
                                                                                                  taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                                              >  2:2606104/1‑148 (MQ=255)
                                                                                                  taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtc                                              <  1:419437/148‑1 (MQ=255)
                                                                                                  taACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATAGAAACCCAtc                                              >  2:2301916/1‑148 (MQ=255)
                                                                                                   aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                             <  2:1130454/148‑1 (MQ=255)
                                                                                                   aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                             >  2:1847238/1‑148 (MQ=255)
                                                                                                   aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                             >  2:1846978/1‑148 (MQ=255)
                                                                                                   aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                             <  1:1150908/148‑1 (MQ=255)
                                                                                                   aaCACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtct                                             >  1:1312825/1‑148 (MQ=255)
                                                                                                    acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtctc                                            >  2:2543522/1‑148 (MQ=255)
                                                                                                    acacCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCAtctc                                            >  2:755110/1‑148 (MQ=255)
                                                                                                      acCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTCGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGa                                          >  2:618927/1‑148 (MQ=255)
                                                                                                        cTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTAGGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAg                                        <  2:2322964/148‑1 (MQ=255)
                                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGAc                                      <  1:985359/148‑1 (MQ=255)
                                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGAc                                      >  2:353583/1‑148 (MQ=255)
                                                                                                          tCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCAGGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGAc                                      >  1:1357551/1‑148 (MQ=255)
                                                                                                             ggcggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATg                                   >  2:1944655/1‑148 (MQ=255)
                                                                                                              gcggcCAGTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATgg                                  >  1:3500422/1‑148 (MQ=255)
                                                                                                                     gTCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCa                           >  1:337933/1‑148 (MQ=255)
                                                                                                                      tCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAg                          >  1:1098428/1‑148 (MQ=255)
                                                                                                                      tCCGCCTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAg                          >  2:875537/1‑148 (MQ=255)
                                                                                                                          ccTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTCGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGa                      >  2:3115676/1‑148 (MQ=255)
                                                                                                                           cTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGaa                     >  1:2271138/1‑148 (MQ=255)
                                                                                                                           cTGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGaa                     >  2:2661561/1‑148 (MQ=255)
                                                                                                                            tGATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAAc                    >  1:612679/1‑148 (MQ=255)
                                                                                                                             gATTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAAcg                   >  2:1007291/1‑148 (MQ=255)
                                                                                                                              aTTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAAcgc                  >  1:599431/1‑148 (MQ=255)
                                                                                                                              aTTTCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGAGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAAcgc                  >  1:566567/1‑148 (MQ=255)
                                                                                                                                 tCATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCa               >  1:9007/1‑148 (MQ=255)
                                                                                                                                  cATTTTATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCaa              <  1:2839456/148‑1 (MQ=255)
                                                                                                                                    ttttATGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACt            >  1:2163717/1‑148 (MQ=255)
                                                                                                                                        aTGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCaaa        >  2:770646/1‑148 (MQ=255)
                                                                                                                                        aTGGATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCaaa        >  1:1960176/1‑148 (MQ=255)
                                                                                                                                          ggATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACt      >  1:226524/1‑148 (MQ=255)
                                                                                                                                          ggATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACt      <  2:2798017/148‑1 (MQ=255)
                                                                                                                                          ggATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACt      >  1:226650/1‑148 (MQ=255)
                                                                                                                                           gATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTg     >  2:2702801/1‑148 (MQ=255)
                                                                                                                                           gATAATAATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTg     <  1:2606104/148‑1 (MQ=255)
                                                                                                                                            ataataATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGa    >  2:3480154/1‑148 (MQ=255)
                                                                                                                                             taataaTTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAc   >  1:3250833/1‑148 (MQ=255)
                                                                                                                                             taataaTTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAc   >  2:3338474/1‑148 (MQ=255)
                                                                                                                                             taataaTTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGAc   >  1:3007442/1‑148 (MQ=255)
                                                                                                                                              aataatTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACa  >  2:3056750/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GACTGTCAAAGGCCACGGCCGTTTAGCGCGTGTTTACAACAAGCTGTAAGAGCTTACTGAAAAAATTAACATCTCTTGCTAAGCTTATTAAAGGCTTATAACACCTTCAGGCGGCCAGTCCGCCTGATTTCATTTTATGGATAATCATTATGAATAAATCGCTGGTAGCGGTAGGCGTCATTGTTGCGCTAGGCGTAGTCTGGACAGGCGGCGCATGGTATACAGGCAAGAAGATTGAAACCCATCTCGAAGACATGGTCGCGCAGGCGAACGCGCAACTCAAACTGACA  >  CP009273/1683960‑1684249

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: