Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,598 G→A G61E (GGA→GAA)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,5980GA95.6% 1320.8 / 33.5 407G61E (GGA→GAA) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (6/12);  new base A (206/182);  total (213/194)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 416 total aligned reads displayed.

TCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTA  >  NC_000913/3966450‑3966745
                                                                                                                                                    |                                                                                                                                                   
tCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAggatgaat                                                                                                                                                   >  8:1278889/1‑151 (MQ=255)
   gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAAtt                                                                                                                                                  <  6:73943/149‑1 (MQ=255)
   gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAAtt                                                                                                                                                  <  7:441877/149‑1 (MQ=255)
   gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTg                                                                                                                                                >  7:565916/1‑151 (MQ=255)
   gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTg                                                                                                                                                >  6:1164228/1‑151 (MQ=255)
   gAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGAGATCTTTGGTGCTGGCGTACTGGAGATATTGCAGGATGAATTTg                                                                                                                                                >  2:899935/1‑151 (MQ=255)
       tGATCACTCTCGGCGAAAATATGGGCCAGGAAAACCTGGCTCGTATGCGTAACCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAACAGTGGCGAAGATATCATTGGTGATGGCGTACTCGAGATATTGCACGATGAATATGGttt                                                                                                                                            >  5:1444144/1‑151 (MQ=255)
         aTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcc                                                                                                                                          >  3:1180957/1‑151 (MQ=255)
         aTCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcc                                                                                                                                          <  4:144572/151‑1 (MQ=255)
          tCACTCTCGGCGAAAATATGGGGCTGGAAACCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGATTTTGGTTtcct                                                                                                                                         >  6:1583119/1‑151 (MQ=255)
           cACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCCCGCAAAGAGAGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTtcctc                                                                                                                                        >  8:43454/1‑151 (MQ=255)
              tctcGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGt                                                                                                                                     >  1:1012314/1‑151 (MQ=255)
                 cGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGt                                                                                                                                     <  4:445315/148‑1 (MQ=255)
                 cGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTcc                                                                                                                                  <  8:1160888/151‑1 (MQ=255)
                 cGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTACTCCGTTcc                                                                                                                                  >  3:4550/1‑151 (MQ=255)
                   gCGAAAATATGGGGCTGGAAAAGCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGc                                                                                                                                <  4:1618382/151‑1 (MQ=255)
                   gCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGc                                                                                                                                >  8:664199/1‑151 (MQ=255)
                    cGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCa                                                                                                                               <  4:362961/151‑1 (MQ=255)
                     gAAAATATGGGGCTGGACCACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGGACGGGAGCTATGGCAGGATGAGTTTGGTTTCCTCCGTTCGGc                                                                                                                                >  4:89404/1‑149 (MQ=255)
                      aaaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGc                                                                                                                                >  7:1250349/1‑148 (MQ=255)
                       aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGAc                                                                                                                            <  2:1616247/151‑1 (MQ=255)
                       aaaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGAc                                                                                                                            <  3:64349/151‑1 (MQ=255)
                        aaTATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACa                                                                                                                           <  2:1399263/151‑1 (MQ=255)
                           atGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGc                                                                                                                         <  1:549690/150‑1 (MQ=255)
                           atGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCt                                                                                                                        >  5:296370/1‑151 (MQ=255)
                                 cTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAcc                                                                                                                  >  1:912000/1‑151 (MQ=255)
                                   ggAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTc                                                                                                                >  3:1404538/1‑151 (MQ=255)
                                   ggAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTc                                                                                                                <  6:1065495/151‑1 (MQ=255)
                                     aaaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCTGCAGACAGCTCCTACCTCGc                                                                                                              >  6:544228/1‑151 (MQ=255)
                                      aaaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGcc                                                                                                             >  1:1069154/1‑151 (MQ=255)
                                       aaCCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCg                                                                                                            <  4:775090/151‑1 (MQ=255)
                                          cTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTc                                                                                                         >  5:1025955/1‑151 (MQ=255)
                                           tGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTcc                                                                                                        >  6:1465337/1‑151 (MQ=255)
                                                cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                                                                                   <  4:551063/151‑1 (MQ=255)
                                                cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatg                                                                                                   <  5:338601/151‑1 (MQ=255)
                                                cGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgat                                                                                                    <  6:938109/150‑1 (MQ=255)
                                                 gTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGAGGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCtgatga                                                                                                  <  4:121961/151‑1 (MQ=255)
                                                    tGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACAt                                                                                               >  2:653226/1‑151 (MQ=255)
                                                      cGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCACAGACAGCTCCTACCTCGCCGGTCCTGATGACATCt                                                                                             <  1:628362/151‑1 (MQ=255)
                                                         aaGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACg                                                                                          >  5:68634/1‑151 (MQ=255)
                                                            cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAATAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                                                                       <  7:645348/151‑1 (MQ=255)
                                                            cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                                                                       >  5:493564/1‑151 (MQ=255)
                                                            cAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGttt                                                                                       >  1:231552/1‑151 (MQ=255)
                                                             aGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTc                                                                                      <  3:1348104/151‑1 (MQ=255)
                                                               gACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTccc                                                                                    >  5:814775/1‑151 (MQ=255)
                                                                 cATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                                                  <  1:585958/151‑1 (MQ=255)
                                                                  attattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                                                  >  7:112534/1‑150 (MQ=255)
                                                                  attattTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCt                                                                                  >  1:777018/1‑150 (MQ=255)
                                                                     attTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGcc                                                                              <  6:1674317/151‑1 (MQ=255)
                                                                           gCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGt                                                                                         >  8:1519710/1‑134 (MQ=255)
                                                                           gCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATc                                                                         <  2:1180369/150‑1 (MQ=255)
                                                                             cATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCTTTCCGCAGACAGATCCTACCGCGCCGGTCCGGATTAAATATACTTTTCCCCTTGCCAAATccgc                                                                      >  5:997017/1‑151 (MQ=255)
                                                                              aTCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgcc                                                                     >  7:1359185/1‑151 (MQ=255)
                                                                                ccTCAAGCAGCAGGCATAGAGTGGCTAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGAGGAATTTGGTTTCCTCCGTTCCTCAGACAGCTACTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATccgccg                                                                    <  2:879489/150‑1 (MQ=255)
                                                                                  tGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGttt                                                                 >  4:1678079/1‑151 (MQ=255)
                                                                                  tGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGttt                                                                 >  8:518887/1‑151 (MQ=255)
                                                                                   gAAGCAGCACGCTAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCGGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc                                                                <  7:228946/151‑1 (MQ=255)
                                                                                    aaGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTc                                                                >  3:1662127/1‑150 (MQ=255)
                                                                                      gcagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAc                                                             <  3:250331/151‑1 (MQ=255)
                                                                                       cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                                            <  5:1276842/151‑1 (MQ=255)
                                                                                       cagcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAAcc                                                            >  3:778744/1‑151 (MQ=255)
                                                                                        agcaCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCt                                                           <  7:1593031/151‑1 (MQ=255)
                                                                                           aCGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCg                                                        <  6:302202/151‑1 (MQ=255)
                                                                                               aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACt                                                    >  8:853071/1‑151 (MQ=255)
                                                                                               aaaGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACt                                                    >  2:67049/1‑151 (MQ=255)
                                                                                                 agagTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTgg                                                  <  2:237946/151‑1 (MQ=255)
                                                                                                   agTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTg                                                <  8:186419/151‑1 (MQ=255)
                                                                                                     tGGCGAAGATATCTTTTGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCATACAGCTCCTACCTCGCCGGTCCTTATGACATCTACGTTTCCACTAGCCAAATACGCCGTTTCAACCTACTCACTGTTGAt                                              >  2:1017020/1‑151 (MQ=255)
                                                                                                        cGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCCCACTGGTGATAcc                                           >  4:495484/1‑151 (MQ=255)
                                                                                                         gAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATAcc                                           >  7:895093/1‑150 (MQ=255)
                                                                                                             atatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCAtctc                                      >  2:573749/1‑151 (MQ=255)
                                                                                                              tatCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATctct                                     >  2:458369/1‑151 (MQ=255)
                                                                                                                   ttGGTGATTGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTgg                                   <  8:736307/148‑1 (MQ=255)
                                                                                                                   ttGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                 >  1:1236838/1‑150 (MQ=255)
                                                                                                                   ttGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTAGCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGt                                  <  1:1636635/149‑1 (MQ=255)
                                                                                                                   ttGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTg                                    <  3:1059563/147‑1 (MQ=255)
                                                                                                                   ttGGTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                 <  1:132687/150‑1 (MQ=255)
                                                                                                                     ggTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAg                               <  4:1324912/150‑1 (MQ=255)
                                                                                                                     ggTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGCAAg                               <  1:1398040/150‑1 (MQ=255)
                                                                                                                      gTGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTa                                 >  5:115991/1‑147 (MQ=255)
                                                                                                                       tGATGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCACCTAGCCAAATCCGCCGTTTCAAACTCCGCACTGGTGATACCATCTCTGGTAAGAtt                            >  5:152204/1‑151 (MQ=255)
                                                                                                                         aTGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcg                          >  2:1259417/1‑151 (MQ=255)
                                                                                                                          tGGCGTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgc                         >  3:1520125/1‑151 (MQ=255)
                                                                                                                            gCGTACTGGAGAGAGTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCGTGAGGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACGTCCGCACTGGTGATACCATCTCTGGTAAGATTcgccc                       <  6:205160/151‑1 (MQ=255)
                                                                                                                              gttcTGGAGATATTGCATGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGGCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                     <  7:352692/148‑1 (MQ=255)
                                                                                                                              gTACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgc                     >  3:323495/1‑151 (MQ=255)
                                                                                                                               tACTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcgcccgcc                    >  1:379606/1‑151 (MQ=255)
                                                                                                                                aCTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                   <  8:1442067/151‑1 (MQ=255)
                                                                                                                                aCTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                   <  2:1507083/151‑1 (MQ=255)
                                                                                                                                aCTGGAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCccgccg                   >  4:1288302/1‑151 (MQ=255)
                                                                                                                                   ggAGATATTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGaa                 >  2:969788/1‑150 (MQ=255)
                                                                                                                                         atTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt          <  7:981424/151‑1 (MQ=255)
                                                                                                                                         atTGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt          <  4:620473/151‑1 (MQ=255)
                                                                                                                                         atTGCAGGATGAATTTGGTTCCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCGCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGt          >  6:4627/1‑151 (MQ=255)
                                                                                                                                          ttGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg         >  5:580393/1‑151 (MQ=255)
                                                                                                                                           tGCAGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTg         <  3:881034/150‑1 (MQ=255)
                                                                                                                                              aGGATGAATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACg     >  6:671221/1‑151 (MQ=255)
                                                                                                                                                 atgaatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCt   >  2:114289/1‑150 (MQ=255)
                                                                                                                                                 atgaatTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTa  <  4:158453/151‑1 (MQ=255)
                                                                                                                                                 atgaatTTGGTTTCCTCCGTTCAGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGATTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTCGTGATACCAACTCTGGTAAAATTAGCACGCAGAAATAAGTTGAACggt   >  4:1120851/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCTGAAGCAGCACGCAAAGAGTGGCGAAGATATCTTTGGTGATGGCGTACTGGAGATATTGCAGGATGGATTTGGTTTCCTCCGTTCCGCAGACAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTA  >  NC_000913/3966450‑3966745

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