Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,528,194 A→T 8.1% pseudogene (293/2037 nt) rhsE → pseudogene, Rhs family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,528,1940AT 8.1% 226.3 / 8.2 86pseudogene (293/2037 nt)rhsEpseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base A (47/31);  new base T (5/2);  total (53/33)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

TACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGA  >  NC_000913/1528049‑1528338
                                                                                                                                                 |                                                                                                                                                
ttcaccGACCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCGGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTGCACGCGCGACCGCCTGCACCGTGAGACg                                                                                                                                              <  1:329456/145‑1 (MQ=255)
   aaCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTg                                                                                                                                           >  1:149336/1‑150 (MQ=255)
    aCGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGAGCGACCGCCTGCACCGTGAGACGGTgc                                                                                                                                          >  1:240391/1‑150 (MQ=255)
      gAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTCTCACGGTgcag                                                                                                                                        <  1:540877/150‑3 (MQ=255)
       aaCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCa                                                                                                                                       >  2:1655097/1‑150 (MQ=255)
        aCAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTCTCACGGTgcaggc                                                                                                                                      <  1:1421202/150‑5 (MQ=255)
           ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCtt                                                                                                                                   >  1:651712/1‑150 (MQ=255)
           ggggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCtt                                                                                                                                   >  1:1371325/1‑150 (MQ=255)
             ggCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                                                                                 >  2:1438729/1‑150 (MQ=255)
                  gCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGgcagc                                                                                                                            >  1:1389410/1‑150 (MQ=255)
                    aaaCCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCAt                                                                                                                          <  2:226682/150‑1 (MQ=255)
                      aCCGCGTCACGCCGGACAGTCTGCCGCCAGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATgg                                                                                                                        >  1:1471828/1‑150 (MQ=255)
                       ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGTATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                       >  1:408287/1‑150 (MQ=255)
                       ccGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCACTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                       <  2:829528/150‑1 (MQ=255)
                        cgcgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCCCCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                                                                                      >  1:1679939/1‑150 (MQ=255)
                          cggcACGCCGGACAGGCTGCGGCCGGTGGAGTGGCGGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                                                                                                    <  1:89833/147‑1 (MQ=255)
                          cgTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcagg                                                                                                                    >  1:1481645/1‑150 (MQ=255)
                            tCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggca                                                                                                                  >  2:610405/1‑150 (MQ=255)
                             cACGCCGGACAGTCTGCCGCCGGTGTAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGgcaggcag                                                                                                                 <  1:441473/150‑1 (MQ=255)
                             cACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGgcaggcag                                                                                                                 >  1:346072/1‑150 (MQ=255)
                              aCGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                                                                                <  1:922195/150‑1 (MQ=255)
                              aCGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGt                                                                                                                >  1:1079229/1‑150 (MQ=255)
                               cGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAATa                                                                                                               >  2:16921/1‑150 (MQ=255)
                                    gACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGcc                                                                                                          <  2:1371325/150‑1 (MQ=255)
                                      cAGTCTGCCGACGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATCGCAGGCAGTACTGCCGc                                                                                                        >  1:146170/1‑150 (MQ=255)
                                        gTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCat                                                                                                      >  2:475836/1‑150 (MQ=255)
                                           tGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAtata                                                                                                   <  2:1481645/150‑1 (MQ=255)
                                            gccgccGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaa                                                                                                  <  1:1655097/150‑1 (MQ=255)
                                             ccgccgGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATaaa                                                                                                 <  2:689673/150‑1 (MQ=255)
                                                 cgGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGa                                                                                             <  2:1337245/150‑1 (MQ=255)
                                                   gTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGAcc                                                                                           <  2:583160/150‑1 (MQ=255)
                                                    tGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCa                                                                                          >  2:1350848/1‑150 (MQ=255)
                                                     ggAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAg                                                                                         >  1:1487179/1‑150 (MQ=255)
                                                      gAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                                                                                        <  2:1528260/150‑1 (MQ=255)
                                                      gAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGc                                                                                        <  1:632974/150‑1 (MQ=255)
                                                       aGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGcc                                                                                       >  2:1076403/1‑149 (MQ=255)
                                                        gTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCTGTAATGCCGCATATAAACTGACCAGcac                                                                                      <  2:651712/150‑1 (MQ=255)
                                                          ggCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAt                                                                                    >  2:803890/1‑150 (MQ=255)
                                                            cTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATac                                                                                  <  1:569892/150‑1 (MQ=255)
                                                                 cTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGGGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACAcccc                                                                             <  2:1084016/150‑1 (MQ=255)
                                                                     ggCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGCGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACCCCCCCGcc                                                                         >  1:553535/1‑149 (MQ=255)
                                                                     ggCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCa                                                                         >  1:421992/1‑150 (MQ=255)
                                                                            gCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAAAGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGt                                                                  >  2:259109/1‑150 (MQ=255)
                                                                             cTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGtt                                                                 <  2:149336/150‑1 (MQ=255)
                                                                                 cTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACag                                                             >  1:148534/1‑150 (MQ=255)
                                                                                   tGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACagag                                                           <  2:1487179/150‑1 (MQ=255)
                                                                                     cGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                         <  2:202771/150‑1 (MQ=255)
                                                                                     cGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                         >  1:712741/1‑150 (MQ=255)
                                                                                         cATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGCCCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCCGCCCAAACACCCCCGCAGGCCAGTTACAGAGCCAGCa                                                     >  2:136591/1‑150 (MQ=255)
                                                                                         cATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGcc                                                     >  2:1347463/1‑149 (MQ=255)
                                                                                          ttGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAc                                                    <  2:350938/149‑1 (MQ=255)
                                                                                                         aCGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGgtg                                     <  1:475836/150‑1 (MQ=255)
                                                                                                         aCGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGgtg                                     <  2:1471828/150‑1 (MQ=255)
                                                                                                          cgccgcTGCAGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGAACAGCCTGgtgg                                    >  2:266528/1‑149 (MQ=255)
                                                                                                           gccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGCCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCCGCACCTGAACAGCCTGGTGTa                                   >  1:761924/1‑150 (MQ=255)
                                                                                                           gccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTa                                   <  1:1650099/150‑1 (MQ=255)
                                                                                                           gccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTa                                   >  1:881202/1‑150 (MQ=255)
                                                                                                            ccgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGGGTAt                                  >  1:432529/1‑150 (MQ=255)
                                                                                                             cgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGACCAGCCTGGTGTATg                                 >  1:1533348/1‑150 (MQ=255)
                                                                                                             cgcTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATg                                 >  1:455524/1‑150 (MQ=255)
                                                                                                              gctgctGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCCGGCCAGTTACAGAGCCAGCACCTGAACAGCTTGGTGTATGa                                >  1:809596/1‑150 (MQ=255)
                                                                                                               ctgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAc                               <  2:346072/150‑1 (MQ=255)
                                                                                                               ctgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAc                               <  1:1350848/150‑1 (MQ=255)
                                                                                                                tgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAcc                              >  2:1401914/1‑150 (MQ=255)
                                                                                                                tgctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGAcc                              >  1:1275710/1‑150 (MQ=255)
                                                                                                                  ctgGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCCGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTCTGCCCGt                            >  2:268616/1‑150 (MQ=255)
                                                                                                                    ggAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACCCCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGa                          >  2:464341/1‑150 (MQ=255)
                                                                                                                     gAGTTCACGCGCGACCGCCTGCACCGTTAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACAAACACCCCCGCAGGCCAGTTACAGAGCCAGCACCAGAACATCCTGGTGTATGACCGAGAc                         >  2:1034889/1‑150 (MQ=17)
                                                                                                                         tCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACg                     >  1:1060920/1‑150 (MQ=255)
                                                                                                                           aCGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCTGCAGCATGTCATGCAGTAATGCCTCATATATACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACggg                   <  2:263195/150‑1 (MQ=14)
                                                                                                                            cgcgcgACCGCCTGCTCCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGt                  >  1:981901/1‑150 (MQ=255)
                                                                                                                            cgcgcgACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACCTACACCCCCGCAGGCCAGTTAGAGAGCCCGCACCTGAACAGCCTGGTGTATGACCGTGACTACgggg                  >  1:1413385/1‑149 (MQ=255)
                                                                                                                              cgcgACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTgg                <  2:633948/150‑1 (MQ=255)
                                                                                                                                cgACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGaa              >  2:442326/1‑150 (MQ=255)
                                                                                                                                  aCCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATg            >  2:880963/1‑150 (MQ=255)
                                                                                                                                  aCCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATg            <  2:421992/150‑1 (MQ=255)
                                                                                                                                    cGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGAc          >  2:605760/1‑150 (MQ=255)
                                                                                                                                     gCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:338111/1‑150 (MQ=40)
                                                                                                                                     ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:577686/5‑150 (MQ=255)
                                                                                                                                     ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:422716/5‑150 (MQ=255)
                                                                                                                                     ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:584393/5‑150 (MQ=255)
                                                                                                                                     ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:1247999/5‑150 (MQ=255)
                                                                                                                                     ctgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACa         >  2:1506779/5‑150 (MQ=255)
                                                                                                                                      ccTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACaa        <  1:1401914/150‑1 (MQ=40)
                                                                                                                                           aCCGTGAGACGGTGCGCAGATTCGGCAGCATGGAAGTCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCAAGTTACAGAGCCAGAACCTGAACATCCTGGTGTATGACCGTTACTACGGGTGGAATGACAACGGCg   <  2:1046170/150‑1 (MQ=14)
                                                                                                                                            ccGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGa  >  1:744543/1‑150 (MQ=39)
                                                                                                                                                 |                                                                                                                                                
TACAACGAACAGGGGCTGGCAAACCGCGTCACGCCGGACAGTCTGCCGCCGGTGGAGTGGCTGACCTACGGCAGCGGCTATCTTGCGGGCATGAAGCTGGGCGGGACGCCGCTGCTGGAGTTCACGCGCGACCGCCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATATAAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTGTATGACCGTGACTACGGGTGGAATGACAACGGCGA  >  NC_000913/1528049‑1528338

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