Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,211,803 G→A 100% P205S (CCT→TCT)  arpA ← ankyrin repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,211,8030GA99.6% 782.8 / ‑4.5 237P205S (CCT→TCT) arpAankyrin repeat protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (134/102);  total (134/103)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 237 total aligned reads displayed.

TTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTG  >  CP009273/4211658‑4211947
                                                                                                                                                 |                                                                                                                                                
tttATTCTTAGCCGGACCGTCAAAGCATGGTTGTCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAAc                                                                                                                                                >  2:195434/1‑148 (MQ=255)
tttATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAAc                                                                                                                                                >  2:194888/1‑148 (MQ=255)
     tCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTg                                                                                                                                           >  1:2473653/1‑148 (MQ=255)
     tCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTg                                                                                                                                           >  2:360098/1‑148 (MQ=255)
      cTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTgg                                                                                                                                          <  1:823260/148‑1 (MQ=255)
       ttAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGa                                                                                                                                         >  1:3059292/1‑148 (MQ=255)
         aGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGaaa                                                                                                                                       <  1:2158579/148‑1 (MQ=255)
             ggACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCaa                                                                                                                                   <  1:1550153/148‑1 (MQ=255)
                       aGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTc                                                                                                                         >  1:1091137/1‑148 (MQ=255)
                        gCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCt                                                                                                                        >  1:736209/1‑148 (MQ=255)
                         cATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTg                                                                                                                       <  1:2878439/148‑1 (MQ=255)
                          aTGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTgg                                                                                                                      >  1:3420001/1‑148 (MQ=255)
                          aTGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTgg                                                                                                                      <  2:3176091/148‑1 (MQ=255)
                           tGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGt                                                                                                                     >  2:1102339/1‑148 (MQ=255)
                           tGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGt                                                                                                                     >  2:1265351/1‑148 (MQ=255)
                            ggTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtt                                                                                                                    >  2:3264690/1‑148 (MQ=255)
                              ttATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtttt                                                                                                                  <  2:1819461/148‑1 (MQ=255)
                              ttATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGtttt                                                                                                                  >  2:528532/1‑148 (MQ=255)
                                aTCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTg                                                                                                                <  1:3335327/148‑1 (MQ=255)
                                 tCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGc                                                                                                               >  2:681465/1‑148 (MQ=255)
                                 tCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGc                                                                                                               <  2:2166295/148‑1 (MQ=255)
                                   caccaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAt                                                                                                             <  2:2293683/148‑1 (MQ=255)
                                    accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtt                                                                                                            >  1:3101004/1‑148 (MQ=255)
                                    accaGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtt                                                                                                            <  1:3453176/148‑1 (MQ=255)
                                     ccaGTAATATGGCTGCCATAACACTGATGCCATTAGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttt                                                                                                           <  2:1064402/148‑1 (MQ=255)
                                       aGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAttttt                                                                                                         <  2:2800352/148‑1 (MQ=255)
                                        gTAATATGGCTGCCATAACACTGATGCCTTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCAtttttt                                                                                                        <  1:843617/148‑1 (MQ=255)
                                         tAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCGCCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                       >  2:2143155/1‑148 (MQ=255)
                                         tAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                       <  2:852526/148‑1 (MQ=255)
                                         tAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                       >  1:1792153/1‑148 (MQ=255)
                                         tAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttg                                                                                                       >  1:1794740/1‑148 (MQ=255)
                                          aaTATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgt                                                                                                      >  2:1403971/1‑148 (MQ=255)
                                            tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgttg                                                                                                    >  1:2948884/1‑148 (MQ=255)
                                            tatGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTttgttg                                                                                                    >  2:2324647/1‑148 (MQ=255)
                                               ggCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCCGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                                                                                                 >  2:2043960/1‑148 (MQ=255)
                                               ggCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAg                                                                                                 >  2:1569749/1‑148 (MQ=255)
                                                gCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGt                                                                                                <  2:388398/148‑1 (MQ=255)
                                                   gCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCg                                                                                             <  1:2731654/148‑1 (MQ=255)
                                                       tAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTgag                                                                                         >  2:547554/1‑148 (MQ=255)
                                                            cTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGACCGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAtt                                                                                    >  2:2543631/1‑148 (MQ=255)
                                                             tGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                   <  2:1520499/148‑1 (MQ=255)
                                                             tGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAAttt                                                                                   <  1:1656504/148‑1 (MQ=255)
                                                                  ccATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAg                                                                              <  1:487213/148‑1 (MQ=255)
                                                                   cATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGc                                                                             <  1:2843824/148‑1 (MQ=255)
                                                                    aTTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCt                                                                            >  1:19246/1‑148 (MQ=255)
                                                                      ttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGa                                                                          >  2:1527940/1‑148 (MQ=255)
                                                                      ttGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGa                                                                          >  1:1708048/1‑148 (MQ=255)
                                                                       tGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGaa                                                                         >  2:1691365/1‑148 (MQ=255)
                                                                        gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAag                                                                        >  2:3026485/1‑148 (MQ=255)
                                                                        gATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAag                                                                        >  2:1038957/1‑148 (MQ=255)
                                                                         aTTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAaga                                                                       <  2:736209/148‑1 (MQ=255)
                                                                          ttCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAagag                                                                      <  2:1736458/148‑1 (MQ=255)
                                                                           tCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGt                                                                     <  1:2513305/148‑1 (MQ=255)
                                                                             cATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTaa                                                                   <  2:753557/148‑1 (MQ=255)
                                                                              aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTTat                                                                  >  2:2558590/1‑148 (MQ=255)
                                                                              aTCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAat                                                                  <  1:1964589/148‑1 (MQ=255)
                                                                               tCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatc                                                                 <  1:2043572/148‑1 (MQ=255)
                                                                                  tGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAatcatc                                                              <  1:229564/148‑1 (MQ=255)
                                                                                    tCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCtt                                                            >  1:3463013/1‑148 (MQ=255)
                                                                                      aaaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGa                                                          <  1:299645/148‑1 (MQ=255)
                                                                                       aaaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGTAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGAt                                                         >  2:2698181/1‑148 (MQ=255)
                                                                                        aaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATg                                                        <  1:3257366/148‑1 (MQ=255)
                                                                                        aaTATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATg                                                        >  1:643058/1‑148 (MQ=255)
                                                                                          tatGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATggg                                                      <  2:2073243/148‑1 (MQ=255)
                                                                                           atGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATgggg                                                     <  2:2420808/148‑1 (MQ=255)
                                                                                             gCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGtt                                                   <  2:3059292/148‑1 (MQ=255)
                                                                                             gCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGtt                                                   <  1:1532816/148‑1 (MQ=255)
                                                                                              cTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGttt                                                  >  1:60436/1‑148 (MQ=255)
                                                                                              cTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGttt                                                  >  1:648999/1‑148 (MQ=255)
                                                                                                  cATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTaa                                              <  2:2686623/148‑1 (MQ=255)
                                                                                                    tAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                            >  2:3305599/1‑148 (MQ=255)
                                                                                                    tAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                            >  1:722512/1‑148 (MQ=255)
                                                                                                    tAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                            >  2:892705/1‑148 (MQ=255)
                                                                                                    tAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATc                                            >  1:1033385/1‑148 (MQ=255)
                                                                                                    tAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTAAGCTGAAGAGTAATCATCTGGATGGGGTTTGTAATc                                            >  1:1033400/1‑148 (MQ=255)
                                                                                                      aCTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCaa                                          <  2:639045/148‑1 (MQ=255)
                                                                                                       cTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAAc                                         <  1:169724/148‑1 (MQ=255)
                                                                                                          gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACAtt                                      <  1:1720058/148‑1 (MQ=255)
                                                                                                          gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACAtt                                      >  2:199450/1‑148 (MQ=255)
                                                                                                          gCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACAtt                                      <  2:1210143/148‑1 (MQ=255)
                                                                                                           cccGAAATCACCGCAGGCAGAAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTc                                     <  1:1580169/148‑1 (MQ=255)
                                                                                                             cGAAATCACCTCAGTCATCAACGCAACGCCCCTCAGAACTGCTGTAAACCAAAGACAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTg                                   <  2:2825442/148‑1 (MQ=255)
                                                                                                             cGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTg                                   >  2:829927/1‑148 (MQ=255)
                                                                                                              gAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGa                                  >  1:2650985/1‑148 (MQ=255)
                                                                                                                aaTCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATa                                >  2:2192347/1‑148 (MQ=255)
                                                                                                                aaTCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATa                                >  2:656519/1‑148 (MQ=255)
                                                                                                                  tCACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATaaa                              >  1:3171166/1‑148 (MQ=255)
                                                                                                                   cACCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAat                             >  1:2518538/1‑148 (MQ=255)
                                                                                                                    aCCGCAGGCAGCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAata                            <  2:1342731/148‑1 (MQ=255)
                                                                                                                          ggcagCAACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACACTCTGATAAATATCATTa                      >  1:1932026/1‑148 (MQ=255)
                                                                                                                            cagcaACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAat                    <  2:425063/148‑1 (MQ=255)
                                                                                                                            cagcaACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAat                    <  2:2948758/148‑1 (MQ=255)
                                                                                                                            cagcaACGCAACGCCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAat                    <  2:2487208/148‑1 (MQ=255)
                                                                                                                            cagcaACGCAACGCCCCTCACAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAat                    <  2:2190453/148‑1 (MQ=255)
                                                                                                                                  cgcaacgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATtg              >  2:1449106/1‑148 (MQ=255)
                                                                                                                                    caacgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATtgtg            >  2:2014108/1‑148 (MQ=255)
                                                                                                                                     aacgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAAAATTCTGATAAATATCATTAATATATTGTgc           >  2:1366546/1‑148 (MQ=255)
                                                                                                                                       cgcCCCTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcag         >  2:3327520/1‑148 (MQ=255)
                                                                                                                                         ccccTCAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTgcagca       <  1:1437478/148‑1 (MQ=255)
                                                                                                                                              cAGAACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAATAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAAtgtg  >  2:2135254/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TTTATTCTTAGCCGGACCGTCAAAGCATGGTTATCCACCAGTAATATGGCTGCCATAACACTGATGCCATTTGATTCCATCTTGTCAAAATATGCTTTCATAACTAGCCCGAAATCACCGCAGGCAGCAACGCAACGCCCCTCAGGACGGCTGGAAACCAACGCCAAAGTCTGGTTTTTGCATTTTTTGTTGAAGTGCGTGAGAAATTTTTCAGCTGAAGAGTAATCATCTTGATGGGGTTTGTAATCAACATTCTGATAAATATCATTAATATATTGTGCAGCAATGTG  >  CP009273/4211658‑4211947

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