Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,758,577 G→T N8K (AAC→AAA)  yfjH ← CP4‑57 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,758,5770GT98.0% 487.3 / ‑2.0 148N8K (AAC→AAA) yfjHCP4‑57 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (66/78);  total (69/79)
Fisher's exact test for biased strand distribution p-value = 5.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 154 total aligned reads displayed.

TTTTCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATAC  >  NC_000913/2758478‑2758671
                                                                                                   |                                                                                              
ttttCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATggg                                                                                                >  2:1494409/1‑100 (MQ=255)
   tCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttt                                                                                             <  2:1525522/100‑1 (MQ=255)
     tACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGtttttt                                                                                           >  2:831795/1‑100 (MQ=255)
     tACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttt                                                                                            >  1:362478/1‑99 (MQ=255)
      aCAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                          <  1:2861365/100‑1 (MQ=255)
      aCAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                          >  2:244906/1‑100 (MQ=255)
      aCAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                          >  1:221600/1‑100 (MQ=255)
      aCAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTtttttt                                                                                          >  2:3301706/1‑100 (MQ=255)
       cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                         >  2:1451756/1‑100 (MQ=255)
       cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                         <  2:3329725/100‑1 (MQ=255)
       cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                         >  2:480318/1‑100 (MQ=255)
       cAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTg                                                                                         <  1:633446/100‑1 (MQ=255)
        aaaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGc                                                                                        <  2:327355/100‑1 (MQ=255)
         aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcn                                                                                       <  1:3057803/100‑2 (MQ=255)
         aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcc                                                                                       <  2:1598828/100‑1 (MQ=255)
         aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcc                                                                                       <  1:313608/100‑1 (MQ=255)
         aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGcc                                                                                       <  2:2985077/100‑1 (MQ=255)
         aaTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACACGGTTTTTTTGcc                                                                                       <  2:2283540/100‑1 (MQ=255)
          aTTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCt                                                                                      <  1:1522552/100‑1 (MQ=255)
           tttGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGTCtt                                                                                     >  1:2746925/1‑100 (MQ=255)
              gTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGttttttagatcggaa                                                                                  >  2:2366097/1‑91 (MQ=255)
              gTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTaaa                                                                                  >  1:3407356/1‑100 (MQ=255)
               ttCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAg                                                                                 <  2:2591691/100‑1 (MQ=255)
                tCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGa                                                                                <  1:1347230/100‑1 (MQ=255)
                 cATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTTTTTTTTTTTTTTTTTGCATTAAAGAt                                                                               >  1:1381928/1‑100 (MQ=255)
                  aTAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                              >  1:101484/1‑100 (MQ=255)
                  aTAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                              >  1:101524/1‑100 (MQ=255)
                  aTAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                              <  2:2305245/100‑1 (MQ=255)
                   tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATa                                                                              >  1:3013522/1‑99 (MQ=255)
                   tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                             >  2:1814870/1‑100 (MQ=255)
                   tAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAc                                                                             >  2:3152621/1‑100 (MQ=255)
                    aaaaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAcc                                                                            <  2:2807585/100‑1 (MQ=255)
                    aaaaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAcc                                                                            <  2:3038744/100‑1 (MQ=255)
                    aaaaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATAcc                                                                            <  2:549099/100‑1 (MQ=255)
                     aaaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTCCTTCAAGATa                                                                              >  1:2353406/1‑97 (MQ=255)
                      aaaTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCa                                                                          <  1:2677179/100‑1 (MQ=255)
                        aTCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCAtt                                                                        <  1:3453404/100‑1 (MQ=255)
                           aTTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAt                                                                     >  2:1972655/1‑100 (MQ=255)
                            ttCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtc                                                                    >  1:1263484/1‑100 (MQ=255)
                              ccgagcaaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAACTTGTATTTTTAGTACATGGTTTTTTTGCCATAAATATATCCCTTATtc                                                                    >  1:616035/7‑98 (MQ=255)
                              cAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtct                                                                   >  2:2536674/1‑99 (MQ=255)
                                gagaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctctc                                                                <  2:627685/100‑1 (MQ=255)
                                gagaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAAtctcta                                                                >  2:1336141/1‑99 (MQ=255)
                                  gaAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTg                                                              <  2:830368/100‑1 (MQ=255)
                                   aaaTCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGc                                                             <  1:1094598/100‑1 (MQ=255)
                                       cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                         >  2:367729/1‑100 (MQ=255)
                                       cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                         >  2:3096078/1‑100 (MQ=255)
                                       cATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGATACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAAcc                                                         >  2:2435369/1‑100 (MQ=255)
                                         tAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCaa                                                       >  2:470955/1‑100 (MQ=255)
                                          aaaaCCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCaaa                                                      >  2:3210300/1‑100 (MQ=255)
                                              ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTTAATTTACCCATTAATCTCTCTGCAACTAAAGt                                                    >  1:1952042/1‑98 (MQ=255)
                                              ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTg                                                   <  1:3351340/99‑1 (MQ=255)
                                              ccaccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGa                                                  <  2:1260647/100‑1 (MQ=255)
                                                accaTGAATCAATTTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                <  2:3249037/100‑1 (MQ=255)
                                                accaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAc                                                <  1:1618416/100‑1 (MQ=255)
                                                 ccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAcc                                               <  2:3439647/100‑1 (MQ=255)
                                                 ccaTAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAAcc                                               <  2:2004764/100‑1 (MQ=255)
                                                     aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATg                                           <  2:1096673/100‑1 (MQ=255)
                                                     aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATg                                           <  1:2549718/100‑1 (MQ=255)
                                                     aaaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATg                                           <  2:2038469/100‑1 (MQ=255)
                                                      aaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTTAACCAATga                                          >  1:2373142/1‑100 (MQ=255)
                                                      aaTCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATga                                          <  2:2586461/100‑1 (MQ=255)
                                                          aaTGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAgg                                      >  2:1067746/1‑100 (MQ=255)
                                                            tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGca                                    <  2:2417797/100‑1 (MQ=255)
                                                            tGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGca                                    <  1:430550/100‑1 (MQ=255)
                                                              tATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGcaac                                  >  2:678175/1‑100 (MQ=255)
                                                                 tAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAg                               >  2:3449155/1‑100 (MQ=255)
                                                                 tAGAGGTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAg                               <  2:582556/100‑1 (MQ=255)
                                                                     ggTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTTTTTTTAACCATTAATCTCTCTGCAACCAACGTGAACCAATGAGAGGCAACAAGAAt                            >  1:3447024/1‑99 (MQ=255)
                                                                     ggTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATg                           <  1:3103122/100‑1 (MQ=255)
                                                                      gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATg                           <  2:1613290/99‑1 (MQ=255)
                                                                      gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGa                          >  2:2542888/1‑100 (MQ=255)
                                                                      gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGa                          <  1:329695/100‑1 (MQ=255)
                                                                      gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGa                          <  1:326712/100‑1 (MQ=255)
                                                                      gTTGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGa                          <  1:1582470/100‑1 (MQ=255)
                                                                       ttGTAAAAATTGTATTTTTGGTACATTGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAATTGAACCAATGAGAGGCAACAAGAATGa                          <  1:3259862/99‑1 (MQ=255)
                                                                        tGTAAAAATTGTATTTTTGGTACATGGTTTTTTTTTTTTTTTTTTACCTATTAATCTCTCTTCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                        >  2:2867562/1‑100 (MQ=255)
                                                                        tGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCACAGTGAACCAATGAGAGGCAACAAGAATGAAc                        <  1:1552505/100‑1 (MQ=255)
                                                                        tGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                        >  2:1871150/1‑100 (MQ=255)
                                                                        tGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                        <  1:3349022/100‑1 (MQ=255)
                                                                        tGTAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAAc                        <  1:1698763/100‑1 (MQ=255)
                                                                          tAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACaa                      >  2:3232403/1‑100 (MQ=255)
                                                                          tAAAAATTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACaa                      >  1:2823347/1‑100 (MQ=255)
                                                                            aaaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                    <  1:1945556/100‑1 (MQ=255)
                                                                            aaaaTTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACTAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatc                    >  1:3298097/1‑100 (MQ=255)
                                                                             aaaTTGTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAatcg                   <  1:1218898/100‑1 (MQ=255)
                                                                                  gTATTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCgg              <  2:2061090/100‑1 (MQ=255)
                                                                                    aTTTTTGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGcc            >  1:823814/1‑100 (MQ=255)
                                                                                    aTTTTTGGTACATGGTTTTTTTGCCTTAAAGATAAACATTAATCTCGCTGCAACCAAAGTGAACAAATGAGAGGCAACAAGAATGAACAATCGATCGGc             >  1:855004/1‑99 (MQ=255)
                                                                                     tttttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCg           <  2:3375/100‑1 (MQ=255)
                                                                                     tttttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCg           <  2:713688/100‑1 (MQ=255)
                                                                                     tttttGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCg           <  1:1494387/100‑1 (MQ=255)
                                                                                     cctttGGTATATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCg           <  2:2415154/98‑1 (MQ=255)
                                                                                      ttttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGt          >  2:746388/1‑100 (MQ=255)
                                                                                       tttGGTACATGGTTTTTTTACTTTAAAAATCCCCATAAATCTCCCTTTATCAAAATTGAACAAATGAAAGGAAAAAAGAATAAACAATCGATCGGCCgag         >  1:2415121/1‑98 (MQ=255)
                                                                                        ttGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTa        >  1:2059737/1‑100 (MQ=255)
                                                                                         tGGTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAg       <  2:974727/100‑1 (MQ=255)
                                                                                          ggTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGa      <  1:2045446/100‑1 (MQ=255)
                                                                                           gTACATGGTTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGaa     <  1:1106115/100‑1 (MQ=255)
                                                                                              cATGGTTTTTTTTTCTTTTATATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAAtac  >  1:1884286/1‑100 (MQ=255)
                                                                                                   |                                                                                              
TTTTCTACAAATTTGTTCATAAAAATCATTCAGAGAAATCATAAAACCACCATAAATCAATGTATTAGAGGTTGTAAAAATTGTATTTTTGGTACATGGGTTTTTTGCCTTAAAGATACCCATTAATCTCTCTGCAACCAAAGTGAACCAATGAGAGGCAACAAGAATGAACAATCGATCGGCCGTTAGAATAC  >  NC_000913/2758478‑2758671

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: ‑