Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 524,334:1 +T 100% pseudogene (239/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,3341.T95.3% 445.7 / 15.2 127pseudogene (239/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base . (4/2);  new base T (73/48);  total (77/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

Only 100 of 149 total aligned reads displayed.

TATTTATACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATA  >  CP009273/524186‑524478
                                                                                                                                                     |                                                                                                                                                
tatttatACAAAGATGAATCTCTTGTAAGTAAACTTTACATAAACTATATTGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGattat                                                                                                                                                 <  2:86410/150‑1 (MQ=255)
tatttatACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGattat                                                                                                                                                 <  1:2641438/150‑1 (MQ=255)
tatttatACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGattat                                                                                                                                                 >  2:1153005/1‑150 (MQ=255)
 atttatACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGACGATGattatt                                                                                                                                                >  1:3684446/1‑150 (MQ=255)
 atttatACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGattatt                                                                                                                                                >  2:1305183/1‑150 (MQ=255)
   ttatACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAtttt                                                                                                                                              >  2:1681452/1‑150 (MQ=255)
   ttatACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAtttt                                                                                                                                              <  1:3125656/150‑1 (MQ=255)
    tataCAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttt                                                                                                                                             <  2:1193932/150‑1 (MQ=255)
           aGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttttttctta                                                                                                                                      >  2:469035/1‑145 (MQ=255)
           aGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttttttctta                                                                                                                                      >  2:469006/1‑145 (MQ=255)
           aGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttttttctta                                                                                                                                      >  1:1777799/1‑145 (MQ=255)
           aGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttttttctta                                                                                                                                      >  1:1778254/1‑145 (MQ=255)
           aGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttttttctta                                                                                                                                      <  1:3763261/150‑6 (MQ=255)
           aGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAttttttttctta                                                                                                                                      <  1:1352605/150‑6 (MQ=255)
                aaTCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTAttttt                                                                                                                                  >  1:1606187/1‑150 (MQ=255)
                 aTCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTcc                                                                                                                                >  1:3022483/1‑150 (MQ=255)
                  tctcTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTcc                                                                                                                                >  2:832370/1‑150 (MQ=255)
                   ctctTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTccc                                                                                                                               >  1:3282226/1‑150 (MQ=255)
                     ctTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTa                                                                                                                             >  1:2868473/1‑150 (MQ=255)
                     ctTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTa                                                                                                                             <  1:365460/150‑1 (MQ=255)
                     ctTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTa                                                                                                                             <  1:3040037/150‑1 (MQ=255)
                     ctTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTa                                                                                                                             <  1:3026699/150‑1 (MQ=255)
                        gtaagtaaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACa                                                                                                                          >  2:1602010/1‑150 (MQ=255)
                           agtaaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAat                                                                                                                       >  1:3863230/1‑150 (MQ=255)
                           agtaaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAat                                                                                                                       <  1:1372174/150‑1 (MQ=255)
                            gtaaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAatg                                                                                                                      >  2:790583/1‑150 (MQ=255)
                             taaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTAtttt                                                                                                                                              <  2:3073944/124‑1 (MQ=255)
                             taaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAatgat                                                                                                                    <  1:3120491/150‑1 (MQ=255)
                             taaACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAatgat                                                                                                                    <  1:1131187/150‑1 (MQ=255)
                                 cATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAatgatgatg                                                                                                                >  1:4671053/1‑150 (MQ=255)
                                  aTTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAatgatgatg                                                                                                                >  1:960660/1‑150 (MQ=255)
                                    tACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATggg                                                                                                              <  1:3625796/150‑1 (MQ=255)
                                      cATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATggggg                                                                                                            >  2:105298/1‑150 (MQ=255)
                                          aaCTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAgg                                                                                                        >  1:4345109/1‑150 (MQ=255)
                                            cTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGtt                                                                                                      <  1:3979568/150‑1 (MQ=255)
                                             tataTGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGttt                                                                                                     <  2:1499219/150‑1 (MQ=255)
                                                 tGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGa                                                                                                 >  2:523254/1‑150 (MQ=255)
                                                  ggCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGAGGGGGGTAGGTTTAAGaa                                                                                                >  2:2865895/1‑150 (MQ=255)
                                                     aaTACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAAc                                                                                             <  2:4171095/150‑1 (MQ=255)
                                                     aaTACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAAc                                                                                             <  2:1947481/150‑1 (MQ=255)
                                                      aTACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACa                                                                                            >  1:3186386/1‑150 (MQ=255)
                                                       tACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACAt                                                                                           >  1:2417691/1‑150 (MQ=255)
                                                        aCCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACAtt                                                                                          <  2:1068389/150‑1 (MQ=255)
                                                        aCCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACAtt                                                                                          <  2:3472904/150‑1 (MQ=255)
                                                         ccAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTg                                                                                         <  2:2114475/150‑1 (MQ=255)
                                                            gAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACa                                                                                      >  2:1986882/1‑150 (MQ=255)
                                                                atgatgGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGa                                                                                  >  2:2612803/1‑150 (MQ=255)
                                                                 tgatgGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAg                                                                                 <  2:2058463/150‑1 (MQ=255)
                                                                     ggAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGcc                                                                             >  2:3343902/1‑150 (MQ=255)
                                                                      gAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCa                                                                            <  1:723621/150‑1 (MQ=255)
                                                                      gAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCa                                                                            >  2:47840/1‑150 (MQ=255)
                                                                          ttttACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGa                                                                        >  1:3269801/1‑150 (MQ=255)
                                                                           tttACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGaa                                                                       >  2:662283/1‑150 (MQ=255)
                                                                           tttACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGaa                                                                       <  1:832370/150‑1 (MQ=255)
                                                                             tACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACa                                                                     >  1:1230696/1‑150 (MQ=255)
                                                                             tACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACa                                                                     >  1:572377/1‑150 (MQ=255)
                                                                              aCATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGTACaa                                                                    >  2:4233359/1‑150 (MQ=255)
                                                                                     tCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGa                                                             <  1:3794710/150‑1 (MQ=255)
                                                                                      cccAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGaa                                                            <  1:1145802/150‑1 (MQ=255)
                                                                                      cccAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGaa                                                            <  1:1145803/150‑1 (MQ=255)
                                                                                       ccAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGaaa                                                           >  2:1521645/1‑150 (MQ=255)
                                                                                          gATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATc                                                        >  2:3834145/1‑150 (MQ=255)
                                                                                               ttttCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATc                                                   <  1:3885190/150‑1 (MQ=255)
                                                                                                tttCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCa                                                  >  2:1610782/1‑150 (MQ=255)
                                                                                                tttCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCa                                                  <  2:2746385/150‑1 (MQ=255)
                                                                                                  tCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCaaa                                                <  2:4549536/150‑1 (MQ=255)
                                                                                                      cACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAAt                                            <  1:4290958/150‑1 (MQ=255)
                                                                                                          tAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATc                                        <  1:4637385/150‑1 (MQ=255)
                                                                                                               gTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGAt                                   >  1:722031/1‑150 (MQ=255)
                                                                                                               gTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGAt                                   >  1:722178/1‑150 (MQ=255)
                                                                                                                  gATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGTACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACa                                <  2:2368249/150‑1 (MQ=255)
                                                                                                                  gATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGTACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACa                                <  2:2369176/150‑1 (MQ=255)
                                                                                                                  gATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaa                               >  2:2449161/1‑150 (MQ=255)
                                                                                                                   atatatCAATATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACaa                               >  2:1828632/1‑150 (MQ=255)
                                                                                                                          aaTATACACAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAgt                        >  1:1629017/1‑150 (MQ=255)
                                                                                                                               acacAAATGTAGAATATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTgtggt                   >  1:1407189/1‑150 (MQ=255)
                                                                                                                               acacAAATGTAGAATATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTgtggt                   >  1:1407029/1‑150 (MQ=255)
                                                                                                                               acacAAATGTAGAATATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTgtggt                   >  1:1406955/1‑150 (MQ=255)
                                                                                                                               acacAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTgtggt                   >  2:1250963/1‑150 (MQ=255)
                                                                                                                               acacAAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTgtggt                   >  2:3310810/1‑150 (MQ=255)
                                                                                                                                cacaAATGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGACCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTgtggtg                  >  2:3413176/1‑150 (MQ=255)
                                                                                                                                   aaaTGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGTACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTaa               <  1:4474945/150‑1 (MQ=255)
                                                                                                                                    aaTGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTaaa              >  1:1970783/1‑150 (MQ=255)
                                                                                                                                     aTGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAat             >  1:3172181/1‑150 (MQ=255)
                                                                                                                                      tGTAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAata            >  2:1935130/1‑150 (MQ=255)
                                                                                                                                        tAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATg          >  1:1271192/1‑150 (MQ=255)
                                                                                                                                        tAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATg          >  2:3823717/1‑150 (MQ=255)
                                                                                                                                        tAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATg          >  2:1413915/1‑150 (MQ=255)
                                                                                                                                        tAGAAGATGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATg          <  1:599910/150‑1 (MQ=255)
                                                                                                                                         agaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATgg         >  1:458499/1‑150 (MQ=255)
                                                                                                                                         agaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATgg         >  2:2268883/1‑150 (MQ=255)
                                                                                                                                         agaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATgg         >  2:2268877/1‑150 (MQ=255)
                                                                                                                                          gaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGTACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGt        <  1:4233359/150‑1 (MQ=255)
                                                                                                                                          gaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGt        >  1:1397300/1‑150 (MQ=255)
                                                                                                                                          gaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGt        >  2:1542197/1‑150 (MQ=255)
                                                                                                                                          gaagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGt        <  2:1260566/150‑1 (MQ=255)
                                                                                                                                           aagaTGATTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTa       >  1:2512514/1‑150 (MQ=255)
                                                                                                                                             gatgatTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAg     <  2:4564006/150‑1 (MQ=255)
                                                                                                                                              atgatTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGa    <  1:512892/150‑1 (MQ=255)
                                                                                                                                                gatTATTTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACCAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATa  >  2:2804832/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
TATTTATACAAAGATGAATCTCTTGTAAGTAAACATTACATAAACTATATGGCAATACCAGAAAATGATGGAGTTTTTACATGGCTCCCAGATTTTTTTCCGCACGTAGCGGTGGATATATCAATATACACAAATGTAGAAGATGATTATTTTTTTCTTATTTTTCCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATA  >  CP009273/524186‑524478

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: