Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,819,035 T→C Y54Y (TAT→TAC)  yicG → UPF0126 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,819,0350TC98.4% 230.0 / NA 63Y54Y (TAT→TAC) yicGUPF0126 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (22/40);  minor base G (0/1);  total (22/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCTT  >  NC_000913/3818891‑3819172
                                                                                                                                                |                                                                                                                                         
gTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTc                                                                                                                                     <  1:558093/151‑1 (MQ=255)
gTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTc                                                                                                                                     <  1:44668/151‑1 (MQ=255)
 tATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCg                                                                                                                                    >  1:683611/1‑151 (MQ=255)
  aTTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCg                                                                                                                                    <  2:33113/150‑1 (MQ=255)
    ttggttggCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCt                                                                                                                                 >  2:965341/1‑151 (MQ=255)
      ggttggCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTg                                                                                                                                <  2:358393/150‑1 (MQ=255)
          ggCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCa                                                                                                                           <  1:1002002/151‑1 (MQ=255)
            cATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCaa                                                                                                                          >  2:60213/1‑150 (MQ=255)
               tACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACGACTCGCTCGGCTGGGTCAAaca                                                                                                                       <  2:919519/150‑1 (MQ=255)
                aCTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACAcc                                                                                                                     >  1:68428/1‑151 (MQ=255)
                      gAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAat                                                                                                               >  2:883966/1‑151 (MQ=255)
                      gAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAat                                                                                                               >  2:277408/1‑151 (MQ=255)
                      gAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAat                                                                                                               >  2:947265/1‑151 (MQ=255)
                       aaGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAat                                                                                                               >  2:929308/1‑150 (MQ=255)
                        aGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAatat                                                                                                             <  1:944818/151‑1 (MQ=255)
                          ccATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATAtgt                                                                                                           <  1:990480/151‑1 (MQ=255)
                            aTGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATAtgtg                                                                                                          >  1:55421/1‑150 (MQ=255)
                            aTGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATGATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCCAACACCCGGAATATGTGa                                                                                                         <  2:727001/151‑1 (MQ=255)
                                cGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTAt                                                                                                     <  1:865259/151‑1 (MQ=255)
                                      cgcTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGc                                                                                               <  1:932957/151‑1 (MQ=255)
                                       gcTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCt                                                                                              >  1:265496/1‑151 (MQ=255)
                                             ggcCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTAccgcc                                                                                        <  2:962434/151‑1 (MQ=255)
                                                 ggACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAg                                                                                    <  1:1061191/151‑1 (MQ=255)
                                                    cgacgGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTAc                                                                                 >  2:186372/1‑151 (MQ=255)
                                                     gacgGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACt                                                                                >  2:539261/1‑151 (MQ=255)
                                                     gacgGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACt                                                                                >  1:816654/1‑151 (MQ=255)
                                                       cgGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTc                                                                               <  2:838379/150‑1 (MQ=255)
                                                         gcgcATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACCCCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTcacc                                                                            <  2:513861/151‑1 (MQ=255)
                                                           gcATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTcaccac                                                                          <  2:225237/151‑1 (MQ=255)
                                                             aTGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTa                                                                        <  1:242734/151‑1 (MQ=255)
                                                                    cATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGc                                                                 >  2:959868/1‑151 (MQ=255)
                                                                    cATTTGGCGGAATTATAATTGGGACCGCCGCCGAAATTGGCGGAGGTTCGGTGCGCGATAGTGTGCGGGGCCGCTACCCGCTCGGCTGGGTCAACCACCCGGAGTATGTGATTATGGTCGCTCCCGCCGCAGTACTCACCACTATCg                                                                     <  2:565363/147‑1 (MQ=255)
                                                                      tttGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCt                                                                <  2:134582/150‑1 (MQ=255)
                                                                       ttGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTcc                                                              >  2:402829/1‑151 (MQ=255)
                                                                       ttGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTc                                                               <  1:482890/150‑1 (MQ=255)
                                                                         ggCGTAATTATTATTGCGCCCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTca                                                              <  2:549374/149‑2 (MQ=255)
                                                                              aaTTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATg                                                       <  1:759206/151‑1 (MQ=255)
                                                                                  attattGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCtt                                                   <  1:248685/151‑1 (MQ=255)
                                                                                  attattGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCtt                                                   >  1:706593/1‑151 (MQ=255)
                                                                                  attattGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCtt                                                   <  1:904009/151‑1 (MQ=255)
                                                                                  attattGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCt                                                    <  1:699044/150‑1 (MQ=255)
                                                                                  ataactGCGACAGTCCCCGCAATTGGGTGAGGGTCATGGCGCGTTATTCTGCAGGGCCACCCGCCGCTCGCCTGGGTCAAACACCCGGCATATGTGATTATCGTGGCCACCGCCGCAGTACTCACCACTATCGATGCACCTGTAATGc                                                      <  2:97128/143‑1 (MQ=255)
                                                                                    tattGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTAc                                                 >  1:887134/1‑151 (MQ=255)
                                                                                        gCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTgc                                             <  2:962655/151‑1 (MQ=255)
                                                                                        gCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTg                                              <  1:293282/150‑1 (MQ=255)
                                                                                           aCCGCCCCCGAAATTGGCGGAGGGTCAGTGCGCGATATCCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTgc                                             <  1:669043/148‑1 (MQ=255)
                                                                                             cGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCaa                                         >  1:900793/1‑150 (MQ=255)
                                                                                              gCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCaa                                         >  1:667618/1‑149 (MQ=255)
                                                                                               ccACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAgt                                      <  2:401177/151‑1 (MQ=255)
                                                                                                    gCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGttt                                  <  2:404395/150‑1 (MQ=255)
                                                                                                      aaTTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTg                               <  1:612126/151‑1 (MQ=255)
                                                                                                        ttGGCGGAGGGTCAGTGGGCGATATTCTGCTGGGCCGCTCCCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTc                                 <  2:775011/147‑1 (MQ=255)
                                                                                                          ggCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGc                           <  2:71938/151‑1 (MQ=255)
                                                                                                             ggAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTc                         <  2:809858/150‑1 (MQ=255)
                                                                                                             ggAGGGTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCg                        >  1:184052/1‑151 (MQ=255)
                                                                                                                gggTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACCCCCGGAATATGTGATTTTCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGAt                      <  2:364678/150‑1 (MQ=255)
                                                                                                                gggTCAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGAt                      <  2:799604/150‑1 (MQ=255)
                                                                                                                    cAGTGCGCGATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATgcgc                  <  2:84077/150‑1 (MQ=255)
                                                                                                                         cgcgATATTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGAc            <  2:788908/151‑1 (MQ=255)
                                                                                                                               atTCTGCTGGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACtggtgg      >  2:281792/1‑151 (MQ=255)
                                                                                                                                  ctgctgGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCt   >  1:876350/1‑151 (MQ=255)
                                                                                                                                    gctgGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTc    <  1:723657/148‑1 (MQ=255)
                                                                                                                                    gctgGGCCACTACCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCtt  <  1:330069/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                         
GTATTTGGTTGGCATTACTGCCGAAGCCATGACGGGGGCGCTGGCGGCCGGACGACGGCGCATGGATACATTTGGCGTAATTATTATTGCGACCGCCACCGCAATTGGCGGAGGGTCAGTGCGCGATATTCTGCTGGGCCACTATCCGCTCGGCTGGGTCAAACACCCGGAATATGTGATTATCGTCGCTACCGCCGCAGTACTCACCACTATCGTTGCTCCTGTAATGCCTTACCTGCGCAAAGTGTTTCTGGTGCTCGATGCGCTCGGACTGGTGGTCTT  >  NC_000913/3818891‑3819172

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑