Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 581,068 C→T P79S (CCA→TCA)  nohD → DLP12 prophage; putative DNA‑packaging protein NohD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913581,0680CT100.0% 796.9 / NA 241P79S (CCA→TCA) nohDDLP12 prophage; putative DNA‑packaging protein NohD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (126/115);  total (126/115)

Only 100 of 255 total aligned reads displayed.

TCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAAC  >  NC_000913/580970‑581165
                                                                                                  |                                                                                                 
tCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTccagtca                                                                                                 >  1:771046/1‑101 (MQ=255)
 cATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTccagtcag                                                                                                >  2:771048/1‑101 (MQ=255)
 cATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTccagtca                                                                                                 <  2:770755/100‑1 (MQ=255)
  aTAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAgg                                                                                               >  1:771050/1‑101 (MQ=255)
   tAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAggg                                                                                              <  1:770724/101‑1 (MQ=255)
      aaTGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACt                                                                                           >  1:771054/1‑101 (MQ=38)
      aaTGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACt                                                                                           <  2:770730/101‑1 (MQ=255)
       aTGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTa                                                                                          >  1:771057/1‑101 (MQ=25)
       aTGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTa                                                                                          >  2:771056/1‑101 (MQ=25)
        tGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTAt                                                                                         >  1:771062/1‑101 (MQ=25)
         ggTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTAtt                                                                                        >  1:771065/1‑101 (MQ=25)
          gTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTg                                                                                       >  2:771066/1‑101 (MQ=25)
           tATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGa                                                                                      <  2:770746/101‑1 (MQ=255)
           tATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGa                                                                                      >  1:771069/1‑101 (MQ=25)
                 gAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGa                                                                                >  1:771073/1‑101 (MQ=38)
                   aaGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAAc                                                                              >  1:771076/1‑101 (MQ=25)
                   aaGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAAc                                                                              >  2:771074/1‑101 (MQ=25)
                     gggATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGc                                                                            <  1:770573/101‑1 (MQ=255)
                      ggATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGcc                                                                           <  1:770526/101‑1 (MQ=255)
                      ggATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGcc                                                                           >  1:771079/1‑101 (MQ=38)
                       gATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCa                                                                          <  2:770871/101‑1 (MQ=255)
                       gATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCa                                                                          <  2:770901/101‑1 (MQ=25)
                         tGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATc                                                                        <  2:771004/101‑1 (MQ=255)
                          gCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCg                                                                       >  1:771082/1‑101 (MQ=25)
                           cTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGa                                                                      >  1:771085/1‑101 (MQ=38)
                            tGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGAc                                                                     >  2:771087/1‑101 (MQ=255)
                              aaaTTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACtt                                                                   >  1:771088/1‑101 (MQ=38)
                              aaaTTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACtt                                                                   <  1:770864/101‑1 (MQ=255)
                                aTTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTAc                                                                 >  2:771090/1‑101 (MQ=38)
                                aTTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTAc                                                                 <  1:770778/101‑1 (MQ=25)
                                  tGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGa                                                                      >  1:771095/1‑94 (MQ=255)
                                  tGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTAcgc                                                               <  1:770649/101‑1 (MQ=25)
                                  tGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTAcgc                                                               <  1:770389/101‑1 (MQ=25)
                                     gaACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTg                                                            >  2:771096/1‑101 (MQ=25)
                                      aaCGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTgc                                                           <  1:770586/101‑1 (MQ=25)
                                      aaCGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTgc                                                           >  1:771098/1‑101 (MQ=255)
                                      aaCGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTgc                                                           >  2:771099/1‑101 (MQ=255)
                                        cGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTgcga                                                         <  2:770777/101‑2 (MQ=25)
                                         gAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCa                                                        <  2:770662/101‑1 (MQ=25)
                                         gAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCa                                                        >  2:771103/1‑101 (MQ=25)
                                          aaaaGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAg                                                       <  1:770692/101‑1 (MQ=25)
                                            aaGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGc                                                     >  2:771106/1‑101 (MQ=25)
                                             aGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGcc                                                    >  1:771107/1‑101 (MQ=25)
                                             aGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGcc                                                    <  1:770825/101‑1 (MQ=25)
                                               cTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACg                                                                             >  2:771111/1‑74 (MQ=255)
                                                tGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTAc                                                                                  >  1:771113/1‑68 (MQ=255)
                                                   gcCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGca                                              <  2:770801/101‑1 (MQ=255)
                                                    ccGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGcac                                             >  1:771117/1‑101 (MQ=255)
                                                     cGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGcaca                                            >  1:771121/1‑101 (MQ=255)
                                                     cGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGcaca                                            >  2:771119/1‑101 (MQ=255)
                                                     cGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGcaca                                            <  2:770556/101‑1 (MQ=255)
                                                      gggAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAg                                           <  1:770636/101‑1 (MQ=255)
                                                       ggAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAgg                                          <  1:770785/101‑1 (MQ=255)
                                                       ggAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAgg                                          >  1:771124/1‑101 (MQ=255)
                                                        gAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGa                                         >  1:771127/1‑101 (MQ=255)
                                                         aGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAg                                        >  2:771128/1‑101 (MQ=255)
                                                           gTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCt                                      >  2:771131/1‑101 (MQ=38)
                                                           gTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCt                                      <  2:770648/101‑1 (MQ=25)
                                                            ttGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTg                                     >  1:771134/1‑101 (MQ=25)
                                                             tGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGa                                    >  2:771136/1‑101 (MQ=25)
                                                              gaagaaCTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGaa                                   <  2:770393/101‑1 (MQ=38)
                                                               aagaaCTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGaaa                                  >  1:771138/1‑101 (MQ=25)
                                                                 gaaCTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGaaaaa                                >  1:771141/1‑101 (MQ=25)
                                                                   aCTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATg                              <  1:770819/101‑1 (MQ=38)
                                                                   aCTGCTGCAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATg                              <  1:770681/101‑1 (MQ=25)
                                                                     tgctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGcc                            <  2:770766/101‑1 (MQ=38)
                                                                      gctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGa                                    <  2:770669/92‑1 (MQ=37)
                                                                      gctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCa                           <  1:770672/101‑1 (MQ=25)
                                                                      gctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCa                           <  1:770754/101‑1 (MQ=25)
                                                                       ctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAg                                           <  1:770767/84‑1 (MQ=37)
                                                                       ctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCag                          <  2:770710/101‑1 (MQ=25)
                                                                       ctgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCag                          <  1:770707/101‑1 (MQ=25)
                                                                        tgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCaga                         <  2:770942/101‑1 (MQ=25)
                                                                        tgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCaga                         >  2:771146/1‑101 (MQ=18)
                                                                        tgcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCaga                         >  1:771144/1‑101 (MQ=18)
                                                                         gcAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCagag                        >  2:771148/1‑101 (MQ=17)
                                                                          cAGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAgaga                       >  1:771150/1‑101 (MQ=17)
                                                                           aGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGAc                      >  1:771154/1‑101 (MQ=17)
                                                                           aGGCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGAc                      <  1:770853/101‑1 (MQ=18)
                                                                             gCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTc                    <  2:770791/101‑1 (MQ=18)
                                                                             gCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTc                    >  2:771159/1‑101 (MQ=17)
                                                                             gCCAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTc                    >  1:771158/1‑101 (MQ=17)
                                                                              ccAGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTcc                   <  2:770721/101‑1 (MQ=18)
                                                                                aGCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGc                 >  1:771161/1‑101 (MQ=17)
                                                                                 gCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTcc                   <  1:770798/98‑1 (MQ=18)
                                                                                 gCGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCt                <  1:770807/101‑1 (MQ=18)
                                                                                  cGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTg               <  2:770841/101‑1 (MQ=18)
                                                                                  cGAGACAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTg               <  1:770833/101‑1 (MQ=17)
                                                                                   gagaCAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGa              >  2:771166/1‑101 (MQ=17)
                                                                                   gagaCAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGa              >  1:771165/1‑101 (MQ=17)
                                                                                      aCAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAgt           <  2:770604/101‑1 (MQ=18)
                                                                                       cAGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAgtg          <  1:770889/101‑1 (MQ=18)
                                                                                        aGATCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAgtgg         >  1:771170/1‑101 (MQ=17)
                                                                                          aTCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAgtggtg       >  1:771173/1‑101 (MQ=17)
                                                                                          aTCTCCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAgtggtg       <  2:770690/101‑1 (MQ=18)
                                                                                           tctcCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGtggtgg      <  2:770743/101‑1 (MQ=18)
                                                                                            ctcCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGa     >  2:771174/1‑101 (MQ=17)
                                                                                             tcCAGTCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGaa    <  1:770914/101‑1 (MQ=18)
                                                                                              ccagtcagGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGaaa   <  2:770970/101‑1 (MQ=18)
                                                                                               cagtcagGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAAc  <  2:770893/101‑1 (MQ=18)
                                                                                                  |                                                                                                 
TCATAAAATGGTATGCCGAAAGGGATGCTGAAATTGAGAACGAAAAGCTGCGCCGGGAGGTTGAAGAACTGCTGCAGGCCAGCGAGACAGATCTCCAGCCAGGGACTATTGAGTACGAACGCCATCGACTTACGCGTGCGCAGGCCGATGCACAGGAGCTGAAAAATGCCAGAGACTCCGCTGAAGTGGTGGAAAC  >  NC_000913/580970‑581165

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