Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,427,647 G→A 100% D42N (GAT→AAT)  ytfI → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,427,6470GA98.1% 819.0 / 4.8 256D42N (GAT→AAT) ytfIuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (4/1);  new base A (132/119);  total (136/120)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 261 total aligned reads displayed.

CCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCC  >  CP009273/4427502‑4427791
                                                                                                                                                 |                                                                                                                                                
  aaataaaataAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTg                                                                                                                                              >  1:1855809/1‑148 (MQ=255)
  aaataaaataAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTg                                                                                                                                              >  1:2080660/1‑148 (MQ=255)
  aaataaaataAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTg                                                                                                                                              <  2:2032608/148‑1 (MQ=255)
   aataaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTgg                                                                                                                                             >  2:850618/1‑148 (MQ=255)
     taaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGAt                                                                                                                                           >  1:2125150/1‑148 (MQ=255)
      aaaataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGAtt                                                                                                                                          <  2:3478124/148‑1 (MQ=255)
        aataaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTc                                                                                                                                        >  2:48676/1‑148 (MQ=255)
          taaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGa                                                                                                                                      <  1:1295235/148‑1 (MQ=255)
           aaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGaa                                                                                                                                     >  2:2673729/1‑148 (MQ=255)
           aaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGaa                                                                                                                                     >  2:3153596/1‑148 (MQ=255)
           aaGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGaa                                                                                                                                     >  2:3153565/1‑148 (MQ=255)
              gTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGACTTa                                                                                                                                  >  1:1817107/1‑148 (MQ=255)
                 gTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGa                                                                                                                               >  1:445357/1‑148 (MQ=255)
                  tCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGAt                                                                                                                              >  2:1231327/1‑148 (MQ=255)
                   cATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATa                                                                                                                             <  2:127890/148‑1 (MQ=255)
                      atGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGcc                                                                                                                          >  2:9172/1‑148 (MQ=255)
                      atGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGcc                                                                                                                          >  2:15734/1‑148 (MQ=255)
                        gTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTg                                                                                                                        <  2:2577454/148‑1 (MQ=255)
                         ttACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGt                                                                                                                       <  2:1140221/148‑1 (MQ=255)
                                 gATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTg                                                                                                               >  1:2977844/1‑148 (MQ=255)
                                  aTCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGa                                                                                                              >  1:3355758/1‑148 (MQ=255)
                                   tCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAg                                                                                                             <  1:1643321/148‑1 (MQ=255)
                                        cacaATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGt                                                                                                        >  2:1760253/1‑148 (MQ=255)
                                         acaATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGtt                                                                                                       <  1:2013858/148‑1 (MQ=255)
                                          caATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTc                                                                                                      >  1:2812365/1‑148 (MQ=255)
                                           aataatTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCt                                                                                                     <  1:510050/148‑1 (MQ=255)
                                            ataatTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCtt                                                                                                    <  1:3168853/148‑1 (MQ=255)
                                             taatTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCttt                                                                                                   <  2:395398/148‑1 (MQ=255)
                                                 tttATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGtt                                                                                               >  2:744175/1‑148 (MQ=255)
                                                  ttATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTg                                                                                              <  2:330149/148‑1 (MQ=255)
                                                       ggTGCGGTGGAGTTATTTGTAAAGTATTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCaa                                                                                         <  1:2712014/148‑1 (MQ=255)
                                                        gTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCTAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAAc                                                                                        >  1:735072/1‑148 (MQ=255)
                                                         tGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACt                                                                                       >  2:642501/1‑148 (MQ=255)
                                                         tGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACt                                                                                       >  2:3104307/1‑148 (MQ=255)
                                                          gCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTc                                                                                      <  2:71379/148‑1 (MQ=255)
                                                              tGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGAt                                                                                  >  1:3226465/1‑148 (MQ=255)
                                                                gAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGTTAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTg                                                                                >  2:679037/1‑148 (MQ=255)
                                                                gAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTg                                                                                <  1:1007018/148‑1 (MQ=255)
                                                                gAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTg                                                                                <  1:1006166/148‑1 (MQ=255)
                                                                 aGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGa                                                                               >  2:1423724/1‑148 (MQ=255)
                                                                   ttATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATa                                                                             >  1:1725352/1‑148 (MQ=255)
                                                                       ttGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTTATAAAAg                                                                         <  2:492744/148‑1 (MQ=255)
                                                                          tAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTgg                                                                      >  1:738071/1‑148 (MQ=255)
                                                                            aaGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGtt                                                                    >  1:1503029/1‑148 (MQ=255)
                                                                            aaGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGtt                                                                    >  2:1370260/1‑148 (MQ=255)
                                                                             aGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAACTGGTTa                                                                   >  1:3263186/1‑148 (MQ=255)
                                                                              gTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTAt                                                                  <  1:2046962/148‑1 (MQ=255)
                                                                               tCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATc                                                                 <  1:1585880/148‑1 (MQ=255)
                                                                                cTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCa                                                                >  1:2567430/1‑148 (MQ=255)
                                                                                 ttCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAg                                                               <  2:688075/148‑1 (MQ=255)
                                                                                    gTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTg                                                            <  1:2994468/148‑1 (MQ=255)
                                                                                         cAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtc                                                       >  1:2493020/1‑148 (MQ=255)
                                                                                          aaaaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtca                                                      >  2:2397868/1‑148 (MQ=255)
                                                                                          aaaaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtca                                                      >  1:796513/1‑148 (MQ=255)
                                                                                          aaaaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtca                                                      >  1:650264/1‑148 (MQ=255)
                                                                                           aaaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtcat                                                     >  1:3441656/1‑148 (MQ=255)
                                                                                           aaaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtcat                                                     >  1:564550/1‑148 (MQ=255)
                                                                                            aaaCACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTtcatc                                                    >  1:3498232/1‑148 (MQ=255)
                                                                                             aaaacaTTCAAACAATTTTTTTAACAACATACATCATTCAATTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTTTTAACATTGAGGGATTTCTTTCGTTGCGCAACTCAGATTGCTAAAAGTGGTTATCAGGTGCTTtcatca                                                   <  2:1819970/145‑1 (MQ=255)
                                                                                             aaCACATTCAAACAATTTTTTTAAAAACATCCATCATGCATTTAAGAAAAAAAATTGGATTACGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCATACTTTGCGCAACTCAGATTGATAAAAGAGTATATCAGGTGCTTtcatca                                                   >  1:1368350/1‑148 (MQ=255)
                                                                                                 caTTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCaaaaa                                               >  2:2232699/1‑148 (MQ=255)
                                                                                                 caTTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCaaaaa                                               >  2:3261697/1‑148 (MQ=255)
                                                                                                    tCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATg                                            >  2:1722871/1‑148 (MQ=255)
                                                                                                    tCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATg                                            <  2:609971/148‑1 (MQ=255)
                                                                                                     cAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGa                                           >  1:943559/1‑148 (MQ=255)
                                                                                                       aaCAAATTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTCATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGaaa                                         <  1:3256410/148‑1 (MQ=255)
                                                                                                         cAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACt                                       <  1:2195550/148‑1 (MQ=255)
                                                                                                         cAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACt                                       <  1:1720498/148‑1 (MQ=255)
                                                                                                         cAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACt                                       >  2:2344016/1‑148 (MQ=255)
                                                                                                           aTTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGt                                     >  1:2890869/1‑148 (MQ=255)
                                                                                                            tttttttAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGtt                                    <  1:1454749/148‑1 (MQ=255)
                                                                                                             ttttttAACAACATCCATCATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTc                                   <  1:1075955/148‑1 (MQ=255)
                                                                                                                    acaacaTCCATCATGAATTTAAGAAAAAAAATTGTATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCAt                            <  1:2938594/148‑1 (MQ=255)
                                                                                                                       acaTCCATCATGAATTTAAGAAAAAAAATTTGATTTCTAATTATGATATCCTGTTAACTTTGTGTGTGTTCTTTCGTTGCTCAAATCATATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAAATGTTCACGACATGGa                         <  2:2664317/148‑1 (MQ=255)
                                                                                                                         atccatcATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATa                       <  1:15734/148‑1 (MQ=255)
                                                                                                                           ccatcATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAg                     >  2:1912099/1‑148 (MQ=255)
                                                                                                                           ccatcATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAg                     >  2:652538/1‑148 (MQ=255)
                                                                                                                           ccatcATGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAg                     <  2:2139909/148‑1 (MQ=255)
                                                                                                                            catcatGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGt                    <  1:2673219/148‑1 (MQ=255)
                                                                                                                              tcatGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGttt                  >  2:121104/1‑148 (MQ=255)
                                                                                                                              tcatGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGttt                  <  2:3343075/148‑1 (MQ=255)
                                                                                                                               catGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTc                 <  1:776233/148‑1 (MQ=255)
                                                                                                                               catGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTGGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGGTTc                 >  2:2883911/1‑148 (MQ=255)
                                                                                                                                atGCATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCt                <  1:1231327/148‑1 (MQ=255)
                                                                                                                                   cATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTAt             >  1:690725/1‑148 (MQ=255)
                                                                                                                                   cATTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTAt             <  2:738071/148‑1 (MQ=255)
                                                                                                                                    aTTTAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATTAAACTGTTAACGCCATGGATAAGTTTCTTATc            >  2:112667/1‑148 (MQ=255)
                                                                                                                                     tttAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCa           <  1:927840/148‑1 (MQ=255)
                                                                                                                                      ttAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAg          <  1:1163864/148‑1 (MQ=255)
                                                                                                                                      ttAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAg          >  1:3077486/1‑148 (MQ=255)
                                                                                                                                      ttAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAg          <  2:1857197/148‑1 (MQ=255)
                                                                                                                                      ttAAGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGATCACGCCATGGATAAGTTTCTTATCAg          >  2:1658958/1‑148 (MQ=255)
                                                                                                                                         aGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGttt       >  1:1106513/1‑148 (MQ=255)
                                                                                                                                         aGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGttt       >  2:2170369/1‑148 (MQ=255)
                                                                                                                                         aGAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGttt       <  1:274088/148‑1 (MQ=255)
                                                                                                                                          gAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAAc                                        <  2:3478856/114‑1 (MQ=255)
                                                                                                                                          gAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAAc                                        >  1:3478856/1‑114 (MQ=255)
                                                                                                                                          gAAAAAAAATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTc      <  2:872610/148‑1 (MQ=255)
                                                                                                                                           aaaaaaaaTTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCa     >  2:2528446/1‑148 (MQ=255)
                                                                                                                                              aaaaATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGcc  <  1:850581/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CCAAATAAAATAAGGTTGTCATATGTTACCCAGGATCAGACACAATAATTTTATTGGTGCGGTGGAGTTATTTGTAAAGTCTTCGTATACAAAAACACATTCAAACAATTTTTTTAACAACATCCATCATGCATTTAAGAAAAAAGATTGGATTTCGAATTATGATAGCCTGTTAACTTTGAGGGAGTTCTTTCGTTGCGCAACTCAGATTGATAAAAGTGGTTATCAGGTGCTTTCATCAAAAAATGAAACTGTTCACGCCATGGATAAGTTTCTTATCAGTTTCAGCC  >  CP009273/4427502‑4427791

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: