Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,434,318 C→T W200* (TGG→TGA b1377 ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,434,3180CT88.4% 775.2 / 59.1 258W200* (TGG→TGAb1377putative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base C (21/9);  new base T (118/110);  total (139/119)
Fisher's exact test for biased strand distribution p-value = 7.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.29e-01

Only 100 of 262 total aligned reads displayed.

GTATGGTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCCCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACTGCAC  >  NC_000913/1434226‑1434410
                                                                                            |                                                                                            
gTATGGTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGAccc                                                                                              >  2:1967435/1‑93 (MQ=255)
  aTGGTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCa                                                                                            >  1:3704148/1‑93 (MQ=255)
    ggTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGc                                                                                          <  2:6960300/93‑1 (MQ=255)
    ggTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGc                                                                                          <  2:571952/93‑1 (MQ=255)
     gTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGcc                                                                                         <  2:6919933/93‑1 (MQ=255)
     gTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCCCAGcc                                                                                         >  2:2138691/1‑93 (MQ=255)
         aCCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCa                                                                                     >  2:2271751/1‑93 (MQ=255)
          ccTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCAc                                                                                    <  1:284052/93‑1 (MQ=255)
           cTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCAcc                                                                                   >  2:6573148/1‑93 (MQ=255)
             ggTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGt                                                                                 <  2:6600754/93‑1 (MQ=255)
              gTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGtt                                                                                <  1:7287761/93‑1 (MQ=255)
                 aTTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTc                                                                             >  1:2159827/1‑93 (MQ=255)
                 aTTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTc                                                                             <  1:1123330/93‑1 (MQ=255)
                  ttGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCa                                                                            <  2:3933482/93‑1 (MQ=255)
                  ttGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCa                                                                            >  2:3693885/1‑93 (MQ=255)
                  ttGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCa                                                                            <  1:378398/93‑1 (MQ=255)
                    ggTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATg                                                                          <  1:7378748/93‑1 (MQ=255)
                    ggTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATg                                                                          <  2:5297401/93‑1 (MQ=255)
                    ggTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCCCAGCCGTCACCGTTTTCATg                                                                          >  1:1203515/1‑93 (MQ=255)
                     gTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATgg                                                                         >  1:903791/1‑93 (MQ=255)
                     gTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATgg                                                                         <  1:1701096/93‑1 (MQ=255)
                      tGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGc                                                                        >  1:3036110/1‑93 (MQ=255)
                          gTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAAc                                                                    >  1:308/1‑93 (MQ=255)
                           tCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACa                                                                   >  2:4503987/1‑93 (MQ=255)
                            cAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACAt                                                                  <  1:6894710/93‑1 (MQ=255)
                            cAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACAt                                                                  >  1:65447/1‑93 (MQ=255)
                              gaagaGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCa                                                                <  1:7509334/93‑1 (MQ=255)
                              gaagaGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCa                                                                >  2:854576/1‑93 (MQ=255)
                              gaagaGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCa                                                                >  1:6250197/1‑93 (MQ=255)
                                 gagGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGt                                                             >  2:7347738/1‑93 (MQ=255)
                                 gagGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCCCAGCCGTCACCGTTTTCATGGCGAACATCACGt                                                             >  2:2648891/1‑93 (MQ=255)
                                    gtgtATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCg                                                          >  1:7022544/1‑93 (MQ=255)
                                     tgtATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCgt                                                         <  1:5594711/93‑1 (MQ=255)
                                      gtATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCgtt                                                        <  1:1585441/93‑1 (MQ=255)
                                       tATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCgttg                                                       <  2:1572020/93‑1 (MQ=255)
                                         tGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCgttgtt                                                     <  2:625542/93‑1 (MQ=255)
                                           ccGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTgg                                                   <  2:4310139/93‑1 (MQ=255)
                                            cGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGt                                                  <  2:1271962/93‑1 (MQ=255)
                                             gCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTg                                                 >  2:3055608/1‑93 (MQ=255)
                                               aCCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGcc                                               >  2:6308946/1‑93 (MQ=255)
                                               aCCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGcc                                               >  2:1389227/1‑93 (MQ=255)
                                                ccAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCt                                              >  2:2501442/1‑93 (MQ=255)
                                                  aGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTc                                            >  1:44192/1‑93 (MQ=255)
                                                    cgcTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCt                                          <  2:4663360/93‑1 (MQ=255)
                                                     gcTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTg                                         <  2:3784613/93‑1 (MQ=255)
                                                     gcTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTg                                         >  2:4760472/1‑93 (MQ=255)
                                                     gcTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTg                                         >  2:32305/1‑93 (MQ=255)
                                                      cTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGa                                        >  1:5705496/1‑93 (MQ=255)
                                                       tAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGAc                                       >  2:4894711/1‑93 (MQ=255)
                                                        aaaGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGAcc                                      <  2:5341138/93‑1 (MQ=255)
                                                        aaaGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGAcc                                      <  2:3572398/93‑1 (MQ=255)
                                                         aaGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCa                                     >  2:662574/1‑93 (MQ=255)
                                                         aaGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCa                                     >  2:4240923/1‑93 (MQ=255)
                                                         aaGCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCa                                     <  1:1360337/93‑1 (MQ=255)
                                                           gCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCAtt                                   >  1:2097456/1‑93 (MQ=255)
                                                           gCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCAtt                                   >  1:588577/1‑93 (MQ=255)
                                                           gCCCATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCAtt                                   <  2:2315135/93‑1 (MQ=255)
                                                             ccATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCgttgt                                                      <  1:1440768/72‑1 (MQ=255)
                                                             ccATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTAc                                 >  1:4442274/1‑93 (MQ=255)
                                                             ccATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTAc                                 <  1:2210133/93‑1 (MQ=255)
                                                              cATGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACt                                >  1:3589167/1‑93 (MQ=255)
                                                               aTGCCTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTg                               >  2:5834707/1‑93 (MQ=255)
                                                                  ccTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGGCTTCCTGACCATTACTGGCt                            >  1:7177130/1‑93 (MQ=255)
                                                                  ccTAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCt                            >  2:3868225/1‑93 (MQ=255)
                                                                    tAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTcc                          >  1:5819850/1‑93 (MQ=255)
                                                                    tAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTcc                          >  1:6054547/1‑93 (MQ=255)
                                                                    tAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTcc                          >  1:6128729/1‑93 (MQ=255)
                                                                    tAAATCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTcc                          >  1:4177158/1‑93 (MQ=255)
                                                                     aaaTCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCt                         <  1:1849188/93‑1 (MQ=255)
                                                                      anTCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCtt                        <  2:5932741/93‑1 (MQ=255)
                                                                      aaTCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCtt                        >  1:2945725/1‑93 (MQ=255)
                                                                       aTCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTc                       <  2:2997849/93‑1 (MQ=255)
                                                                       aTCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTc                       >  1:6754816/1‑93 (MQ=255)
                                                                        tCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCg                      >  2:761504/1‑93 (MQ=255)
                                                                        tCATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCg                      <  2:4987262/93‑1 (MQ=255)
                                                                         cATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgt                     >  2:3732737/1‑93 (MQ=255)
                                                                         cATATGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgt                     <  2:1157521/93‑1 (MQ=255)
                                                                          atatGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgtt                    >  2:217976/1‑93 (MQ=255)
                                                                          atatGTTGTGGAAAGACCCCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgtt                    >  1:1420577/1‑93 (MQ=255)
                                                                           tatGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgttg                   <  1:5505593/93‑1 (MQ=255)
                                                                           tatGTTGTGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgttg                   >  1:6835011/1‑93 (MQ=255)
                                                                           tatGTTGTGGAAAGACCCCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCgttg                   >  1:3208941/1‑93 (MQ=255)
                                                                               ttgtGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTAc               >  2:562620/1‑93 (MQ=255)
                                                                               ttgtGGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTAc               <  2:2747339/93‑1 (MQ=255)
                                                                                tgtgGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTAcc              <  1:4759190/93‑1 (MQ=255)
                                                                                 gtgGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCt             <  2:7310741/93‑1 (MQ=255)
                                                                                 gtgGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCt             <  1:1914096/93‑1 (MQ=255)
                                                                                  tgGAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCtt            >  2:7186673/1‑93 (MQ=255)
                                                                                    gAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGa          <  2:6396298/93‑1 (MQ=255)
                                                                                    gAAAGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGa          <  1:5351973/93‑1 (MQ=255)
                                                                                     aaaGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGAt         >  2:1500753/1‑93 (MQ=255)
                                                                                      aaGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATa        <  2:4571420/93‑1 (MQ=255)
                                                                                      aaGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATa        <  2:1879708/93‑1 (MQ=255)
                                                                                      aaGACCCCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATa        >  2:7510718/1‑93 (MQ=255)
                                                                                       aGACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATAc       >  1:4655541/1‑93 (MQ=255)
                                                                                        gACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACt      >  1:543612/1‑93 (MQ=255)
                                                                                        gACCTCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACt      >  1:6264023/1‑93 (MQ=255)
                                                                                         aCCCCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACTg     >  2:5961768/1‑93 (MQ=255)
                                                                                            ccAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACTGCAc  <  2:4069159/93‑1 (MQ=255)
                                                                                            ccAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACTGCAc  <  2:36039/93‑1 (MQ=255)
                                                                                            |                                                                                            
GTATGGTTAACCTGGTCATTGGTGCGGTCAGAAGAGGTGTATGCCGCACCAGCGCTAAAGCCCATGCCTAAATCATATGTTGTGGAAAGACCCCAGCCGTCACCGTTTTCATGGCGAACATCACGTCCGTTGTTGGTGCCTTCCTGACCATTACTGGCTCCTTCGTTGTTACCTTGATACTGCAC  >  NC_000913/1434226‑1434410

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