Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,883,686 G→A V328I (GTT→ATT)  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,6860GA99.2% 1414.2 / ‑3.5 400V328I (GTT→ATT) yidC60 KD inner‑membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base A (179/218);  total (179/221)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

Only 100 of 409 total aligned reads displayed.

GGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCC  >  NC_000913/3883596‑3883775
                                                                                          |                                                                                         
 gTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTg                                                                                        >  2:4004694/1‑93 (MQ=255)
 gTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTg                                                                                        >  1:5353661/1‑93 (MQ=255)
  tCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGa                                                                                       <  2:1261384/93‑1 (MQ=255)
   cAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGAt                                                                                      <  1:7499588/93‑1 (MQ=255)
     gACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTa                                                                                    <  2:5179472/93‑1 (MQ=255)
     gACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTa                                                                                    <  2:6965736/93‑1 (MQ=255)
      aCTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTAc                                                                                   <  2:7059389/93‑1 (MQ=255)
       cTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACg                                                                                  <  2:2200436/93‑1 (MQ=255)
        tGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACgt                                                                                 >  2:2680441/1‑92 (MQ=255)
         ggCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCATCTGGATCTGACCATTGATTACggt                                                                                >  2:907478/1‑93 (MQ=255)
         ggCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggt                                                                                >  2:3750344/1‑93 (MQ=255)
         ggCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggt                                                                                <  2:316791/93‑1 (MQ=255)
          gcgcGATGAACAGCACCCTGTGGGTTGGCCCGGAGATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggtt                                                                               <  1:6626562/93‑1 (MQ=255)
           cgcgATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggttg                                                                              <  2:4129640/93‑1 (MQ=255)
             cgATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggttggt                                                                            <  2:6490170/93‑1 (MQ=255)
              gATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggttggtt                                                                           <  2:720242/93‑1 (MQ=255)
              gATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACggttggtt                                                                           >  2:2303198/1‑93 (MQ=255)
                tGAACAGCACCCTGTGGGTTGTCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTtgt                                                                         <  2:5500343/93‑1 (MQ=255)
                tGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGgttggttgg                                                                         >  1:5051041/1‑92 (MQ=255)
                tGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTtgt                                                                         <  1:7587961/93‑1 (MQ=255)
                 gAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTtgtg                                                                        >  1:5808380/1‑93 (MQ=255)
                   aCAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGt                                                                      <  2:2384822/93‑1 (MQ=255)
                     aGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGGTGTGGTTc                                                                    >  2:663972/1‑93 (MQ=255)
                      gCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCa                                                                   <  1:6123302/93‑1 (MQ=255)
                      gCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCa                                                                   <  1:6297191/93‑1 (MQ=255)
                      gCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCa                                                                   <  2:2812640/93‑1 (MQ=255)
                        aCCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCAtc                                                                 >  2:5089091/1‑93 (MQ=255)
                         cccTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCAtct                                                                >  2:7530268/1‑93 (MQ=255)
                           cTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCAtctct                                                              <  2:3894316/93‑1 (MQ=255)
                            tgtgGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCAtctctc                                                             >  1:6988392/1‑93 (MQ=255)
                            tgtgGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCAtctctc                                                             >  2:3500961/1‑93 (MQ=255)
                             gtgGGTTGGCCCGGAAATCCGGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCa                                                            >  1:165603/1‑93 (MQ=255)
                               gggTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGc                                                          <  1:6780315/93‑1 (MQ=255)
                               gggTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGc                                                          <  1:5866266/93‑1 (MQ=255)
                                ggTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGcc                                                         >  2:7878772/1‑93 (MQ=255)
                                   tGGCGCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCt                                                      <  1:7368222/93‑1 (MQ=255)
                                    ggCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTg                                                     >  1:3233760/1‑93 (MQ=255)
                                    ggCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTg                                                     >  2:164781/1‑93 (MQ=255)
                                     gCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGt                                                    <  1:7232271/93‑1 (MQ=255)
                                      cccGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGGCCATTGATTACGGTTGGTTGGGGTTCATCTCTCAGCCGCTGtt                                                   <  2:4854734/93‑1 (MQ=255)
                                       ccGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTc                                                  <  2:6780723/93‑1 (MQ=255)
                                        cGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCa                                                 >  1:1718762/1‑93 (MQ=255)
                                          gAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCaaa                                               <  1:820789/93‑1 (MQ=255)
                                          gAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCaaa                                               >  2:4641428/1‑93 (MQ=255)
                                          gAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCaaa                                               >  1:3764265/1‑93 (MQ=255)
                                             aTCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAActg                                            >  1:2082939/1‑93 (MQ=255)
                                             aTCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAActg                                            >  1:3637986/1‑93 (MQ=255)
                                             aTCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAActg                                            <  2:1018018/93‑1 (MQ=255)
                                              tCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAActgc                                           >  1:4551460/1‑93 (MQ=255)
                                                cAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAActgctg                                         >  2:1897536/1‑93 (MQ=255)
                                                 aGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGa                                        <  1:7658448/93‑1 (MQ=255)
                                                 aGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGa                                        >  1:4725762/1‑93 (MQ=255)
                                                 aGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCACCCGCTGTTCAAAATGCTGa                                        >  2:4782988/1‑93 (MQ=255)
                                                 aGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACAATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGa                                        >  1:1837552/1‑93 (MQ=255)
                                                   gACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGaaa                                      >  2:7347918/1‑93 (MQ=255)
                                                   gACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGaaa                                      >  2:298617/1‑93 (MQ=255)
                                                   gACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGaaa                                      <  2:449166/93‑1 (MQ=255)
                                                     cAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATg                                    <  1:5440589/93‑1 (MQ=255)
                                                     cAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATg                                    <  2:2080259/93‑1 (MQ=255)
                                                     cAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATg                                    <  1:6376934/93‑1 (MQ=255)
                                                      aaaaTGGCAGCTGTTTCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATgg                                   <  1:1953931/93‑1 (MQ=255)
                                                      aaaaTGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATgg                                   >  1:5379522/1‑93 (MQ=255)
                                                        aaTGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGAt                                 >  1:3414006/1‑93 (MQ=255)
                                                         aTGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATc                                >  1:3007228/1‑93 (MQ=255)
                                                          tGGCAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATcc                               >  1:1294687/1‑93 (MQ=255)
                                                             cAGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGCGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATa                            <  1:806514/93‑1 (MQ=255)
                                                              aGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAg                           <  2:6282909/93‑1 (MQ=255)
                                                              aGCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAg                           >  2:2570771/1‑93 (MQ=255)
                                                               gCTGTTGCTCCGCACCTGGATCTGGCCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGc                          >  1:1377918/1‑93 (MQ=255)
                                                               gCTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGc                          >  1:2800867/1‑93 (MQ=255)
                                                                cTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCt                         >  1:7343668/1‑93 (MQ=255)
                                                                cTGTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCt                         >  2:3648513/1‑93 (MQ=255)
                                                                  gTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCttt                       <  2:1144163/93‑1 (MQ=255)
                                                                  gTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCttt                       >  1:2519636/1‑93 (MQ=255)
                                                                  gTTGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCttt                       >  1:5943343/1‑93 (MQ=255)
                                                                   ttGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTtg                      >  1:7761875/1‑93 (MQ=255)
                                                                   ttGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTtg                      <  2:1529856/93‑1 (MQ=255)
                                                                    tGCTCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTtgt                     <  2:5705138/93‑1 (MQ=255)
                                                                     gCTCCGAACCTCGCTCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTtgtg                    <  2:2462296/93‑1 (MQ=255)
                                                                       tCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTggg                  <  1:1510638/93‑1 (MQ=255)
                                                                       tCCGCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTggg                  <  2:1957752/93‑1 (MQ=255)
                                                                          gCACCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTaa               >  2:337939/1‑93 (MQ=255)
                                                                            aCCTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACt             >  2:5878391/1‑93 (MQ=255)
                                                                             ccTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTg            <  2:4745595/93‑1 (MQ=255)
                                                                             ccTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTg            >  1:265317/1‑93 (MQ=255)
                                                                             ccTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTg            <  1:3984707/93‑1 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATCGCTTTGTGGGTAACTgg           >  1:4864142/1‑93 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg           >  1:1193745/1‑93 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg           <  1:2140169/93‑1 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg           <  1:4264923/93‑1 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg           <  2:5353660/93‑1 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg           <  1:5617291/93‑1 (MQ=255)
                                                                              cTGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTgg           <  1:6483371/93‑1 (MQ=255)
                                                                               tGGATCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTggg          <  1:2295592/93‑1 (MQ=255)
                                                                                  aTCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCt       <  2:3048656/93‑1 (MQ=255)
                                                                                   tCTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCtt      <  1:1410716/93‑1 (MQ=255)
                                                                                   tCTGACCATTGATTACGGGTGGTTGTGGGTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCtt      >  1:4403328/1‑93 (MQ=255)
                                                                                    cTGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTtc     <  1:7729327/93‑1 (MQ=255)
                                                                                     tGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGTTAACTGGGGCTtct    >  1:1079091/1‑93 (MQ=255)
                                                                                     tGACCATTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTtct    >  2:6910719/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGTCAGACTGGCGCGATGAACAGCACCCTGTGGGTTGGCCCGGAAATCCAGGACAAAATGGCAGCTGTTGCTCCGCACCTGGATCTGACCGTTGATTACGGTTGGTTGTGGTTCATCTCTCAGCCGCTGTTCAAACTGCTGAAATGGATCCATAGCTTTGTGGGTAACTGGGGCTTCTCC  >  NC_000913/3883596‑3883775

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