Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,989,295 C→A 69.9% pseudogene (85/261 nt) ygeN ← pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,989,2950CA69.9% 103.6 / 71.4 76pseudogene (85/261 nt)ygeNpseudogene
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base A (32/21);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTCAACAA  >  CP009273/2989146‑2989442
                                                                                                                                                     |                                                                                                                                                   
cAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAtttc                                                                                                                                                     >  2:22510/1‑150 (MQ=255)
 aaGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAtttct                                                                                                                                                    >  1:283285/1‑150 (MQ=255)
   gggTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTATGAGAAACTAtttcttt                                                                                                                                                  >  2:4778370/1‑150 (MQ=255)
   gggTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTatttattt                                                                                                                                                  <  1:64169/150‑1 (MQ=255)
    ggTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTAtttatttc                                                                                                                                                 >  1:4171299/1‑150 (MQ=255)
     gTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTATGTGAAAGTATTTATTTCa                                                                                                                                                >  1:3649699/1‑150 (MQ=255)
          aaaTTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCa                                                                                                                                           <  2:568501/150‑1 (MQ=255)
          aaaTTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCa                                                                                                                                           <  2:2828232/150‑1 (MQ=255)
            aTTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATc                                                                                                                                         >  2:2163048/1‑150 (MQ=255)
              tATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAg                                                                                                                                       >  2:4676887/1‑150 (MQ=255)
              tATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAg                                                                                                                                       >  2:4641166/1‑150 (MQ=255)
                       catcaATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCATTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTc                                                                                                                              <  1:920168/150‑1 (MQ=255)
                           aataatGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCtttt                                                                                                                          >  2:2573753/1‑150 (MQ=255)
                           aataatGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCtttt                                                                                                                          >  2:2574053/1‑150 (MQ=255)
                                           taGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAttt                                                                                                          <  1:3423837/150‑1 (MQ=255)
                                           taGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCAttt                                                                                                          <  1:3424148/150‑1 (MQ=255)
                                           taGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTAAATACATAAGAAATTTTTCTTTTAATAACATAACAAttt                                                                                                          >  2:1332212/1‑150 (MQ=255)
                                             gCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGc                                                                                                        >  2:2761374/1‑150 (MQ=255)
                                             gCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGc                                                                                                        >  2:2761955/1‑150 (MQ=255)
                                               catcatAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAg                                                                                                      >  2:4009448/1‑150 (MQ=255)
                                                  catAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGgag                                                                                                   <  1:3366629/150‑1 (MQ=255)
                                                  catAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGgag                                                                                                   <  1:1261910/150‑1 (MQ=255)
                                                         acattGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTAt                                                                                            <  2:3336164/147‑1 (MQ=255)
                                                         acattGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTAt                                                                                            <  2:4055117/147‑1 (MQ=255)
                                                         aCCTTTAATAAATTTTTCTTATCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTAt                                                                                            <  2:4055013/150‑1 (MQ=255)
                                                              gAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATAACAACACCATTTGCAGGAGATTTTATCAAGc                                                                                       >  2:1309272/1‑150 (MQ=255)
                                                              gAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGc                                                                                       >  2:921292/1‑150 (MQ=255)
                                                               aaTAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCt                                                                                      <  1:1058752/150‑1 (MQ=255)
                                                                aTAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCa                                                                                    <  2:4829890/150‑1 (MQ=255)
                                                                aTAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCa                                                                                    <  2:4829836/150‑1 (MQ=255)
                                                                aTAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCa                                                                                    <  2:2960128/150‑1 (MQ=255)
                                                                 tAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCa                                                                                    <  2:4713162/150‑1 (MQ=255)
                                                                   aaTTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCAtc                                                                                  <  1:3019663/150‑1 (MQ=255)
                                                                   aaTTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCAtc                                                                                  <  1:3019650/150‑1 (MQ=255)
                                                                   aaTTTTTCTTCTCTCTATTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTATAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCAtc                                                                                  <  1:3019816/150‑1 (MQ=255)
                                                                    aTTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACAATTTGCAGGAGATTTTATCAAGCTCAtct                                                                                 >  2:2110630/1‑150 (MQ=255)
                                                                      ttttCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATctct                                                                               >  1:608923/1‑150 (MQ=255)
                                                                      ttttCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATctct                                                                               >  2:1957816/1‑150 (MQ=255)
                                                                        ttcttcTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAt                                                                             >  1:497990/1‑150 (MQ=255)
                                                                        atctactCACTTTTCCTCAACACTATGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAt                                                                             <  1:1332212/149‑1 (MQ=255)
                                                                            tctctcTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAttttt                                                                         <  2:4076078/150‑1 (MQ=255)
                                                                            tctctcTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAttttt                                                                         <  2:4075831/150‑1 (MQ=255)
                                                                            tctctcTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAttttt                                                                         <  2:807366/150‑1 (MQ=255)
                                                                             ctctctTTTCCTCACTACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAtttttt                                                                        <  1:3066517/150‑1 (MQ=255)
                                                                             ctctctTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTAtttttt                                                                        >  1:4090621/1‑150 (MQ=255)
                                                                                tctTTTCCTCAATACTCTGAATTTTTGCGGTTGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTc                                                                     <  2:585578/150‑1 (MQ=255)
                                                                                         cAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAAt                                                            <  1:4578210/150‑1 (MQ=255)
                                                                                         cAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTGTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAAt                                                            >  1:2656092/1‑150 (MQ=255)
                                                                                          aaTACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAAtt                                                           >  1:2722955/1‑149 (MQ=255)
                                                                                             aCTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTAGCATAAAATGTaa                                                        >  2:2502916/1‑150 (MQ=255)
                                                                                             aCTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATAACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTaa                                                        >  2:2502907/1‑150 (MQ=255)
                                                                                                 tGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATtct                                                    <  2:1797115/150‑1 (MQ=255)
                                                                                                  gAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATtctc                                                   >  2:2760664/1‑150 (MQ=255)
                                                                                                  gAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATtctc                                                   >  2:2760729/1‑150 (MQ=255)
                                                                                                    aTTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATtctctc                                                 >  2:38632/1‑150 (MQ=255)
                                                                                                           cGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAgt                                          >  2:2694796/1‑150 (MQ=255)
                                                                                                            ggtggtTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAgtg                                         >  1:4050259/1‑150 (MQ=255)
                                                                                                                    cAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCa                                 >  1:3356384/1‑150 (MQ=255)
                                                                                                                    cAAGAAGGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCa                                 >  1:2523049/1‑150 (MQ=255)
                                                                                                                     aagaagGATTTTTTTTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTGGGATAAAATTTAATTCTCTGTTAAAGTGTTGTTGGaa                                >  2:2939068/1‑150 (MQ=255)
                                                                                                                      agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAg                               >  1:397983/1‑150 (MQ=255)
                                                                                                                      agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAg                               >  2:483476/1‑150 (MQ=255)
                                                                                                                      agaagGATTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAg                               >  1:1462187/1‑150 (MQ=255)
                                                                                                                          ggATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGt                           >  1:4586387/1‑150 (MQ=255)
                                                                                                                            aTTTTTTCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTaa                         >  2:1403606/1‑150 (MQ=255)
                                                                                                                              tttttCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACg                       <  1:2710178/150‑1 (MQ=255)
                                                                                                                                 ttCTTTAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCg                    >  2:1259837/1‑150 (MQ=255)
                                                                                                                                  tCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGa                   <  1:2573753/150‑1 (MQ=255)
                                                                                                                                  tCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGa                   <  1:2574053/150‑1 (MQ=255)
                                                                                                                                    tttAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGAtt                 <  2:4171299/150‑1 (MQ=255)
                                                                                                                                     ttAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTg                <  1:3639519/150‑1 (MQ=255)
                                                                                                                                      tAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGt               <  1:4009448/150‑1 (MQ=255)
                                                                                                                                      tAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTACTTCTCTCTCAAAGTGTTGATGCAAGTGGTAACGGCGATTGt               >  1:997228/1‑150 (MQ=255)
                                                                                                                                      tAAGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTCATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGt               >  1:3444402/1‑150 (MQ=255)
                                                                                                                                       aaGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTa              >  1:353859/1‑150 (MQ=255)
                                                                                                                                       aaGAGAAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTa              <  1:4641166/150‑1 (MQ=255)
                                                                                                                                           gaAACTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTAtgtg          >  1:145532/1‑150 (MQ=255)
                                                                                                                                            aaaCTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTgg         >  2:4698125/1‑150 (MQ=255)
                                                                                                                                             aaCTATTTATTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGt        >  2:773197/1‑150 (MQ=255)
                                                                                                                                                tatttattTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTcaa     >  2:2351970/1‑150 (MQ=255)
                                                                                                                                                tatttattTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTcaa     >  2:2351967/1‑150 (MQ=255)
                                                                                                                                                tATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTcaa     >  2:537990/1‑150 (MQ=255)
                                                                                                                                                tATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTcaa     >  2:538003/1‑150 (MQ=255)
                                                                                                                                                   ttctttcAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTcaacaa  <  2:2898985/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAGGGTATAAAATTATCCATCTCATCAATAATGCCTTTTGTATAGCCATCATAATAACCTTGAATAAATTTTTCTTCTCTCTTTTCCTCAATACTCTGAATTTTTGCGGTGGTTTCAAGAAGGATTTTTTCTTTAAGAGAAACTATTTCTTTCAATCCATCAGAAATTTTTCTTTTAATCACAACACCATTTGCAGGAGATTTTATCAAGCTCATCTCTATTTTTTTTCGCATAAAATGTAATTCTCTCTAAAAGTGTTGATGCAAGTGGTAACGGCGATTGTATGTGGTCAACAA  >  CP009273/2989146‑2989442

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: