Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 525,424 T→C 9.5% pseudogene (49/357 nt) ylbG ← pseudogene, DNA‑binding transcriptional regulator homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273525,4240TC 9.5% 828.9 / 62.2 264pseudogene (49/357 nt)ylbGpseudogene, DNA‑binding transcriptional regulator homology
Reads supporting (aligned to +/- strand):  ref base T (125/114);  new base C (16/9);  total (141/123)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

Only 100 of 269 total aligned reads displayed.

TTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTCGCGACTAATA  >  CP009273/525276‑525572
                                                                                                                                                    |                                                                                                                                                    
ttCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGAtt                                                                                                                                                      >  1:4505022/1‑149 (MQ=255)
      aGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCaa                                                                                                                                                >  1:2544682/1‑149 (MQ=255)
       ggTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCaaa                                                                                                                                               >  2:3295240/1‑149 (MQ=255)
               gcgcGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCttt                                                                                                                                       <  2:3383068/149‑1 (MQ=255)
                 gcggcgACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTtgt                                                                                                                                     <  1:3962608/149‑1 (MQ=255)
                 gcggcgACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTtgt                                                                                                                                     <  2:2939224/149‑1 (MQ=255)
                   ggcgACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTgtgt                                                                                                                                   >  2:3341032/1‑149 (MQ=255)
                    gcgACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGtt                                                                                                                                  <  2:1579540/149‑1 (MQ=255)
                     cgACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGttt                                                                                                                                 <  1:1376880/149‑1 (MQ=255)
                      gACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTtg                                                                                                                                >  1:1160348/1‑149 (MQ=255)
                       aCTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTtga                                                                                                                               >  1:1302548/1‑149 (MQ=255)
                       aCTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTtga                                                                                                                               <  2:1989500/149‑1 (MQ=255)
                       aCTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTtga                                                                                                                               <  1:2562350/149‑1 (MQ=255)
                        cTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTtgat                                                                                                                              >  2:1767514/1‑149 (MQ=255)
                         tACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTtgatg                                                                                                                             >  1:4416421/1‑149 (MQ=255)
                             aTTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACgg                                                                                                                         >  2:805421/1‑149 (MQ=255)
                                    tGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGc                                                                                                                  <  1:2708326/149‑1 (MQ=255)
                                    tGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGc                                                                                                                  >  2:443376/1‑149 (MQ=255)
                                       gcATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGt                                                                                                               <  2:151817/149‑1 (MQ=255)
                                         aTCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTAt                                                                                                             <  2:4013861/149‑1 (MQ=255)
                                         aTCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATAAGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGTCTTTGTGTTTGATGACGGGATTGGCAGTAt                                                                                                             >  2:1929309/1‑149 (MQ=255)
                                             acacCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAg                                                                                                         >  1:3262239/1‑149 (MQ=255)
                                              cacCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGc                                                                                                        <  1:3469785/149‑1 (MQ=255)
                                                  gCCTTCTGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGGTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATga                                                                                                    <  1:683622/149‑1 (MQ=255)
                                                  gCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATga                                                                                                    >  1:3689352/1‑149 (MQ=255)
                                                  gCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATga                                                                                                    >  2:312892/1‑149 (MQ=255)
                                                   ccTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATgag                                                                                                   <  2:73783/149‑1 (MQ=255)
                                                    cTTCAGCGACCAGTTCGCTGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATgaga                                                                                                  <  2:2655003/149‑1 (MQ=255)
                                                    cTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATgaga                                                                                                  >  1:3876045/1‑149 (MQ=255)
                                                      tCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGt                                                                                                <  2:2283732/149‑1 (MQ=255)
                                                      tCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGt                                                                                                <  1:493307/149‑1 (MQ=255)
                                                      tCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGt                                                                                                <  1:324050/149‑1 (MQ=255)
                                                         gCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTAc                                                                                             <  1:4327528/149‑1 (MQ=255)
                                                            aCCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACctc                                                                                          >  2:3278850/1‑149 (MQ=255)
                                                            aCCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACctc                                                                                          >  2:121359/1‑149 (MQ=255)
                                                            aCCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACctc                                                                                          <  2:4238363/149‑1 (MQ=255)
                                                             ccAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACctct                                                                                         <  1:1004280/149‑1 (MQ=255)
                                                               aGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTg                                                                                       >  1:2387533/1‑149 (MQ=255)
                                                                gTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGt                                                                                      >  1:2770786/1‑149 (MQ=255)
                                                                gTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGt                                                                                      <  2:4581665/149‑1 (MQ=255)
                                                                 ttCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGtt                                                                                     <  2:4113546/149‑1 (MQ=255)
                                                                     cgGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTgg                                                                                 >  2:4543711/1‑149 (MQ=255)
                                                                     cgGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTgg                                                                                 <  1:4645723/149‑1 (MQ=255)
                                                                        tAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATa                                                                              <  2:2263456/149‑1 (MQ=255)
                                                                         aaCGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATaa                                                                             >  2:3448945/1‑149 (MQ=255)
                                                                          aCGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAg                                                                            >  2:3900698/1‑149 (MQ=255)
                                                                          aCGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAg                                                                            <  1:2364673/149‑1 (MQ=255)
                                                                            gATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGa                                                                          <  2:1307021/149‑1 (MQ=255)
                                                                            gATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCCCTTGTTTTGGATAAGGa                                                                          >  1:3135748/1‑149 (MQ=255)
                                                                              tAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGAtt                                                                        >  2:3603637/1‑149 (MQ=255)
                                                                              tAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGAtt                                                                        <  1:3079005/149‑1 (MQ=255)
                                                                              tAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGAtt                                                                        >  1:1711810/1‑149 (MQ=255)
                                                                              tAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGAtt                                                                        <  2:2039558/149‑1 (MQ=255)
                                                                              tAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGAtt                                                                        <  1:767585/149‑1 (MQ=255)
                                                                               aaaaCGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTc                                                                       <  1:2656689/149‑1 (MQ=255)
                                                                               aaaaCGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTc                                                                       <  1:3928701/149‑1 (MQ=255)
                                                                                 aaCGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGa                                                                     <  2:2580205/149‑1 (MQ=255)
                                                                                   cGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGaca                                                                   >  2:4416106/1‑149 (MQ=255)
                                                                                    gTATCACGCTAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGacac                                                                  <  1:3329463/149‑1 (MQ=255)
                                                                                    gTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGacac                                                                  >  2:1470766/1‑149 (MQ=255)
                                                                                    gTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGacac                                                                  <  2:3544924/149‑1 (MQ=255)
                                                                                     tATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACt                                                                 >  2:2211733/1‑149 (MQ=255)
                                                                                      aTCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTc                                                                >  2:1381144/1‑149 (MQ=255)
                                                                                       tCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCa                                                               <  2:2610191/149‑1 (MQ=255)
                                                                                         aCGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCata                                                             <  1:1746007/149‑1 (MQ=255)
                                                                                         aCGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCata                                                             <  2:3714353/149‑1 (MQ=255)
                                                                                               aaCGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCaaaa                                                       <  2:1120754/149‑1 (MQ=255)
                                                                                                aCGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAAc                                                      <  2:487416/149‑1 (MQ=255)
                                                                                                   cccATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCgg                                                   >  1:1401006/1‑149 (MQ=255)
                                                                                                   cccATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATCGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCgg                                                   >  1:3045148/1‑149 (MQ=255)
                                                                                                     cATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTa                                                 >  1:3616653/1‑149 (MQ=255)
                                                                                                       tGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAAc                                               >  2:2497026/1‑149 (MQ=255)
                                                                                                          ttttACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAActct                                            >  1:2492843/1‑149 (MQ=255)
                                                                                                             tACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCaa                                         <  2:1433360/149‑1 (MQ=255)
                                                                                                             tACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCaa                                         >  2:3106264/1‑149 (MQ=255)
                                                                                                              aCAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAAc                                        >  1:918356/1‑149 (MQ=255)
                                                                                                              aCAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAAc                                        >  1:918347/1‑149 (MQ=255)
                                                                                                                aGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCt                                      <  1:1797359/149‑1 (MQ=255)
                                                                                                                 ggCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCtt                                     >  2:2332824/1‑149 (MQ=255)
                                                                                                                     tttGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCaa                                 <  2:2770786/149‑1 (MQ=255)
                                                                                                                        gATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGc                              >  2:3873030/1‑149 (MQ=255)
                                                                                                                          tACGTTGCTGTGTTCTTCAGCCAGATTGAGCAATCCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGTTAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCa                            <  1:2617220/149‑1 (MQ=255)
                                                                                                                           aCGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCAt                           <  2:2335475/149‑1 (MQ=255)
                                                                                                                             gTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCAtgt                         >  2:616252/1‑149 (MQ=255)
                                                                                                                              ttGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCAtgtg                        >  1:4115983/1‑149 (MQ=255)
                                                                                                                               tGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATgtgt                       <  1:2857973/149‑1 (MQ=255)
                                                                                                                                gCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTc                      <  1:4192336/149‑1 (MQ=255)
                                                                                                                                  tGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAACCCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAg                    >  1:4568364/1‑149 (MQ=255)
                                                                                                                                     gTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGTGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATc                 <  2:2644836/149‑1 (MQ=255)
                                                                                                                                     gTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATc                 >  2:3420421/1‑149 (MQ=255)
                                                                                                                                     gTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATc                 <  1:2125071/149‑1 (MQ=255)
                                                                                                                                         ttcAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGt             <  1:3271109/149‑1 (MQ=255)
                                                                                                                                          tcAGCCAGATTGAGCAAACCGGCGTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTc            <  2:3190419/149‑1 (MQ=255)
                                                                                                                                            aGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTcgc          >  1:2184861/1‑149 (MQ=255)
                                                                                                                                               cAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATAAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTCGCGAc       <  2:1679825/149‑1 (MQ=255)
                                                                                                                                                 gATCGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGAAGTATGAAGCATGAGAGAAAACTCTTGTTTTGGATAAGGATTCGACACTCAAATCAAAACCGGTAACTCTCAAACTTTCAATGCCATGTGTCAGATCAATTAGAGACTa     >  2:1045965/1‑149 (MQ=255)
                                                                                                                                                  aTTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTCGCGACTaa    <  1:2229589/149‑1 (MQ=255)
                                                                                                                                                  aTTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTCGCGACTaa    <  1:1018084/149‑1 (MQ=255)
                                                                                                                                                   ttGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTCGCGACTAAt   >  2:2819601/1‑149 (MQ=255)
                                                                                                                                                    tGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTAGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCATACAAGGCCAAGTGTCAGATCAAGACGCGACTAATa  >  2:1182020/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCTTAAGGTTAGGAGCGCGGCGACTACGATTAATCTGCGCATCCACACCGCCTTCAGCGACCAGTTCGCGGTAACGATAAAACGTATCACGCGAAACGCCCATGATTTTACAGGCTTTTGATACGTTGCTGAGTTCTTCAGCCAGATTGAGCAAACCGGCTTTGTGTTTGATGACGGGATTGGCAGTATGAAGCATGAGAGTTACCTCTTGTTTTGGATAAGGATTCGACACTCATATCAAAACCGGTAACTCTCAACCTTTCAAGGCCATGTGTCAGATCAAGTCGCGACTAATA  >  CP009273/525276‑525572

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: