Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 732,503 G→A 5.8% R974Q (CGG→CAG)  rhsC → Rhs protein with 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_000913732,5030GA 5.8% 315.9 / 8.2 103R974Q (CGG→CAG) rhsCRhs protein with putative toxin domain; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base G (49/48);  new base A (3/3);  total (52/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 105 total aligned reads displayed.

ACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCG  >  NC_000913/732383‑732648
                                                                                                                        |                                                                                                                                                 
aCTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAata                                                                                                                        <  3:135359/148‑1 (MQ=32)
  tATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                <  2:251386/138‑1 (MQ=32)
  caTCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                >  1:251386/2‑138 (MQ=17)
       ttACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGt                                                                                                                                             <  3:1096929/120‑1 (MQ=32)
       ttACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGt                                                                                                                                             >  4:1096929/1‑120 (MQ=32)
         aCCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGcc                                                                                                               >  2:250836/1‑148 (MQ=32)
            ggTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCt                                                                                                            >  1:613419/1‑148 (MQ=32)
                  aCCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGATCGa                                                                                                      >  3:428191/1‑148 (MQ=17)
                      tCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                <  2:823827/118‑1 (MQ=33)
                      tCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCgg                                                                                                                                >  1:823827/1‑118 (MQ=32)
                       cACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTgg                                                                                                                                              <  2:61773/103‑1 (MQ=32)
                       cACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTgg                                                                                                                                              >  1:61773/1‑103 (MQ=32)
                         cGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGATCGAAAGTCGc                                                                                               <  4:428191/148‑1 (MQ=25)
                                cTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCttt                                                                                        <  4:860878/148‑1 (MQ=35)
                                cTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCttt                                                                                        >  3:860878/1‑148 (MQ=35)
                                  gACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                                                            <  4:270434/94‑1 (MQ=34)
                                  gACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTg                                                                                                                                            >  3:270434/1‑94 (MQ=34)
                                  gACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTAc                                                                                      >  4:765556/1‑148 (MQ=35)
                                    cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGTACACAATATGCAGAGCCGCTGGTCGAAAGTCGCt                                                                                              <  3:207168/138‑1 (MQ=18)
                                    cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGa                                                                                    <  2:613419/148‑1 (MQ=33)
                                    cAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCAGGTCGAAAGTCGCt                                                                                              >  4:207168/1‑138 (MQ=18)
                                      gagaAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCtt                                                                                         <  4:718788/141‑1 (MQ=32)
                                      gagaAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCtt                                                                                         >  3:718788/1‑141 (MQ=32)
                                        gaAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGATATCTTTACGACCCg                                                                                >  3:139490/1‑148 (MQ=17)
                                         aaaaCCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGc                                                                               >  1:333446/1‑148 (MQ=32)
                                                aCCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCg                                                                        >  3:214347/1‑148 (MQ=32)
                                                    cATCCCGGAAAGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAgg                                                                    >  3:650923/1‑148 (MQ=17)
                                                      tCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                  <  3:765556/148‑1 (MQ=32)
                                                        ccGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                                                                <  1:250836/148‑1 (MQ=32)
                                                         cGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                  <  2:835572/145‑1 (MQ=32)
                                                         cGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                  >  1:835572/1‑145 (MQ=32)
                                                                ggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCt                                                                              <  2:504575/126‑1 (MQ=33)
                                                                ggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCt                                                                              >  1:504575/1‑126 (MQ=33)
                                                                ggTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        <  4:650923/148‑1 (MQ=32)
                                                                     tCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATg                                                   <  4:214347/148‑1 (MQ=34)
                                                                         cACGGATAATGAGCGCACCCACCGATATCATTACGGCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg                                               <  4:139490/148‑1 (MQ=255)
                                                                             gatgatGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCt                                                                                              >  1:125541/1‑97 (MQ=33)
                                                                             gatgatGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCt                                                                                              <  2:125541/97‑1 (MQ=33)
                                                                              atgatgAGCGCACCCACCGGTACCATTACGACATTCAGCACCGGCTGGTGCACTACAAGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg                                               >  2:593614/1‑143 (MQ=25)
                                                                              atgatgAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg                                               <  1:593614/143‑1 (MQ=35)
                                                                               tgatgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTAt                                                    <  4:563423/137‑1 (MQ=33)
                                                                                 atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCt                                                                                              <  4:881894/93‑1 (MQ=34)
                                                                                 atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCt                                                                                              >  3:881894/1‑93 (MQ=34)
                                                                                 atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaaaa                                                           >  3:1089516/1‑128 (MQ=32)
                                                                                 atgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCaaaa                                                           <  4:1089516/128‑1 (MQ=33)
                                                                                  tgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGCCGAAAGTCGCTATCTTTACGAcc                                                                                  <  2:584129/104‑1 (MQ=18)
                                                                                  tgaGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGCCGAAAGTCGCTATCTTTACGAcc                                                                                  >  1:584129/1‑104 (MQ=18)
                                                                                      cgcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    <  3:1130793/146‑1 (MQ=39)
                                                                                      cgcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    >  4:1130793/1‑146 (MQ=38)
                                                                                       gcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGgcggcg                                             >  3:633613/1‑136 (MQ=39)
                                                                                       gcACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGgcggcg                                             <  4:633613/136‑1 (MQ=37)
                                                                                        cacccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCttt                                                                                        <  4:381220/92‑1 (MQ=33)
                                                                                        cacccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCttt                                                                                        >  3:381220/1‑92 (MQ=33)
                                                                                        cacccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        <  4:1018515/124‑1 (MQ=32)
                                                                                        cacccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        >  3:1018515/1‑124 (MQ=32)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                  >  4:865909/1‑112 (MQ=33)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGt                                                                  <  3:865909/112‑1 (MQ=32)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg                                               >  1:149646/1‑131 (MQ=35)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATggcgg                                               <  2:149646/131‑1 (MQ=35)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    >  4:293159/1‑142 (MQ=35)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    <  3:293159/142‑1 (MQ=38)
                                                                                          cccaccGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGAc                              >  2:321608/1‑148 (MQ=35)
                                                                                                      atgaCGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAAGATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    <  1:440740/127‑1 (MQ=25)
                                                                                                      aTTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    >  2:440740/1‑130 (MQ=35)
                                                                                                       ttACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        >  4:758571/1‑109 (MQ=32)
                                                                                                       ttACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        <  3:758571/109‑1 (MQ=33)
                                                                                                        tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                                                                <  1:247164/100‑1 (MQ=33)
                                                                                                        tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                                                                >  2:247164/1‑100 (MQ=33)
                                                                                                        tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACggg                                     >  2:514319/1‑127 (MQ=35)
                                                                                                        tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACggg                                     <  1:514319/127‑1 (MQ=37)
                                                                                                        tACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGc                <  2:333446/148‑1 (MQ=34)
                                                                                                            aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                                                                <  3:507788/96‑1 (MQ=33)
                                                                                                            aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTgg                                                                >  4:507788/1‑96 (MQ=33)
                                                                                                            aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg                                                       >  3:923378/1‑105 (MQ=32)
                                                                                                            aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGa                                          >  4:763472/1‑118 (MQ=37)
                                                                                                            aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGa                                          <  3:763472/118‑1 (MQ=38)
                                                                                                            aCAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTc            >  2:725471/1‑148 (MQ=34)
                                                                                                               gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACggg                                     <  3:40685/120‑1 (MQ=38)
                                                                                                               gTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACggg                                     >  4:40685/1‑120 (MQ=38)
                                                                                                                tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAg                                                                     >  3:560177/1‑87 (MQ=33)
                                                                                                                tCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAg                                                                     <  4:560177/87‑1 (MQ=35)
                                                                                                                 cAGCACCGGCTGGTGCACTACACGCGGACACAATGTGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    <  3:315834/119‑1 (MQ=25)
                                                                                                                 cAGCACCGGCTGGTACACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGAGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    >  4:315834/1‑119 (MQ=255)
                                                                                                                  tgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGc                                                                                               >  3:285197/2‑59 (MQ=11)
                                                                                                                  tgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGc                                                                                               <  4:285197/58‑1 (MQ=11)
                                                                                                                  tgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACg                                                         >  1:172540/2‑97 (MQ=255)
                                                                                                                  tgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACg                                                         <  2:172540/96‑1 (MQ=255)
                                                                                                                  tgCACCAGCCGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTg              >  1:90094/2‑140 (MQ=255)
                                                                                                                    cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        <  2:1016025/96‑1 (MQ=34)
                                                                                                                    cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACgg                                                        >  1:1016025/1‑96 (MQ=34)
                                                                                                                    cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACggg                                     <  3:857080/115‑1 (MQ=38)
                                                                                                                    cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACggg                                     >  4:857080/1‑115 (MQ=38)
                                                                                                                    cACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACgg                                      >  2:968875/1‑114 (MQ=37)
                                                                                                                      ccGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg                                                       <  4:374992/95‑1 (MQ=34)
                                                                                                                      ccGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACggg                                                       >  3:374992/1‑95 (MQ=34)
                                                                                                                      ccGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACg                                       <  2:773779/111‑1 (MQ=37)
                                                                                                                      ccGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACg                                       >  1:773779/1‑111 (MQ=37)
                                                                                                                      ccGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCg  <  1:321608/148‑1 (MQ=35)
                                                                                                                        ggCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    <  3:165506/112‑1 (MQ=37)
                                                                                                                        ggCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGa                                    >  4:165506/1‑112 (MQ=37)
                                                                                                                        |                                                                                                                                                 
ACTATCTTTACCGGTATGACCGTCACGGCAGGCTGACAGAGAAAACCGACCTCATCCCGGAAGGGGTTATCCGCACGGATGATGAGCGCACCCACCGGTACCATTACGACAGTCAGCACCGGCTGGTGCACTACACGCGGACACAATATGCAGAGCCGCTGGTCGAAAGTCGCTATCTTTACGACCCGCTGGGCCGCAGGGTGGCAAAACGGGTATGGCGGCGTGAACGGGACCTGACGGGCTGGATGTCGCTGTCACGGAAACCG  >  NC_000913/732383‑732648

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: