Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,371 T→C Y223H (TAC→CAC)  yidC → 60 KD inner‑membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,883,3710TC99.6% 858.0 / NA 267Y223H (TAC→CAC) yidC60 KD inner‑membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (119/147);  minor base A (1/0);  total (120/147)
Fisher's exact test for biased strand distribution p-value = 4.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

TTTGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGTACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTGG  >  NC_000913/3883281‑3883460
                                                                                          |                                                                                         
 ttGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACt                                                                                        >  2:2489062/1‑93 (MQ=255)
 ttGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACt                                                                                        <  2:653835/93‑1 (MQ=255)
   ggTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTcc                                                                                      >  1:5156113/1‑93 (MQ=255)
   ggTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTcc                                                                                      <  1:5544675/93‑1 (MQ=255)
    gTCTGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCa                                                                                     <  2:810937/93‑1 (MQ=255)
    gTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCa                                                                                     <  2:352516/93‑1 (MQ=255)
     tCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCAc                                                                                    <  1:475333/93‑1 (MQ=255)
     tCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCAc                                                                                    <  1:2713855/93‑1 (MQ=255)
       aGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGc                                                                                  >  2:754840/1‑93 (MQ=255)
       aGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGc                                                                                  >  1:1791987/1‑93 (MQ=255)
       aGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGc                                                                                  >  2:3263330/1‑93 (MQ=255)
       aGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGc                                                                                  <  2:3814151/93‑1 (MQ=255)
        gTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACAGCTACCGTGGCGCCGCGCACTCCACGcc                                                                                 >  2:1040197/1‑93 (MQ=255)
        gTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGcc                                                                                 >  2:4833229/1‑93 (MQ=255)
        gTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGcc                                                                                 <  1:3348776/93‑1 (MQ=255)
         ttGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCt                                                                                >  2:2131316/1‑93 (MQ=255)
         ttGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCt                                                                                >  1:4836830/1‑93 (MQ=255)
         ttGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCt                                                                                <  1:3387875/93‑1 (MQ=255)
         ttGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCt                                                                                <  1:2238345/93‑1 (MQ=255)
          tGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTg                                                                               <  1:343230/93‑1 (MQ=255)
           gAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGa                                                                              >  1:6136807/1‑93 (MQ=255)
           gAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGa                                                                              <  1:3533884/93‑1 (MQ=255)
              gCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACga                                                                           <  1:5067754/93‑1 (MQ=255)
              gCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACga                                                                           <  1:2501702/93‑1 (MQ=255)
               cAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgag                                                                          >  1:6005418/1‑93 (MQ=255)
                aaTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgaga                                                                         >  2:3081380/1‑93 (MQ=255)
                aaTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgaga                                                                         >  1:3551753/1‑93 (MQ=255)
                aaTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACgaga                                                                         >  1:237216/1‑93 (MQ=255)
                 aTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGaa                                                                        <  2:4745036/93‑1 (MQ=255)
                 aTCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGaa                                                                        <  1:1466431/93‑1 (MQ=255)
                  tCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAg                                                                       <  1:2849325/93‑1 (MQ=255)
                  tCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAg                                                                       <  2:2102607/93‑1 (MQ=255)
                   ccATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGTCGCGGCGCACTCCACGCCTGACGAGAAGt                                                                      >  1:767782/1‑93 (MQ=255)
                   ccATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGt                                                                      >  2:5900367/1‑93 (MQ=255)
                    cATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTa                                                                     <  1:5463123/93‑1 (MQ=255)
                     aTCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTAt                                                                    >  1:912952/1‑93 (MQ=255)
                     aTCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTAt                                                                    <  1:5897112/93‑1 (MQ=255)
                      tCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATg                                                                   >  2:1395335/1‑93 (MQ=255)
                      tCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATg                                                                   <  1:1629406/93‑1 (MQ=255)
                         ctctGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATgaga                                                                >  2:2315786/1‑93 (MQ=255)
                           ctGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGaaa                                                              <  1:5619058/93‑1 (MQ=255)
                           ctGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGaaa                                                              <  2:5749591/93‑1 (MQ=255)
                           ctGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGaaa                                                              <  1:1991893/93‑1 (MQ=255)
                            tGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAAt                                                             >  2:5665735/1‑93 (MQ=255)
                             gCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATa                                                            <  2:5466001/93‑1 (MQ=255)
                                cccGCATCTCGATACCGGAAGCAGCAGCTTCGCACTGCACCTCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACaa                                                         <  1:5309037/92‑1 (MQ=255)
                                aCCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACaa                                                         >  2:3633259/1‑93 (MQ=255)
                                aCCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACaa                                                         >  2:1835726/1‑93 (MQ=255)
                                aCCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACaa                                                         <  2:4689409/93‑1 (MQ=255)
                                  cGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGt                                                       >  1:4735366/1‑93 (MQ=255)
                                    cATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTc                                                     >  1:6278757/1‑93 (MQ=255)
                                      tctcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGa                                                   >  2:374420/1‑93 (MQ=255)
                                        tcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATa                                                 <  1:2951822/93‑1 (MQ=255)
                                        tcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATa                                                 >  1:323532/1‑93 (MQ=255)
                                        tcGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATa                                                 >  2:115294/1‑93 (MQ=255)
                                         cGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATAc                                                >  2:2025041/1‑93 (MQ=255)
                                           aTACCGGAAGCACCAACTCCGCACTGCACCCCTCCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCa                                              <  1:3564947/93‑1 (MQ=255)
                                            tACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCAt                                             <  2:2541252/93‑1 (MQ=255)
                                            tACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCAt                                             >  2:178758/1‑93 (MQ=255)
                                             aCCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCAtg                                            >  1:5811011/1‑92 (MQ=255)
                                              ccGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTg                                           <  2:5066581/93‑1 (MQ=255)
                                               cGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGc                                          <  1:4001123/93‑1 (MQ=255)
                                               cGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGc                                          >  1:2271962/1‑93 (MQ=255)
                                                ggAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGcc                                         >  2:2825145/1‑93 (MQ=255)
                                                   agcagcAACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGAt                                      <  2:6303175/93‑1 (MQ=255)
                                                    gcagcaACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATa                                     <  2:236947/93‑1 (MQ=255)
                                                     cagcagCTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATaa                                    <  2:942560/93‑1 (MQ=255)
                                                      agcaACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAAc                                   >  1:1207860/1‑93 (MQ=255)
                                                        caACTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGa                                 <  2:1981345/93‑1 (MQ=255)
                                                         aaCTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaa                                <  1:1521385/93‑1 (MQ=255)
                                                          aCTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaaa                               >  1:4182422/1‑93 (MQ=255)
                                                           cTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaaaa                              <  2:2651811/93‑1 (MQ=255)
                                                           cTTCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGaaaa                              >  1:2986402/1‑93 (MQ=255)
                                                             tCGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAAcc                            <  2:5374607/93‑1 (MQ=255)
                                                              cGCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCt                           <  1:5009018/93‑1 (MQ=255)
                                                               gCACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTg                          <  1:84906/93‑1 (MQ=255)
                                                                cACTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATGCCATTGCCGATAACGAAAACCTGa                         <  1:3179186/93‑1 (MQ=255)
                                                                 aCTGCACCCCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGaa                        <  2:4454682/93‑1 (MQ=255)
                                                                  cTGCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAAc                       <  1:183799/93‑1 (MQ=255)
                                                                    gCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAt                     >  2:3189674/1‑93 (MQ=255)
                                                                    gCACACCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAt                     <  1:4613370/93‑1 (MQ=255)
                                                                      acacCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAtct                   <  2:4154122/93‑1 (MQ=255)
                                                                       cacCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACAtctc                  >  1:4559236/1‑93 (MQ=255)
                                                                        cccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATctct                 <  1:3970781/92‑1 (MQ=255)
                                                                        acCTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATctct                 >  1:3829961/1‑93 (MQ=255)
                                                                         ccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCtt                >  2:3915664/1‑93 (MQ=255)
                                                                         ccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCtt                <  1:1718453/93‑1 (MQ=255)
                                                                         ccTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCtt                <  1:5265158/93‑1 (MQ=255)
                                                                          cTTCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTc               >  1:4356627/1‑93 (MQ=255)
                                                                           ttCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCg              >  1:1257350/1‑93 (MQ=255)
                                                                           ttCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCg              >  2:172042/1‑93 (MQ=255)
                                                                            tCCGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGa             >  2:3770728/1‑93 (MQ=255)
                                                                              cGTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGaaa           >  2:149184/1‑93 (MQ=255)
                                                                               gTGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAg          <  1:1143398/93‑1 (MQ=255)
                                                                                tGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAgg         >  2:4609779/1‑93 (MQ=255)
                                                                                tGGCGCGGCGCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAgg         <  1:1940500/93‑1 (MQ=255)
                                                                                  gcgcGGCGCACTCCACTCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAggtg       <  2:1314777/93‑1 (MQ=255)
                                                                                     cggcgCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGtt    <  1:5692864/93‑1 (MQ=255)
                                                                                      ggcgCACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTg   >  1:2676063/1‑93 (MQ=255)
                                                                                       gcgcACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTgg  <  2:2904352/93‑1 (MQ=255)
                                                                                          |                                                                                         
TTTGGTCAGTTGAAGCAATCCATCACTCTGCCACCGCATCTCGATACCGGAAGCAGCAACTTCGCACTGCACACCTTCCGTGGCGCGGCGTACTCCACGCCTGACGAGAAGTATGAGAAATACAAGTTCGATACCATTGCCGATAACGAAAACCTGAACATCTCTTCGAAAGGTGGTTGG  >  NC_000913/3883281‑3883460

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