Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 525,739 G→C 7.3% E827Q (GAG→CAG) ‡ rhsD → Rhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913525,7390GC 7.3% 308.8 / 11.3 95E827Q (GAG→CAG) ‡rhsDRhs protein with DUF4329 family putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (31/57);  new base C (3/4);  total (34/61)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTG  >  NC_000913/525590‑525888
                                                                                                                                                     |                                                                                                                                                     
acGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGGGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTg                                                                                                                                                       <  2:239930/150‑1 (MQ=255)
   cATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTgagc                                                                                                                                                    >  2:1599854/1‑149 (MQ=255)
    aTACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGAc                                                                                                                                                   >  1:1073290/1‑150 (MQ=255)
    aTACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGAc                                                                                                                                                   <  2:540475/150‑1 (MQ=255)
         ccGAGCAGGGGCTGGCAAACCGCGGCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgc                                                                                                                                              <  1:1692238/149‑1 (MQ=255)
         aCGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTgc                                                                                                                                              <  1:1707247/150‑1 (MQ=255)
           gAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcag                                                                                                                                            <  1:401452/150‑3 (MQ=255)
           gAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcag                                                                                                                                            <  2:582595/150‑3 (MQ=255)
            aGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCa                                                                                                                                           <  2:963349/150‑1 (MQ=255)
             gCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTCTCACGGTgcagcc                                                                                                                                          <  2:290211/150‑5 (MQ=255)
              cAGGGGCTCGCAAACCGCGTCACTCCGGACATCCTTCCGCCTGTGGAGTGGCTGACGTATGGCTGCGGTTACCTGGCGGGAATGACGCTGGGGGGGACGCCGCTGGTCGAGTATACGCGTGACAGGCTGCACCGTCTCACGGTgcagcct                                                                                                                                         <  1:1112852/150‑6 (MQ=255)
                ggggCTGGCAAACCGCGCCACGCCGCACAGCCGGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCCGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCtt                                                                                                                                       <  2:1434790/150‑1 (MQ=255)
                 gggCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGCGACGGTGCGCAGCTTc                                                                                                                                      >  2:335476/1‑150 (MQ=255)
                  ggCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                                                                                     <  1:454827/150‑1 (MQ=255)
                  ggCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                                                                                     >  1:203720/1‑150 (MQ=255)
                  ggCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCg                                                                                                                                     <  1:766026/150‑1 (MQ=255)
                    ccggaaaaCCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgc                                                                                                                                   <  2:1372463/145‑1 (MQ=255)
                    cTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgc                                                                                                                                   >  1:632383/1‑150 (MQ=255)
                    cTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgc                                                                                                                                   <  1:1028373/150‑1 (MQ=255)
                    cTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgc                                                                                                                                   <  2:1560041/150‑1 (MQ=255)
                     gggCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGGGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgca                                                                                                                                  <  2:721898/149‑1 (MQ=255)
                      ggCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGgcag                                                                                                                                 >  2:1611685/1‑150 (MQ=255)
                          aaCCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATg                                                                                                                             >  2:534660/1‑150 (MQ=255)
                            ccGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggc                                                                                                                           <  1:1679217/150‑1 (MQ=255)
                             cgcgTCACCCCGGACAGCCGGCCGCCGGCGGAGTGGCTGACGTATGGCAGCGGTAACCTGGCGGGAATGACGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggca                                                                                                                          <  1:698268/150‑1 (MQ=18)
                                gTCACGCCGGCCAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                       <  1:1726889/150‑1 (MQ=255)
                                gTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                       <  1:459893/150‑1 (MQ=255)
                                gTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                       <  2:1508743/150‑1 (MQ=255)
                                gTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATggcaggc                                                                                                                       <  2:1520599/150‑1 (MQ=255)
                                    cGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTa                                                                                                                   >  1:1059461/1‑150 (MQ=255)
                                       cGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGTGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATg                                                                                                                <  1:991983/150‑1 (MQ=255)
                                       cGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATg                                                                                                                <  1:400030/150‑1 (MQ=255)
                                       cGGACAGCCGGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATg                                                                                                                <  2:1162448/150‑1 (MQ=255)
                                         gACAGCCTGCCGCCGGTGGAGTGGCTGCCGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGcc                                                                                                              <  2:1065069/150‑1 (MQ=255)
                                          aCAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCg                                                                                                             <  2:1639778/150‑1 (MQ=255)
                                            cgCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTGACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCa                                                                                                           <  1:488279/149‑1 (MQ=255)
                                             gCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCAt                                                                                                          <  2:517813/150‑1 (MQ=255)
                                               cTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATAc                                                                                                        >  2:979729/1‑150 (MQ=255)
                                                tGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACg                                                                                                       >  1:96534/1‑150 (MQ=255)
                                                   cgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGaa                                                                                                     <  1:925295/149‑1 (MQ=255)
                                                   cgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGaa                                                                                                     >  2:925295/1‑149 (MQ=255)
                                                   cgccgGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAAc                                                                                                    <  2:878538/150‑1 (MQ=255)
                                                        gTGGAGTGGCGGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCGGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGAcc                                                                                               <  1:1536652/150‑1 (MQ=255)
                                                         tGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCa                                                                                              >  1:591095/1‑150 (MQ=255)
                                                         tGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCa                                                                                              <  2:178020/150‑1 (MQ=255)
                                                         tGGAGAGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCACTAATGCCGCATACGAACTGACCa                                                                                              >  1:228796/1‑150 (MQ=255)
                                                         gggAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCa                                                                                              <  1:502736/149‑1 (MQ=255)
                                                           gAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGTCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGc                                                                                            >  2:242409/1‑150 (MQ=255)
                                                           gAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGc                                                                                            >  1:542248/1‑150 (MQ=255)
                                                               ggCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACAt                                                                                        <  1:1206147/150‑1 (MQ=255)
                                                                 cTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATac                                                                                      <  1:1649941/150‑1 (MQ=255)
                                                                  tGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATaca                                                                                     <  2:423667/150‑1 (MQ=255)
                                                                  ggACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATaca                                                                                     <  1:1056818/149‑1 (MQ=255)
                                                                    aCGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACAcc                                                                                   <  2:1084926/150‑1 (MQ=255)
                                                                          ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGCCGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAATTGACCAGCACATACACCCCCGCa                                                                             >  2:644111/1‑150 (MQ=255)
                                                                          ggCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGcc                                                                             >  1:1578828/1‑149 (MQ=255)
                                                                                ggTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAg                                                                       <  1:766170/150‑1 (MQ=255)
                                                                                ggTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAg                                                                       <  1:335476/150‑1 (MQ=255)
                                                                                 gTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGt                                                                      >  1:1240180/1‑150 (MQ=255)
                                                                                 gTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGt                                                                      <  2:1059461/150‑1 (MQ=255)
                                                                                    aCCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTAc                                                                   <  1:129860/150‑1 (MQ=255)
                                                                                     ccTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACa                                                                  <  1:1311831/150‑1 (MQ=255)
                                                                                         gCGGGAAAGAAGCTGGGCGGGACGCCGCTGGACGAGTATAAGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCACTAATGCAGAATACGAACTGCCCCGCAAAAACACCCCCGCAGCCCAGTAACaactc                                                              >  1:235268/1‑146 (MQ=17)
                                                                                          cGGGAATGACGCGGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGcc                                                             <  1:299249/150‑1 (MQ=255)
                                                                                           gggAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCa                                                            <  2:1330181/150‑1 (MQ=255)
                                                                                             gAATGATGCTGGGCGGGTCGCCGCTGGTGGAGTATACGCGGGTCAGGCTGCACCGGGAGAAGGTGCGCAGCTTCGGCAGCATGGCATGCAGTAATGCCGCATACGAACTGACCAGCACATACTCCCCCACAGTCCAGTTACAGAGCCAGc                                                          <  2:1165541/150‑1 (MQ=17)
                                                                                             gAATGAAGCTGGGCGGGACGCCGCTGGTCGAGGATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGCCCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCCGc                                                          >  1:348079/1‑150 (MQ=255)
                                                                                                tGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCAcc                                                       >  2:1324993/1‑150 (MQ=255)
                                                                                                tGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAcc                                                       <  1:766513/150‑1 (MQ=255)
                                                                                                tGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAcc                                                       >  1:522343/1‑150 (MQ=255)
                                                                                                tGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAcc                                                       <  2:691987/150‑1 (MQ=255)
                                                                                                tGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCAcc                                                       <  2:201826/150‑1 (MQ=255)
                                                                                                 gAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCt                                                      >  1:1697146/1‑150 (MQ=255)
                                                                                                        ggCGGGACGCCGCTGGTCGAGTATACGCGGGTCAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGc                                               <  1:1472852/150‑1 (MQ=255)
                                                                                                        ggCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGc                                               <  1:534660/150‑1 (MQ=255)
                                                                                                             gACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGt                                          <  1:629455/150‑1 (MQ=255)
                                                                                                                      gggcGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCg                                 <  1:1684880/147‑1 (MQ=255)
                                                                                                                      ggTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCg                                 >  1:341559/1‑150 (MQ=255)
                                                                                                                      ggTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCg                                 >  2:773233/1‑150 (MQ=255)
                                                                                                                       gTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCCGCCCCTGAACAGCCTGGTATATGACCGt                                >  1:1111546/1‑150 (MQ=255)
                                                                                                                       gTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGt                                >  1:166991/1‑150 (MQ=255)
                                                                                                                         cGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGAACAGCCTGGAATATCACCGTGa                              >  1:151200/1‑150 (MQ=255)
                                                                                                                          gAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGAc                             <  1:103651/150‑1 (MQ=255)
                                                                                                                           cgTATACGCGGGACAGGCTGCACCGTGAGTCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACt                            <  1:2115/149‑1 (MQ=255)
                                                                                                                                 cgcgGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATCGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGt                      <  2:1073290/150‑1 (MQ=255)
                                                                                                                                  gcgGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTg                     >  2:1211695/1‑150 (MQ=255)
                                                                                                                                     ggACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAg                  >  1:1411184/1‑150 (MQ=255)
                                                                                                                                       aCAGCCTGCACCGGGAGACGGCGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGGTACAGAGCCAGCATCTGACCAGCCGGGTATATGACCGTGACTACGGGGGGAGTg                <  2:1021251/150‑1 (MQ=14)
                                                                                                                                          ggCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa             <  2:203720/150‑1 (MQ=40)
                                                                                                                                          ggCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa             >  1:773120/1‑150 (MQ=40)
                                                                                                                                          ggCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa             <  1:883/150‑1 (MQ=40)
                                                                                                                                          ggCTGCACCGGGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACa             <  1:818319/150‑1 (MQ=255)
                                                                                                                                           tgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCCTACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGAACAGCCTGGTATATGACCGTGACGACGGGTGGAGTGACaa            >  1:77675/4‑150 (MQ=14)
                                                                                                                                           tgcgCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGAACAGCCTGGTATATGACCCTGACTACGGGTGGAGTGACaa            >  1:1686627/4‑150 (MQ=14)
                                                                                                                                            cTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACaag           <  1:1046682/150‑2 (MQ=39)
                                                                                                                                              gCACCGTCTCACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACACCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACgg         >  1:912873/1‑150 (MQ=255)
                                                                                                                                                 ccGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGa      >  2:550348/1‑150 (MQ=39)
                                                                                                                                                     gagaCGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGCGCCAGCACCTGAACAGCCTGGTATATGACCGTGCCTACGGGTGGAGTGCCAACGGCGACCTg  >  1:1582921/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGCATACAACGAGCAGGGGCTGGCAAACCGCGTCACGCCGGACAGCCTGCCGCCGGTGGAGTGGCTGACGTATGGCAGCGGTTACCTGGCGGGAATGAAGCTGGGCGGGACGCCGCTGGTCGAGTATACGCGGGACAGGCTGCACCGTGAGACGGTGCGCAGCTTCGGCAGCATGGCAGGCAGTAATGCCGCATACGAACTGACCAGCACATACACCCCCGCAGGCCAGTTACAGAGCCAGCACCTGAACAGCCTGGTATATGACCGTGACTACGGGTGGAGTGACAACGGCGACCTG  >  NC_000913/525590‑525888

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