Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,277,706 (A)7→6 coding (370/465 nt) yhaV → toxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,277,7000A.98.5% 579.0 / 1.3 129coding (364/465 nt)yhaVtoxin of the SohB(PrlF)‑YhaV toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (74/53);  total (75/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 153 total aligned reads displayed.

GTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATT  >  NC_000913/3277602‑3277799
                                                                                                  |                                                                                                   
gTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaa                                                                                                    >  1:14556/1‑100 (MQ=255)
 tGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaa                                                                                                   >  1:1006554/1‑100 (MQ=255)
 tGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaa                                                                                                   >  1:649897/1‑100 (MQ=255)
  gCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaa                                                                                                  >  1:2482147/1‑100 (MQ=255)
  gCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaa                                                                                                  >  2:288913/1‑100 (MQ=255)
  gCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaa                                                                                                  >  1:154294/1‑100 (MQ=255)
  gCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaa                                                                                                  >  1:154370/1‑100 (MQ=255)
   cTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaa                                                                                                 <  2:3118887/100‑1 (MQ=255)
     ggTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTTTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaa                                                                                                >  1:2619922/1‑99 (MQ=255)
     ggTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaac                                                                                               <  2:2748907/100‑2 (MQ=255)
      gTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaaca                                                                                              <  1:1898257/100‑3 (MQ=255)
      gTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaaca                                                                                              <  1:730923/100‑3 (MQ=255)
      gTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaaca                                                                                              <  2:2472394/100‑3 (MQ=255)
      gTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaaca                                                                                              <  1:12819/100‑3 (MQ=255)
       tCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaacag                                                                                             <  1:3277425/100‑4 (MQ=255)
       tCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaacag                                                                                             >  2:1660503/1‑97 (MQ=255)
       tCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaacag                                                                                             <  2:2489076/100‑4 (MQ=255)
       tCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaacag                                                                                             <  2:64321/100‑4 (MQ=255)
        cGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTaaaaaacaga                                                                                            <  2:892066/100‑5 (MQ=255)
           tATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGc                                                                                        >  1:683498/1‑100 (MQ=255)
           tATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGc                                                                                        >  1:2224099/1‑100 (MQ=255)
            aTCGTCTCTTCTTTCGTTAGAGGGAAAAAGAGAAAGTCATCATGCTGGGATGGATGAACGATGAAAAGACTCTGCGCACCTAGGGTaaaaaacag                                                                                             <  2:2639695/95‑4 (MQ=255)
             tCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCt                                                                                      >  1:2709301/1‑100 (MQ=255)
              cGTCTCTTCTTTCGTTATAGTGAAAAAGCGAAAGTCATCATTCTGGTATGGATGAACGATGAAAATACTCTGCGCACCTACGGTAGAAACCAGATGCCt                                                                                      >  2:3038077/1‑99 (MQ=255)
               gTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCtat                                                                                    >  1:2847213/1‑100 (MQ=255)
               gTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCtat                                                                                    <  1:2781648/100‑1 (MQ=255)
                  tcttctTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATAcc                                                                                 <  2:1568020/100‑1 (MQ=255)
                   cttcttTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCg                                                                                <  1:1126316/100‑1 (MQ=255)
                      cttTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTAt                                                                             >  1:2302447/1‑100 (MQ=255)
                          cGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAg                                                                         <  2:77329/100‑1 (MQ=255)
                          cGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAg                                                                         <  2:3216345/100‑1 (MQ=255)
                           gTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGc                                                                        >  2:3277042/1‑100 (MQ=255)
                             tataGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTTCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCaa                                                                      <  2:3339187/100‑1 (MQ=255)
                             tataGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCATCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCaa                                                                      <  2:2077876/100‑1 (MQ=255)
                               taGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCaaaa                                                                    <  1:2881280/100‑1 (MQ=255)
                               taGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCaaaa                                                                    >  2:2430550/1‑100 (MQ=255)
                                 gTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAAGACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCaaaacg                                                                  >  2:1035319/1‑98 (MQ=255)
                                 gTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAAt                                                                   >  2:662552/1‑99 (MQ=255)
                                 gTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATg                                                                  >  2:2217612/1‑100 (MQ=255)
                                  tGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGt                                                                 >  2:61512/1‑100 (MQ=255)
                                   gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGtt                                                                <  1:2695248/100‑1 (MQ=255)
                                   gAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGtt                                                                >  2:1008113/1‑100 (MQ=255)
                                    aaaaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTa                                                               >  1:2169241/1‑100 (MQ=255)
                                    aaaaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTa                                                               >  1:191290/1‑100 (MQ=255)
                                     aaaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaa                                                              <  1:380155/100‑1 (MQ=255)
                                     aaaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaa                                                              >  2:2795598/1‑100 (MQ=255)
                                      aaaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaaa                                                             >  1:3009238/1‑100 (MQ=255)
                                       aaGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaaaa                                                            <  1:3288455/100‑1 (MQ=255)
                                        agagAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTaaaaa                                                           >  2:1472746/1‑100 (MQ=255)
                                         gagaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAag                                                          >  2:3345822/1‑100 (MQ=255)
                                         gagaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAag                                                          >  1:166467/1‑100 (MQ=255)
                                         gagaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAag                                                          >  1:1012140/1‑100 (MQ=255)
                                          agaAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACAGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAaga                                                         <  2:649896/100‑1 (MQ=255)
                                            aaaGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAgg                                                       <  1:309628/100‑1 (MQ=255)
                                              aGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGAc                                                     >  2:3309914/1‑100 (MQ=255)
                                              aGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGAc                                                     <  1:2963870/100‑1 (MQ=255)
                                               gTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACa                                                    <  1:1660507/100‑1 (MQ=255)
                                                tcatcaTTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACAt                                                   >  2:157549/1‑100 (MQ=255)
                                                 catcatTCTGGGACGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACAtc                                                  >  1:341449/1‑100 (MQ=255)
                                                  atcatTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACAtcc                                                 >  2:447116/1‑100 (MQ=255)
                                                       tCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTg                                            >  1:1290341/1‑100 (MQ=255)
                                                        cTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGc                                           <  1:546896/100‑1 (MQ=255)
                                                         tcggatggatGAACGATGAAAACACTCTGCGAACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGcc                                          <  2:2038641/98‑1 (MQ=255)
                                                           ggatggatGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGa                                        >  2:1131517/1‑100 (MQ=255)
                                                            gatggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGAc                                       <  1:860935/100‑1 (MQ=255)
                                                             atggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAGCAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACt                                      <  2:2352197/100‑1 (MQ=255)
                                                              tggatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTg                                     >  1:290483/1‑100 (MQ=255)
                                                                gatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTggg                                   >  1:96908/1‑100 (MQ=255)
                                                                gatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTggg                                   >  1:2182142/1‑100 (MQ=255)
                                                                gatgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTgg                                    >  2:678960/1‑99 (MQ=255)
                                                                 atgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGa                                  >  2:1861580/1‑100 (MQ=255)
                                                                 atgAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGa                                  <  1:1904557/100‑1 (MQ=255)
                                                                   gAACGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGaan                                <  1:1267347/100‑2 (MQ=255)
                                                                    aaCGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAc                               <  2:3181840/100‑1 (MQ=255)
                                                                    aaCGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAc                               >  1:1761874/1‑100 (MQ=255)
                                                                    aaCGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAc                               >  1:2334362/1‑100 (MQ=255)
                                                                     aCGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAcc                              >  2:1492273/1‑100 (MQ=255)
                                                                     aCGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAcc                              <  1:308056/100‑1 (MQ=255)
                                                                      cGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAccc                             >  2:3102415/1‑100 (MQ=255)
                                                                      cGATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAccc                             >  1:380412/1‑100 (MQ=255)
                                                                      cGATGAAAACACTCTGCGAACCAACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAccn                             <  2:2183428/100‑2 (MQ=255)
                                                                      cGATGAAAACACTCGGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAAccc                             >  1:624344/1‑100 (MQ=255)
                                                                       gATGAAAACACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCt                            >  2:3070252/1‑100 (MQ=255)
                                                                           aaaaCACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCAcc                        <  1:2218957/100‑1 (MQ=255)
                                                                           aaaaCACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCAcc                        >  2:1743676/1‑100 (MQ=255)
                                                                            aaaCACTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCAccc                       >  1:1754965/1‑100 (MQ=255)
                                                                               cacTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCgag                    <  1:1008115/100‑1 (MQ=255)
                                                                                acTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCgaga                   <  1:447117/100‑1 (MQ=255)
                                                                                acTCTGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCgaga                   <  1:2986598/100‑1 (MQ=255)
                                                                                 ctctGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGaa                  <  2:2302436/100‑1 (MQ=255)
                                                                                  tctGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGaaa                 <  2:2334351/100‑1 (MQ=255)
                                                                                   ctGCGCACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAAc                <  2:191289/100‑1 (MQ=255)
                                                                                       gcACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACagaa            >  1:2527619/1‑100 (MQ=255)
                                                                                        cACCTACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACagaag           <  2:290482/100‑1 (MQ=255)
                                                                                            tACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAc       >  2:3040738/1‑100 (MQ=255)
                                                                                            tACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAc       >  1:2498497/1‑100 (MQ=255)
                                                                                            tACGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAc       >  2:906658/1‑100 (MQ=255)
                                                                                             aCGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAcc      <  2:380411/100‑1 (MQ=255)
                                                                                             aCGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAcc      >  2:602620/1‑100 (MQ=255)
                                                                                             aCGGTAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAAcc      >  1:335467/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTGCTGGTCGTTATCGTCTCTTCTTTCGTTATAGTGAAAAAGAGAAAGTCATCATTCTGGGATGGATGAACGATGAAAACACTCTGCGCACCTACGGTAAAAAAACAGATGCCTATACCGTATTCAGCAAAATGTTAAAAAGAGGACATCCTCCTGCCGACTGGGAAACCCTCACCCGAGAAACAGAAGAAACCCATT  >  NC_000913/3277602‑3277799

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