Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,596,510 A→G 6.5% pseudogene (1577/5422 nt) yneO ← AIDA‑I family autotransporter YneO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,596,5100AG 6.5% 326.3 / 9.6 107pseudogene (1577/5422 nt)yneOAIDA‑I family autotransporter YneO
Reads supporting (aligned to +/- strand):  ref base A (57/43);  new base G (2/5);  total (59/48)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 108 total aligned reads displayed.

TTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCC  >  NC_000913/1596362‑1596659
                                                                                                                                                    |                                                                                                                                                     
tttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGAt                                                                                                                                                      >  2:75971/1‑150 (MQ=255)
 ttAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTa                                                                                                                                                     <  1:144780/150‑1 (MQ=255)
   aagtataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTACTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATaaa                                                                                                                                                   <  2:1269969/147‑1 (MQ=255)
   aaCTATACCTCCGGTGGTAATATTCAGTTCGCCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATaaa                                                                                                                                                   >  2:1045201/1‑150 (MQ=255)
    aCTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTATCGGAGGAAATGCCTCCGGTaagg                                                                                                                                                  <  1:550051/150‑3 (MQ=255)
    aCTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTaagg                                                                                                                                                  <  1:1408714/150‑3 (MQ=255)
     cTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTaaggg                                                                                                                                                 <  1:818584/150‑4 (MQ=255)
     cTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGAAATGCCTCCGGTaaggg                                                                                                                                                 <  1:824106/150‑4 (MQ=255)
      tataCCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATcc                                                                                                                                                <  1:545898/150‑1 (MQ=255)
                 tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                                     <  1:910730/150‑1 (MQ=255)
                 tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                                     >  1:945581/1‑150 (MQ=255)
                 tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                                     >  1:552220/1‑150 (MQ=255)
                 tggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTg                                                                                                                                     <  1:1162420/150‑1 (MQ=255)
                  ggtAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGa                                                                                                                                    >  1:702585/1‑150 (MQ=255)
                      atatTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTa                                                                                                                                <  1:1627874/150‑1 (MQ=255)
                        atTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                                                                                              >  2:834017/1‑150 (MQ=255)
                        atTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                                                                                              <  2:74857/150‑1 (MQ=255)
                        atTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                                                                                              >  1:74857/1‑150 (MQ=255)
                        atTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAg                                                                                                                              >  1:1373962/1‑150 (MQ=255)
                                cACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCAc                                                                                                                      >  2:633655/1‑150 (MQ=255)
                                cACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCAc                                                                                                                      >  1:987605/1‑150 (MQ=255)
                                  cccGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCa                                                                                                                    >  1:1095119/1‑150 (MQ=255)
                                  cccGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCa                                                                                                                    >  2:1005683/1‑150 (MQ=255)
                                  cccGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCa                                                                                                                    >  1:1275289/1‑150 (MQ=255)
                                    cGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCAtt                                                                                                                  <  2:159378/150‑1 (MQ=255)
                                      ttccTAAAACACAGAACTGTAGCAATTGTGCGTTAGTACATGATTTATTTAAATTAAAGATACTCTCCGTTGAAATATAAACTATCCCCAAACCGGAGGCATATACTCCGAAAAAACCATACACCATTGAATTAAGGTTGTCACCATTaa                                                                                                                >  2:175773/3‑149 (MQ=255)
                                        ccTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTg                                                                                                              <  1:405176/150‑1 (MQ=255)
                                        ccTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTg                                                                                                              <  2:443704/150‑1 (MQ=255)
                                         cTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGa                                                                                                             >  1:1177916/1‑150 (MQ=255)
                                           aaTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACa                                                                                                           >  1:204340/1‑150 (MQ=255)
                                           aaTACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACa                                                                                                           >  1:299616/1‑150 (MQ=255)
                                               cacCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGGGGCACCATTACTGACATTTa                                                                                                       >  1:510113/1‑150 (MQ=255)
                                                acCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTAc                                                                                                      <  2:204340/150‑1 (MQ=255)
                                                     ccTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATAc                                                                                                 >  1:478870/1‑150 (MQ=255)
                                                       tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATAcag                                                                                               >  1:1331566/1‑148 (MQ=255)
                                                       tGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCg                                                                                               >  2:1503608/1‑150 (MQ=255)
                                                        gTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTAAGATACCGc                                                                                              >  2:784530/1‑150 (MQ=255)
                                                          aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTTAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTa                                                                                            >  2:1026753/1‑150 (MQ=255)
                                                          aGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTa                                                                                            >  2:257487/1‑150 (MQ=255)
                                                           gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTaa                                                                                           >  2:205674/1‑149 (MQ=255)
                                                           gCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCACTTACCGGAGGCATTTCCTCCGATAAATCCATAGAACGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTAc                                                                                           >  1:1120573/1‑150 (MQ=255)
                                                                 gtgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAt                                                                                     <  1:1064800/150‑1 (MQ=255)
                                                                  tgCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAAtt                                                                                    >  1:1380115/1‑150 (MQ=255)
                                                                    cGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAc                                                                                  >  1:486592/1‑150 (MQ=255)
                                                                     gTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTAcc                                                                                 >  1:1289561/1‑150 (MQ=255)
                                                                        agtagATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGAc                                                                              <  2:187919/150‑1 (MQ=255)
                                                                           agATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACAtt                                                                           >  1:1214966/1‑150 (MQ=255)
                                                                             atgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGa                                                                         <  1:633655/150‑1 (MQ=255)
                                                                              tgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                                                        <  2:742728/150‑1 (MQ=255)
                                                                              tgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                                                        >  1:1017719/1‑150 (MQ=255)
                                                                              tgatgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGAt                                                                        >  1:642705/1‑150 (MQ=255)
                                                                                 tgTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAcc                                                                     <  1:1005683/150‑1 (MQ=255)
                                                                                  gTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATAccc                                                                    <  2:1403370/150‑1 (MQ=255)
                                                                                   tCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCca                                                                   >  2:735274/1‑150 (MQ=255)
                                                                                    tcttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                                  >  2:561688/3‑150 (MQ=255)
                                                                                    tcttAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                                  <  1:561688/148‑1 (MQ=255)
                                                                                    cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                                  <  1:23878/150‑1 (MQ=255)
                                                                                    cTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCcac                                                                  >  1:732110/1‑150 (MQ=255)
                                                                                          aTTCCAGATACTGTCCTTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTTACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATc                                                            <  2:1120573/150‑1 (MQ=255)
                                                                                               agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACg                                                       >  2:1591795/1‑150 (MQ=255)
                                                                                               agagACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATAGAACg                                                       >  1:455865/1‑150 (MQ=255)
                                                                                                     tGTCCGTTGAAATATTAACTATCCCCATACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGTCATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGatgatg                                                 <  2:1078342/150‑1 (MQ=255)
                                                                                                      gTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACAGACATTGATACCCACAGTATCGAACGATGATGt                                                >  1:1382268/1‑150 (MQ=255)
                                                                                                        ccGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCg                                              >  2:315090/1‑150 (MQ=255)
                                                                                                         cGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                                             >  1:517254/1‑150 (MQ=255)
                                                                                                         cGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGc                                             >  2:1336635/1‑150 (MQ=255)
                                                                                                             gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc                                         >  1:955338/1‑150 (MQ=255)
                                                                                                             gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc                                         >  1:1546058/1‑150 (MQ=255)
                                                                                                             gAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatcc                                         <  1:1285757/150‑1 (MQ=255)
                                                                                                                atatTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGcatccatc                                      >  2:1473444/1‑150 (MQ=255)
                                                                                                                  atTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAg                                    <  2:702585/150‑1 (MQ=255)
                                                                                                                        tatggatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAg                              >  2:1238245/8‑150 (MQ=255)
                                                                                                                          tCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAgtg                            <  2:1546058/150‑1 (MQ=255)
                                                                                                                            gatttATCGGAGGAAATGCCTCCGGTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGtgtg                          >  1:605618/4‑150 (MQ=255)
                                                                                                                             ccTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGc                         >  1:114825/1‑150 (MQ=255)
                                                                                                                              cTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc                        <  2:486592/150‑1 (MQ=255)
                                                                                                                              cTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc                        <  1:735274/150‑1 (MQ=255)
                                                                                                                                 atcGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGAt                     >  1:410890/3‑150 (MQ=255)
                                                                                                                                 atcGGAGGAAATGCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGAt                     >  1:1160281/3‑150 (MQ=255)
                                                                                                                                 aCCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGAt                     <  1:1026753/150‑1 (MQ=255)
                                                                                                                                    ggAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc                  <  1:1084631/150‑1 (MQ=255)
                                                                                                                                    ggAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc                  <  1:901118/150‑1 (MQ=255)
                                                                                                                                    ggAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc                  <  2:642705/150‑1 (MQ=255)
                                                                                                                                    ggAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGCCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc                  <  1:501379/150‑1 (MQ=255)
                                                                                                                                       ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc                        <  1:452107/141‑1 (MQ=255)
                                                                                                                                       ggCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGcc                        >  2:452107/1‑141 (MQ=255)
                                                                                                                                        gCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgctt              <  2:1380115/150‑1 (MQ=255)
                                                                                                                                         cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATAcc                  <  2:529773/145‑1 (MQ=255)
                                                                                                                                         cATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttg             >  2:887307/1‑150 (MQ=255)
                                                                                                                                          aTTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgc            <  1:1503608/150‑1 (MQ=255)
                                                                                                                                           tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTATTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgct           <  1:75971/150‑1 (MQ=255)
                                                                                                                                           tttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgct           >  1:1522034/1‑150 (MQ=255)
                                                                                                                                            ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgctt          <  2:1095119/150‑1 (MQ=255)
                                                                                                                                            ttCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCgcttgctt          <  2:1275289/150‑1 (MQ=255)
                                                                                                                                             tcctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCttt         >  2:1621923/1‑150 (MQ=255)
                                                                                                                                               ctccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTAt       <  2:455865/150‑1 (MQ=255)
                                                                                                                                                 ccGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCg     <  2:1289561/150‑1 (MQ=255)
                                                                                                                                                 ccGATAAATCCATAGCCCCTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCg     <  2:1250259/150‑1 (MQ=255)
                                                                                                                                                  cGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGc    >  1:1210720/1‑150 (MQ=255)
                                                                                                                                                    aTAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGccc  <  1:887307/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTTAACTATACCTCCGGTGGTAATATTCAGTTCACCCGTGCCTAATACACCGACCTGTAGCAATTGTGCGTTAGTAGATGATGTCTTTAAATTCCAGAGACTGTCCGTTGAAATATTAACTATCCCCTTACCGGAGGCATTTCCTCCGATAAATCCATAGCCCGTTGAATTAAGGGTGGCACCATTACTGACATTTACGATACCGCTACCAAAATTACCGACATTGATACCCACAGTATCGAACGATGATGTCGCATCCATCAGGTTCAGTGTGCCGATACCGCTTGCTTTATCGCCC  >  NC_000913/1596362‑1596659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: