Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,362,693 T→C L751P (CTC→CCC)  yhcD → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,362,6930TC99.7% 1021.5 / ‑5.0 308L751P (CTC→CCC) yhcDputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (166/141);  total (166/142)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

Only 100 of 309 total aligned reads displayed.

CCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCAA  >  NC_000913/3362603‑3362782
                                                                                          |                                                                                         
    cGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATaa                                                                                     >  2:5040496/1‑93 (MQ=255)
      gCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACa                                                                                   >  2:3494468/1‑93 (MQ=255)
      gCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACa                                                                                   <  2:3322853/93‑1 (MQ=255)
      gCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACa                                                                                   >  1:4400237/1‑93 (MQ=255)
       cAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACAt                                                                                  <  1:855096/93‑1 (MQ=255)
       cAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACAt                                                                                  >  1:2076556/1‑93 (MQ=255)
        aaaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATc                                                                                 <  1:1072923/93‑1 (MQ=255)
        aaaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATc                                                                                 >  1:755711/1‑93 (MQ=255)
        aaaaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATc                                                                                 <  1:4246649/93‑1 (MQ=255)
          aaaaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCaa                                                                               <  1:5734276/93‑1 (MQ=255)
            aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                             >  2:3829672/1‑93 (MQ=255)
            aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                             <  1:6245116/93‑1 (MQ=255)
            aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                             <  2:6412126/93‑1 (MQ=255)
            aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                             >  2:1718470/1‑93 (MQ=255)
            aaCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAc                                                                             >  2:2480561/1‑93 (MQ=255)
             aCTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAAcc                                                                            >  1:2491626/1‑93 (MQ=255)
              cTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCg                                                                           <  2:260241/93‑1 (MQ=255)
              cTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCg                                                                           <  2:5418215/93‑1 (MQ=255)
               tGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCgg                                                                          >  1:3611225/1‑93 (MQ=255)
               tGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCgg                                                                          >  1:1266083/1‑93 (MQ=255)
                gCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGa                                                                         >  1:5173357/1‑93 (MQ=255)
                 cccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGAc                                                                        <  1:2922125/93‑1 (MQ=255)
                 cccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGAc                                                                        >  2:1871618/1‑93 (MQ=255)
                  ccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACa                                                                       >  1:5017667/1‑93 (MQ=255)
                  ccAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACa                                                                       >  1:1002770/1‑93 (MQ=255)
                   cAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACAt                                                                      >  2:3102598/1‑93 (MQ=255)
                    agagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATg                                                                     <  1:5816234/93‑1 (MQ=255)
                     gagGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATgg                                                                    >  1:419431/1‑93 (MQ=255)
                      agGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGt                                                                   >  1:2971756/1‑93 (MQ=255)
                      agGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGt                                                                   >  1:1789233/1‑93 (MQ=255)
                         cGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCt                                                                >  2:5959708/1‑93 (MQ=255)
                          gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               <  1:4585286/93‑1 (MQ=255)
                          gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               >  1:5490040/1‑93 (MQ=255)
                          gTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTa                                                               <  2:439875/93‑1 (MQ=255)
                            tGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGa                                                             <  2:14078/93‑1 (MQ=255)
                            tGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGa                                                             <  2:5751386/93‑1 (MQ=255)
                            tGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGa                                                             >  1:4657359/1‑93 (MQ=255)
                             gCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGat                                                            <  1:3101545/93‑1 (MQ=255)
                                aTTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                         >  2:1018648/1‑93 (MQ=255)
                                aTTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGAtata                                                         >  2:1524301/1‑93 (MQ=255)
                                 ttGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATaa                                                        <  2:5169438/93‑1 (MQ=255)
                                    aaGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAg                                                     <  1:2577910/93‑1 (MQ=255)
                                     aGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGa                                                    >  2:2231424/1‑93 (MQ=255)
                                       aTAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATg                                                  >  1:663635/1‑93 (MQ=255)
                                        tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                 >  1:6124534/1‑93 (MQ=255)
                                        tAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGa                                                 <  2:4338589/93‑1 (MQ=255)
                                         aGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGag                                                <  2:3068816/93‑1 (MQ=255)
                                         aGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGag                                                >  1:3677713/1‑93 (MQ=255)
                                          gTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagc                                               <  1:4904495/93‑1 (MQ=255)
                                            gAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcag                                             >  2:1996192/1‑93 (MQ=255)
                                             aaGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGagcagc                                            >  1:5897617/1‑93 (MQ=255)
                                              aGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAgcagca                                           >  2:1140073/1‑93 (MQ=255)
                                              aGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAgcagca                                           >  1:2187835/1‑93 (MQ=255)
                                              aGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAgcagca                                           >  2:5633400/1‑93 (MQ=255)
                                                gCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAAt                                         <  2:121282/93‑1 (MQ=255)
                                                  aaCTATCTGACACCATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAACCCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                       <  2:6602402/93‑1 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGGTAGATATAAAAGATGAGCAGCAATCa                                       >  2:3338557/1‑93 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                       >  1:2371275/1‑93 (MQ=255)
                                                  aaCTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCa                                       >  2:6317762/1‑93 (MQ=255)
                                                     tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                    >  1:1211215/1‑93 (MQ=255)
                                                     tATCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGc                                    <  2:3930337/93‑1 (MQ=255)
                                                      aTCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGcc                                   <  2:3078993/93‑1 (MQ=255)
                                                      aTCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGcc                                   >  1:4954933/1‑93 (MQ=255)
                                                       tCTGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCg                                  <  1:4709205/93‑1 (MQ=255)
                                                         tGACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTa                                <  2:453090/93‑1 (MQ=255)
                                                          gACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTAt                               <  1:4024936/93‑1 (MQ=255)
                                                          gACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTAt                               <  2:2726703/93‑1 (MQ=255)
                                                          gACAACATCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTAt                               >  2:6045396/1‑93 (MQ=255)
                                                            caacaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTc                             >  1:3346952/1‑93 (MQ=255)
                                                            caacaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTc                             <  1:2977770/93‑1 (MQ=255)
                                                              acaTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCAc                           <  1:6406585/93‑1 (MQ=255)
                                                                aTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCt                         <  2:5560530/93‑1 (MQ=255)
                                                                aTCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCt                         <  1:5608135/93‑1 (MQ=255)
                                                                 tCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCtt                        >  2:64591/1‑93 (MQ=255)
                                                                 tCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCtt                        >  1:2307252/1‑93 (MQ=255)
                                                                 tCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCtt                        <  1:5098461/93‑1 (MQ=255)
                                                                 tCAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCtt                        >  2:3691174/1‑93 (MQ=255)
                                                                  cAGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTa                       >  1:700495/1‑93 (MQ=255)
                                                                   aGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTAc                      <  1:3260770/93‑1 (MQ=255)
                                                                   aGTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTAc                      >  1:2662730/1‑93 (MQ=255)
                                                                    gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                     >  2:6647246/1‑93 (MQ=255)
                                                                    gTGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACa                     <  1:6039401/93‑1 (MQ=255)
                                                                     tGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACAt                    <  2:4013424/93‑1 (MQ=255)
                                                                     tGGATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACAt                    >  1:1546405/1‑93 (MQ=255)
                                                                      ggATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACAtt                   <  2:2903885/93‑1 (MQ=255)
                                                                       gATGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTc                  <  2:4193434/93‑1 (MQ=255)
                                                                         tGACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTccc                >  2:4965379/1‑93 (MQ=255)
                                                                          gACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCa               <  1:1378302/93‑1 (MQ=255)
                                                                          gACGGTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCa               >  1:1934466/1‑93 (MQ=255)
                                                                             ggTGTTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaa            >  2:1305510/1‑93 (MQ=255)
                                                                              gtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaag           >  2:5259569/1‑93 (MQ=255)
                                                                              gtgtTGTATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaag           >  2:4899891/1‑93 (MQ=255)
                                                                               tgttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCagaaga          <  2:6175055/93‑1 (MQ=255)
                                                                                gttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGAt         >  1:6514481/1‑93 (MQ=255)
                                                                                gttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGAt         >  2:735623/1‑93 (MQ=255)
                                                                                gttgtATTCCCCAATAACATCAAACCCGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGAt         >  1:6551366/1‑93 (MQ=255)
                                                                                 ttgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATg        <  2:2284420/93‑1 (MQ=255)
                                                                                  tgtATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGc       <  1:5442681/93‑1 (MQ=255)
                                                                                    tATTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCAc     >  2:857720/1‑93 (MQ=255)
                                                                                     aTTCCCCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCAcc    >  1:5278762/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGACGGCAAAAAACTGCCCAGAGGCGTTGCGATTGAAGATAGTGAAGGCAACTATCTGACAACATCAGTGGATGACGGTGTTGTATTCCTCAATAACATCAAACCGGACATGGTGCTAGATATAAAAGATGAGCAGCAATCATGCCGTATTCACCTTACATTCCCAGAAGATGCACCAA  >  NC_000913/3362603‑3362782

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