Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,707,989 A→G 43.6% intergenic (‑284/+627) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,707,9890AG43.6% 60.9 / 207.9 175intergenic (‑284/+627)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base A (49/49);  new base G (42/34);  total (91/84)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

TGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAA  >  NC_000913/3707893‑3708087
                                                                                                 |                                                                                                  
tgtgAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAAt                                                                                                 >  2:3904787/1‑100 (MQ=255)
 gtgAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAAtt                                                                                                >  1:502643/1‑100 (MQ=255)
   gAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTtg                                                                                              >  2:852414/1‑100 (MQ=255)
   gAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtg                                                                                              <  2:962392/100‑1 (MQ=255)
   gAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtg                                                                                              <  2:2650464/100‑1 (MQ=255)
   gAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtg                                                                                              <  2:1471088/100‑1 (MQ=255)
    aaaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTtgt                                                                                             <  1:4230814/100‑1 (MQ=255)
    aaaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgt                                                                                             <  1:1056251/100‑1 (MQ=255)
    aaaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgt                                                                                             >  2:3720276/1‑100 (MQ=255)
    aaaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgt                                                                                             <  1:3202480/100‑1 (MQ=255)
     ttccgatcttACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTtgt                                                                                             <  2:2984472/91‑1 (MQ=255)
     aaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTtgtg                                                                                            <  2:3315192/100‑1 (MQ=255)
     aaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTtgtg                                                                                            <  1:2221243/100‑1 (MQ=255)
     aaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTtgtg                                                                                            <  2:1887582/100‑1 (MQ=255)
     aaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgtg                                                                                            <  2:1163836/100‑1 (MQ=255)
     aaaCGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTtgtg                                                                                            <  2:1008162/100‑1 (MQ=255)
      aaCGTCAATACAGCCAACCGGGATTTACACTAGCGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGc                                                                                           <  1:2449870/100‑1 (MQ=255)
        ngTCAATACAGCCAACCGGAATTTACACCAACAGTGAGAATCCACACACAAAGATTAAAAAAACTCCAAACAGCACTTTGCCGCAAAGAACATTGTTCTa                                                                                         >  1:3601937/2‑100 (MQ=255)
         gTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTAc                                                                                        <  2:2177313/100‑1 (MQ=255)
         gTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTAc                                                                                        <  2:116740/100‑1 (MQ=255)
          tCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACt                                                                                       <  1:1226088/100‑1 (MQ=255)
           cAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTc                                                                                      >  2:696455/1‑100 (MQ=255)
            aaTACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTcc                                                                                     <  1:475868/100‑1 (MQ=255)
             aTACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCa                                                                                    >  1:1350417/1‑100 (MQ=255)
              tACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCaa                                                                                   >  1:3950794/1‑100 (MQ=255)
              tACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCaa                                                                                   <  1:3031686/100‑1 (MQ=255)
                cAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACa                                                                                 >  2:2275417/1‑100 (MQ=255)
                  gCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATg                                                                               <  1:1089467/100‑1 (MQ=255)
                   ccAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGa                                                                              >  2:3581501/1‑100 (MQ=255)
                     tccgatcttttACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGAcc                                                                            <  1:2911395/92‑1 (MQ=255)
                       ncGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAg                                                                          >  1:3404627/2‑100 (MQ=255)
                       ccGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAg                                                                          >  1:3806582/1‑100 (MQ=255)
                         gggATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGaa                                                                        >  2:2229423/1‑100 (MQ=255)
                          ggATTTACACCAACGATGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAAc                                                                       <  2:4204461/100‑1 (MQ=255)
                               tACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCa                                                                  <  1:2275411/100‑1 (MQ=255)
                                acacCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAg                                                                 <  2:1616600/100‑1 (MQ=255)
                                acacCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAg                                                                 <  2:1282774/100‑1 (MQ=255)
                                 cacCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGc                                                                >  2:2654676/1‑100 (MQ=255)
                                 cacCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGc                                                                >  2:4146019/1‑100 (MQ=255)
                                   ccAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCtt                                                              >  2:1920712/1‑100 (MQ=255)
                                   ccAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTACAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACagatcggaa                                                              >  2:1760746/1‑92 (MQ=255)
                                      aCGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAat                                                           <  1:377885/100‑1 (MQ=255)
                                        ggTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAatat                                                         >  1:85550/1‑100 (MQ=255)
                                        ggTGAGAATCAAAACACAAAGATTAAAAAAACTTCAGACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAatat                                                         <  1:1455441/100‑1 (MQ=255)
                                          tGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCATAACAATCAGCTTAATAttt                                                       <  2:1049912/100‑1 (MQ=255)
                                           gagaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTa                                                      <  1:237014/100‑1 (MQ=255)
                                            agaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAg                                                     <  1:1732251/100‑1 (MQ=255)
                                            agaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATATGCAGCAGAGGAATTTGTGCTAATCCAACATGACCAGAACAATCAGCTTAATATTTAg                                                     >  1:2090789/1‑100 (MQ=255)
                                             gaATCCATACAAAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGc                                                    <  1:3160897/100‑1 (MQ=255)
                                             gaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGc                                                    <  2:4081194/100‑1 (MQ=255)
                                             gaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGc                                                    <  2:3966238/100‑1 (MQ=255)
                                             gaATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGc                                                    >  1:1840690/1‑100 (MQ=255)
                                              aaTCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTAGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGc                                                    <  2:3841262/99‑1 (MQ=255)
                                                 ccACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACa                                                <  2:3137376/100‑1 (MQ=255)
                                                  cacacacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACAt                                               >  1:1966203/1‑100 (MQ=255)
                                                  cacacacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACAt                                               >  2:29365/1‑100 (MQ=255)
                                                      cacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGATCAGAACAATCAGCTTAATATTTAGCAACATTTgg                                           >  2:3517746/1‑100 (MQ=255)
                                                      cacaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTgg                                           <  1:338642/100‑1 (MQ=255)
                                                       acaAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGt                                          <  1:1047052/100‑1 (MQ=255)
                                                        ctaaGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTg                                         <  1:1903833/98‑1 (MQ=255)
                                                        caAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTg                                         <  2:3103339/100‑1 (MQ=255)
                                                        caAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTg                                         >  2:145445/1‑100 (MQ=255)
                                                         aaaGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTg                                         <  1:391450/99‑1 (MQ=255)
                                                         aaaGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGa                                        >  2:427671/1‑100 (MQ=255)
                                                             aTTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTaa                                    >  2:540370/1‑100 (MQ=255)
                                                             aTTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTaa                                    >  2:1735763/1‑100 (MQ=255)
                                                              ttAAAAAAACTTCAAACATCTATTTGAAGCAGAGAAATATATGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAAGATTTGGTGATTaa                                    >  1:998404/1‑99 (MQ=255)
                                                              ttAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAg                                   <  2:2718258/100‑1 (MQ=255)
                                                               tAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGt                                  >  2:207563/1‑100 (MQ=255)
                                                               tAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGt                                  >  2:2261108/1‑100 (MQ=255)
                                                                 aaaaaaaCTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGtt                                 >  1:3015500/1‑100 (MQ=255)
                                                                  aaaaaaCTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGttt                                >  2:2850105/1‑100 (MQ=255)
                                                                  aaaaaaCTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACTTTTGGTGATTAAGttt                                >  1:1328071/1‑100 (MQ=255)
                                                                  aaaaaaCTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGttt                                >  2:169744/1‑100 (MQ=255)
                                                                  aaaaaaCTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGttt                                >  1:1295176/1‑100 (MQ=255)
                                                                     aaaCTTCAAACAGCTATTTTCAGCTGAGGAATTTGTTCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTa                              <  1:337823/99‑1 (MQ=255)
                                                                      aaCTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATg                            <  2:3444807/100‑1 (MQ=255)
                                                                      aaCTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATg                            <  1:2984661/100‑1 (MQ=255)
                                                                       aattCAAACAGCTATTTGCAGCAGAGAAATTAGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGc                           >  2:2879035/3‑100 (MQ=255)
                                                                       aCTTCAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGc                           <  2:658329/100‑1 (MQ=255)
                                                                           cAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCAttt                       >  2:3948919/1‑100 (MQ=255)
                                                                           cAAACAGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCAttt                       <  2:439962/100‑1 (MQ=255)
                                                                           cAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCAttt                       >  1:328021/1‑100 (MQ=255)
                                                                             aaaaGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTaa                     >  1:1843849/4‑100 (MQ=255)
                                                                             aaCAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTaa                     >  1:1460844/1‑100 (MQ=255)
                                                                              aCAGCTATTTGCAGCAGAGCAAAAAGTACTCCTCCAACATGACAACCACCATAAGCTTAATCTTTAGCAACATTAGGTGATTAAGTTTTATGCATTTAAt                    <  1:353926/100‑1 (MQ=255)
                                                                               cAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATg                   <  1:2682860/100‑1 (MQ=255)
                                                                               cAGCTATTTGCAGCAGAGAAACTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATg                   >  1:364730/1‑100 (MQ=255)
                                                                                aGCTATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGa                  <  1:2228658/100‑1 (MQ=255)
                                                                                   tATTTGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGaaaa               <  2:1114949/100‑1 (MQ=255)
                                                                                   tATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGaaaa               >  1:3753270/1‑100 (MQ=255)
                                                                                   tATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGaaaa               <  2:910678/100‑1 (MQ=255)
                                                                                    actcgcaacaGAGAAATACGTGCTACTCCAACATGATCAGAACAATCAGCTTAATATTTAGCACTATTCGTTGATCAATATTTATGCATTTAATGaaaaa              >  1:1117250/5‑100 (MQ=255)
                                                                                       tGCAGCAGAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAAt           <  1:1511524/100‑1 (MQ=255)
                                                                                        gcagcaGAGGCATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATc          >  1:2309872/1‑100 (MQ=255)
                                                                                         cagcagAGAAATTTGTGCTACTCCAACATGACCATAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATc          <  2:39046/99‑1 (MQ=255)
                                                                                            cagAGGAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTg        <  1:161020/98‑1 (MQ=255)
                                                                                            cagAGAAATTTGTGCTACTCCAACATGACCGGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAg      <  1:1374237/100‑1 (MQ=255)
                                                                                            cagAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGAAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAg      <  1:4029210/100‑1 (MQ=255)
                                                                                              gagaAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGa    <  1:4097569/100‑1 (MQ=255)
                                                                                                 |                                                                                                  
TGTGAAAACGTCAATACAGCCAACCGGGATTTACACCAACGGTGAGAATCCACACACAAAGATTAAAAAAACTTCAAACAGCTATTTGCAGCAGAGAAATTTGTGCTACTCCAACATGACCAGAACAATCAGCTTAATATTTAGCAACATTTGGTGATTAAGTTTTATGCATTTAATGAAAAAAATCTGAGGAAA  >  NC_000913/3707893‑3708087

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