Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,967,132 A→G H239R (CAC→CGC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,967,1320AG100.0% 237.7 / NA 65H239R (CAC→CGC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (30/35);  total (30/35)

CATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGG  >  NC_000913/3966986‑3967276
                                                                                                                                                  |                                                                                                                                                
cATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGtt                                                                                                                                              <  2:265024/151‑1 (MQ=255)
   tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCa                                                                                                                                            <  1:846026/150‑1 (MQ=255)
   tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCa                                                                                                                                            <  1:57115/150‑1 (MQ=255)
       cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGtt                                                                                                                                        <  2:433326/150‑1 (MQ=255)
       cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGtt                                                                                                                                        <  1:845071/150‑1 (MQ=255)
             aTTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGaaa                                                                                                                                 <  1:987710/151‑1 (MQ=255)
                   tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGt                                                                                                                             <  1:358346/149‑1 (MQ=255)
                    aCAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGAt                                                                                                                          >  1:675115/1‑151 (MQ=255)
                     cAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATc                                                                                                                         >  2:274629/1‑151 (MQ=255)
                       accaccCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCga                                                                                                                       <  1:204110/151‑1 (MQ=255)
                        ccaccCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCgag                                                                                                                      <  2:645494/151‑1 (MQ=255)
                        ccaccCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCg                                                                                                                        >  2:215521/1‑149 (MQ=255)
                              ggATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGcc                                                                                                                >  1:413462/1‑151 (MQ=255)
                                aTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaa                                                                                                              <  1:320081/151‑1 (MQ=255)
                                   gtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc                                                                                                            <  1:985858/150‑1 (MQ=255)
                                     gtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGc                                                                                                          <  1:745706/150‑1 (MQ=255)
                                       gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTg                                                                                                       <  2:139133/151‑1 (MQ=255)
                                         tGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGt                                                                                                     >  1:845144/1‑151 (MQ=255)
                                          gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                                                                                    <  1:198682/151‑1 (MQ=255)
                                             ggTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGa                                                                                                  <  2:712911/150‑1 (MQ=255)
                                                   gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGTTGAAGTTGTTGCTTTTACTTTTGACGAACCCGCATCTCGCCGCGTTCTGGTTGCGGAAATGGTGATCGAGAGGGCCAAACGCCTGGTTGAGcact                                                                                             <  2:194077/149‑2 (MQ=255)
                                                   gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac                                                                                              <  1:592162/148‑1 (MQ=255)
                                                   gctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac                                                                                              <  1:724501/148‑1 (MQ=255)
                                                     tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                          >  2:194985/1‑150 (MQ=255)
                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaa                                                                                        <  1:247354/151‑1 (MQ=255)
                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                         >  1:487202/1‑150 (MQ=255)
                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                         >  2:531545/1‑150 (MQ=255)
                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                         <  2:401811/150‑1 (MQ=255)
                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                          <  1:602901/149‑1 (MQ=255)
                                                         cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaaga                                                                                      <  2:390454/150‑1 (MQ=255)
                                                         cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAc                                                                                     >  2:945888/1‑151 (MQ=255)
                                                         cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGAc                                                                                     <  2:866350/151‑1 (MQ=255)
                                                               aCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc                                                                               <  2:802751/151‑1 (MQ=255)
                                                                 gTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCAt                                                                             >  2:145108/1‑151 (MQ=255)
                                                                   ccGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc                                                                           >  2:432850/1‑151 (MQ=255)
                                                                     ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                                                         <  1:577291/151‑1 (MQ=255)
                                                                        agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc                                                                      <  1:521153/151‑1 (MQ=255)
                                                                          aagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGa                                                                    <  2:393225/151‑1 (MQ=255)
                                                                             tAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTc                                                                 >  1:335279/1‑151 (MQ=255)
                                                                                 cGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATc                                                             >  1:901928/1‑151 (MQ=255)
                                                                                  gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATc                                                             >  1:694087/1‑150 (MQ=255)
                                                                                  gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATc                                                             >  2:900920/1‑150 (MQ=255)
                                                                                     aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc                                                         <  2:35745/151‑1 (MQ=255)
                                                                                     aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc                                                         >  1:424629/1‑151 (MQ=255)
                                                                                      tGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCg                                                        >  2:200090/1‑151 (MQ=255)
                                                                                        cAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTc                                                      >  2:233988/1‑151 (MQ=255)
                                                                                          gCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTg                                                    <  2:675039/151‑1 (MQ=255)
                                                                                            gTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgc                                                  >  2:518693/1‑151 (MQ=255)
                                                                                             tCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcg                                                 >  2:383474/1‑151 (MQ=255)
                                                                                                ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcg                                                 <  1:623609/148‑1 (MQ=255)
                                                                                                  tAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcgc                                            >  1:371506/1‑151 (MQ=255)
                                                                                                  tAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcgc                                            >  1:698982/1‑151 (MQ=255)
                                                                                                   aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCt                                           <  1:592048/151‑1 (MQ=255)
                                                                                                   aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCt                                           <  1:711097/151‑1 (MQ=255)
                                                                                                    aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCtt                                          >  1:248726/1‑151 (MQ=255)
                                                                                                                   gCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTc                           >  1:444147/1‑151 (MQ=255)
                                                                                                                     ttCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCg                         >  2:190612/1‑151 (MQ=255)
                                                                                                                      ccTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCgg                        >  2:184595/2‑151 (MQ=255)
                                                                                                                         aCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCgg                        <  1:753071/148‑1 (MQ=255)
                                                                                                                         aCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGc                       <  2:352409/149‑1 (MQ=255)
                                                                                                                                   aaCCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAgt           >  1:252366/1‑151 (MQ=255)
                                                                                                                                         cATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGCTCCGGCGTCAGGTAAAGTGTTGAc     <  2:819337/151‑1 (MQ=255)
                                                                                                                                           tctcGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCg   >  2:281147/1‑151 (MQ=255)
                                                                                                                                            ctcGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCgg  >  2:658517/1‑151 (MQ=255)
                                                                                                                                             tcGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGATAAAGTGTTGACCg   >  1:439150/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGG  >  NC_000913/3966986‑3967276

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