Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 434,599 T→C V80A (GTA→GCA)  nusB → transcription termination

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913434,5990TC99.7% 1209.0 / ‑5.3 346V80A (GTA→GCA) nusBtranscription termination
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (150/195);  minor base A (1/0);  total (151/195)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

Only 100 of 355 total aligned reads displayed.

GCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAC  >  NC_000913/434510‑434688
                                                                                         |                                                                                         
   gctgGCCGGGGTGGCGACTAATACCGCATCCCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGaaaa                                                                                     <  2:2356850/93‑1 (MQ=255)
   gctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGaaaa                                                                                     >  2:39107/1‑93 (MQ=255)
    ctgGCCGGGGTGGGGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGaaaaa                                                                                    <  2:1542972/93‑1 (MQ=255)
    ctgGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGaaaaa                                                                                    <  2:6515524/93‑1 (MQ=255)
      ggCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGc                                                                                  >  2:5077987/1‑93 (MQ=255)
       tccGGGGTGGGGACTAATACCATAAACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCa                                                                                 <  1:8186868/92‑1 (MQ=255)
       gCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCa                                                                                 >  2:2935600/1‑93 (MQ=255)
       gCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGATGGAAGAACTGGGACAGGCAGAAAAAGCa                                                                                 >  2:3195154/1‑93 (MQ=255)
        ccGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGAGGGGCAGCAAAAGCAg                                                                                >  1:370033/1‑93 (MQ=255)
        ccGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAg                                                                                <  1:1066920/93‑1 (MQ=255)
        ccGGGGTGCCGACTAATACCGCATCCCTCGACGGACTGATGAAGCCATCCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAg                                                                                <  1:4199246/93‑1 (MQ=255)
          ggggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTa                                                                              <  1:3911528/93‑1 (MQ=255)
          ggggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTa                                                                              >  1:7071806/1‑93 (MQ=255)
           gggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTAc                                                                             <  1:1669574/93‑1 (MQ=255)
            ggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACt                                                                            >  2:6972037/1‑93 (MQ=255)
            ggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACt                                                                            <  1:6369210/93‑1 (MQ=255)
            ggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACt                                                                            >  1:2376035/1‑93 (MQ=255)
            ggTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACt                                                                            >  1:834806/1‑93 (MQ=255)
              gggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTgc                                                                          <  1:3061148/92‑1 (MQ=255)
               ggCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTgcg                                                                         >  1:4928117/1‑93 (MQ=255)
                  gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAt                                                                      <  2:6740212/93‑1 (MQ=255)
                  gACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAt                                                                      <  1:7373715/93‑1 (MQ=255)
                   aCTAATCCCCCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     <  2:3433273/93‑1 (MQ=255)
                   aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     <  2:6257535/93‑1 (MQ=255)
                   aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     <  1:5078612/93‑1 (MQ=255)
                   aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     <  2:7063766/93‑1 (MQ=255)
                   aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     >  1:3367488/1‑93 (MQ=255)
                   aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     <  1:5628914/93‑1 (MQ=255)
                   aCTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCAtt                                                                     >  2:3692612/1‑93 (MQ=255)
                     tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTgc                                                                   <  1:2503240/93‑1 (MQ=255)
                     tAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTgc                                                                   >  2:382179/1‑93 (MQ=255)
                      aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTgcg                                                                  >  1:4641208/1‑93 (MQ=255)
                      aaTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTgcg                                                                  >  1:2811414/1‑93 (MQ=255)
                       gtACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTgcgc                                                                 >  2:5317404/2‑93 (MQ=255)
                       aTACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTgcgc                                                                 >  1:4696094/1‑93 (MQ=255)
                        tACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCt                                                                >  2:7866708/1‑93 (MQ=255)
                         aCCGCATCCCCCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTg                                                               <  2:8135254/93‑1 (MQ=255)
                         aCCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTg                                                               <  1:4455187/93‑1 (MQ=255)
                          ccGCATCCCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGt                                                              <  1:1985385/93‑1 (MQ=255)
                            gCATCCCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCACTGTAc                                                            <  1:7755441/93‑1 (MQ=255)
                            ccATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTAc                                                            <  1:6772104/92‑1 (MQ=255)
                             cATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACg                                                           <  1:5244368/93‑1 (MQ=255)
                               tACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGaa                                                         >  2:1009215/1‑93 (MQ=255)
                                    cGACGGACGGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTc                                                    <  2:1409097/93‑1 (MQ=255)
                                       cGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaa                                                 <  1:723917/93‑1 (MQ=255)
                                        ggACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                >  2:8440926/1‑93 (MQ=255)
                                        gaccTGATAAAGCCATCCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTaaa                                                <  1:378048/90‑1 (MQ=255)
                                         gACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                               >  2:6796213/1‑93 (MQ=255)
                                         gACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                               >  1:8160624/1‑93 (MQ=255)
                                         gACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAAc                                               <  2:2482950/93‑1 (MQ=255)
                                           cTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGt                                             <  2:4555832/93‑1 (MQ=255)
                                           cTGATGAAGCCATACCTGTCCCACCTGCTGGAAAAACTGGGACCGGCAGAAAAAGCAGGACTGCGCATTGCGCTCTACGAACTGTCTAAACGt                                             >  1:5932562/1‑93 (MQ=255)
                                            tgatgaAGCCATACTTGTCCCGCCTGCTGGAAGGACTGCGACAGGCAGAAAGAGCAGTACTGCGCATTGCGCTGTACGAACTGTATAAACGTa                                            >  2:7439163/1‑93 (MQ=255)
                                            tgatgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTa                                            >  1:1215638/1‑93 (MQ=255)
                                              atgaAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGc                                          <  2:5197756/93‑1 (MQ=255)
                                                 aaGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAt                                       <  1:6902043/93‑1 (MQ=255)
                                                  aGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAtg                                      <  2:3607823/93‑1 (MQ=255)
                                                     cATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGc                                   >  2:4469865/1‑93 (MQ=255)
                                                      aTACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGcc                                  <  2:5581564/93‑1 (MQ=255)
                                                       tACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCa                                 >  1:5293083/1‑93 (MQ=255)
                                                        aCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                <  1:7603990/93‑1 (MQ=255)
                                                        aCCTGTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCAt                                <  1:2923925/93‑1 (MQ=255)
                                                            gTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaa                            <  2:3298671/93‑1 (MQ=255)
                                                            gTCCCGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACaa                            >  1:2579469/1‑93 (MQ=255)
                                                              cccGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAg                          <  1:8503093/93‑1 (MQ=255)
                                                               ccGCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGt                         <  1:6143766/93‑1 (MQ=255)
                                                                 gCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                       >  2:1722539/1‑93 (MQ=255)
                                                                 gCCTGCTGGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTgg                       <  2:4253349/93‑1 (MQ=255)
                                                                   ctgctgGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACTAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                     >  1:5437849/1‑93 (MQ=255)
                                                                   ctgctgGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGcc                     <  2:553626/93‑1 (MQ=255)
                                                                    tgctgGAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCa                    <  1:2259995/93‑1 (MQ=255)
                                                                        ggAAGAACTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGAGGTGCCATACAAAGTGGCCATTaa                >  2:4824028/1‑93 (MQ=255)
                                                                         gaagaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc               >  1:4580765/1‑93 (MQ=255)
                                                                         gaagaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc               <  1:5996832/93‑1 (MQ=255)
                                                                         gaagaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc               <  2:4141130/93‑1 (MQ=255)
                                                                         gaaaaaCTGGGACAGGCAAAAAAAGCAGTACTGCGCATGGCGCTGTACAAACTGTCAAAACGTAGCGATGTGCCATACAAAGTGGCCATTAAc               <  2:1280366/93‑1 (MQ=255)
                                                                          aagaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg              <  2:3137832/93‑1 (MQ=255)
                                                                          aagaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACg              >  2:1946550/1‑93 (MQ=255)
                                                                           agaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGa             <  1:8717574/93‑1 (MQ=255)
                                                                            gaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa            <  1:901651/93‑1 (MQ=255)
                                                                            gaaCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGaa            >  1:312984/1‑93 (MQ=255)
                                                                             aaCTGGGACAGGCAGAAAAAGAAGTACTGCGCATTGCGCTTTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAg           >  2:7220406/1‑93 (MQ=255)
                                                                              aCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc          <  1:1551536/93‑1 (MQ=255)
                                                                              aCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc          <  1:7862783/93‑1 (MQ=255)
                                                                              aCTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGc          <  1:2686390/93‑1 (MQ=255)
                                                                               cTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACACAGTGGCCATTAACGAAGCg         >  1:8422172/1‑93 (MQ=255)
                                                                               cTGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCg         <  1:1737475/93‑1 (MQ=255)
                                                                                tGGGACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGa        >  2:4314207/1‑93 (MQ=255)
                                                                                 gggACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt       <  1:4043161/93‑1 (MQ=255)
                                                                                 gggACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt       <  2:799012/93‑1 (MQ=255)
                                                                                 gggACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt       >  1:4272998/1‑93 (MQ=255)
                                                                                 gggACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGAt       >  2:5745902/1‑93 (MQ=255)
                                                                                  ggACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc      <  2:6804016/93‑1 (MQ=255)
                                                                                  ggACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATc      >  1:7459038/1‑93 (MQ=255)
                                                                                   gACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAACGTGGCCATTAAAGAAGCGATCg     >  1:3476933/1‑93 (MQ=255)
                                                                                   gACAGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg     <  2:5379676/93‑1 (MQ=255)
                                                                                   gACAGGAAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCg     >  2:1183544/1‑93 (MQ=255)
                                                                                     cAGGCAGAAAAAGCAGTACTGCGAATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGag   <  2:6339166/93‑2 (MQ=255)
                                                                                      aGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGagc  <  2:5963899/93‑3 (MQ=255)
                                                                                      aGGCAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGagc  <  2:2778361/93‑3 (MQ=255)
                                                                                         |                                                                                         
GCTGCTGGCCGGGGTGGCGACTAATACCGCATACCTCGACGGACTGATGAAGCCATACCTGTCCCGCCTGCTGGAAGAACTGGGACAGGTAGAAAAAGCAGTACTGCGCATTGCGCTGTACGAACTGTCTAAACGTAGCGATGTGCCATACAAAGTGGCCATTAACGAAGCGATCGAAC  >  NC_000913/434510‑434688

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