Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,426,178 A→T 16.7% intergenic (+96/‑276) ydaY → / → tmpR pseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,426,1780AT16.7% 313.5 / 50.5 125intergenic (+96/‑276)ydaY/tmpRpseudogene, Rac prophage;Phage or Prophage Related/Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base A (49/55);  new base T (13/8);  total (62/63)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

Only 100 of 126 total aligned reads displayed.

TAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAG  >  NC_000913/1426080‑1426276
                                                                                                  |                                                                                                  
tAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                                                   <  1:1155529/100‑1 (MQ=255)
 aGACGCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATGCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAc                                                                                                   <  2:1011315/99‑1 (MQ=255)
  ggcgcGGTTCTAAAGGAGTTGTTGTACGGCTATTCGGTAAGACCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTATATCTCTGACag                                                                                                 <  2:1544328/96‑1 (MQ=255)
    cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                               <  2:1275575/100‑1 (MQ=255)
    cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                               <  1:450826/100‑1 (MQ=255)
    cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                               <  1:86217/100‑1 (MQ=255)
    cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                               >  2:1372246/1‑100 (MQ=255)
    cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                               >  1:2492203/1‑100 (MQ=255)
    cccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTCATCTCTGACagag                                                                                               >  2:2214594/1‑100 (MQ=255)
     ccGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACagag                                                                                               <  1:3032531/99‑1 (MQ=255)
       ggTGCTTAAGTAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCAGGCGTGCCTTTGCCTTTATCTCTGACAGAGATa                                                                                            <  2:2048598/100‑1 (MQ=255)
          gCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaag                                                                                         <  2:1207590/100‑2 (MQ=255)
           cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                        <  1:670108/100‑3 (MQ=255)
           cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                        <  1:2584549/100‑3 (MQ=255)
           cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                        <  2:2807970/100‑3 (MQ=255)
           cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaagg                                                                                        <  2:668392/100‑3 (MQ=255)
           cTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAAtn                                                                                        <  2:591105/100‑2 (MQ=255)
            ttAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaaggc                                                                                       <  2:848404/100‑4 (MQ=255)
            ttAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                       >  1:1071188/1‑100 (MQ=255)
            ttAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                       >  2:635856/1‑100 (MQ=255)
            ttAAGGAGTGGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGTCAGAGATaaaggc                                                                                       <  1:818408/100‑4 (MQ=255)
            ntAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGa                                                                                       >  2:427503/2‑100 (MQ=255)
              aaGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAat                                                                                     >  2:789130/1‑100 (MQ=255)
               aGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAata                                                                                    >  1:2132249/1‑100 (MQ=255)
               aGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAata                                                                                    <  1:1270905/100‑1 (MQ=255)
                   gttgttGTACGCCTATTCGGTAGTATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTa                                                                                >  1:234151/1‑100 (MQ=255)
                    ttgttgTACTCCTATTCGGTAATGTCCCGTGCCCGACGTTATGCTGGTATGTCTGGGGTGCCTTTGCTTTTATCTCTGACAGAGATAAATGAATATTTAg                                                                               >  2:1230939/1‑100 (MQ=255)
                     tgttgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGc                                                                              >  2:1839088/1‑100 (MQ=255)
                        tgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGCATATTTAGcc                                                                             >  1:945817/1‑98 (MQ=255)
                        tgtACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGcctc                                                                           <  1:2107719/100‑3 (MQ=255)
                             gCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATc                                                                      >  2:1402502/1‑100 (MQ=255)
                              ccTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATcc                                                                     >  2:1503440/1‑100 (MQ=255)
                                tATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCgg                                                                   >  2:27870/1‑100 (MQ=255)
                                tATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTGGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCgg                                                                   <  1:986488/100‑1 (MQ=255)
                                 aTTCGGTAATCTCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGt                                                                  >  2:2967236/1‑100 (MQ=255)
                                  ttCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTa                                                                 >  2:2723538/1‑100 (MQ=255)
                                    cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTatt                                                               <  1:1105183/100‑1 (MQ=255)
                                    cGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGAAAGAAATACATGAAGATTTAGCCACTCATCCGGt                                                                  >  2:2552530/1‑97 (MQ=255)
                                     ggTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattg                                                              <  2:2957786/100‑1 (MQ=255)
                                         atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                          <  2:2748332/100‑1 (MQ=255)
                                         atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                          <  1:2888348/100‑1 (MQ=255)
                                         atatCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTattgattg                                                          <  2:1023856/100‑1 (MQ=255)
                                              ccGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAgcg                                                      >  2:2247482/1‑99 (MQ=255)
                                               cGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGcgcg                                                    <  2:2464185/100‑1 (MQ=255)
                                                 tGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGAt                                                  <  1:1674692/100‑1 (MQ=255)
                                                  gCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATg                                                 <  1:575639/100‑1 (MQ=255)
                                                   cccGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGa                                                >  2:75030/1‑100 (MQ=255)
                                                     cGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAAt                                              >  2:338948/1‑100 (MQ=255)
                                                        cGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTg                                           <  1:3071554/100‑1 (MQ=255)
                                                         gTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGa                                          >  2:965263/1‑100 (MQ=255)
                                                         gTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGa                                          >  2:1289137/1‑100 (MQ=255)
                                                           tATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                        <  2:1771853/100‑1 (MQ=255)
                                                           tATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAg                                        <  1:965267/100‑1 (MQ=255)
                                                             tGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                      >  1:1584699/1‑100 (MQ=255)
                                                             tGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCa                                      >  2:525544/1‑100 (MQ=255)
                                                               cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGt                                    <  2:1567300/100‑1 (MQ=255)
                                                               cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGt                                    >  1:372839/1‑100 (MQ=255)
                                                               cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGt                                    <  1:467107/100‑1 (MQ=255)
                                                               cTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCCGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGt                                    >  1:2227737/1‑100 (MQ=255)
                                                                tGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGtn                                   <  2:1268257/100‑2 (MQ=255)
                                                                   aaTGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATc                                <  2:2257005/100‑1 (MQ=255)
                                                                     tGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCtt                              <  1:760055/100‑1 (MQ=255)
                                                                       gCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTg                            >  2:2007756/1‑100 (MQ=255)
                                                                         taggatgcctttgccatTATCTCTGACAGAAACAAATGAATACTTACCCACTAATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGc                           <  2:836406/94‑1 (MQ=255)
                                                                          ggggTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCAc                         <  1:800521/100‑1 (MQ=255)
                                                                            ggTGCCTTTGCCTTTATCCCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTg                       <  1:1857579/100‑1 (MQ=255)
                                                                              tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGa                     <  1:121314/100‑1 (MQ=255)
                                                                              tgcctttgccttTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGa                     <  1:635857/100‑1 (MQ=255)
                                                                                cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATg                   <  2:204994/100‑1 (MQ=255)
                                                                                cctttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATg                   <  1:437626/100‑1 (MQ=255)
                                                                                  tttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAc                 <  2:940005/100‑1 (MQ=255)
                                                                                  tattcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATg                   >  2:550169/7‑98 (MQ=255)
                                                                                   ttgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTCAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                >  2:2184068/1‑100 (MQ=255)
                                                                                   attcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                >  2:882884/6‑100 (MQ=255)
                                                                                   attcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGAcc                >  1:2246327/6‑100 (MQ=255)
                                                                                    ttcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               >  2:540658/5‑100 (MQ=255)
                                                                                    ttcatttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               >  1:2132938/5‑100 (MQ=255)
                                                                                    ttcatttATCTCTGTCAGAGATAAATGAATATTTAACCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               >  2:1305870/5‑100 (MQ=255)
                                                                                    ttcatttATCTCTGTCAGAGATAAAGGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               >  2:666659/5‑100 (MQ=39)
                                                                                    tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               <  2:3004408/100‑1 (MQ=255)
                                                                                    tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               <  1:1250915/100‑1 (MQ=255)
                                                                                    tgcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCa               >  1:2338952/1‑100 (MQ=255)
                                                                                     gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAg              <  1:655040/100‑1 (MQ=255)
                                                                                     gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAg              <  1:1288804/100‑1 (MQ=255)
                                                                                     gcctttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAg              >  1:132681/1‑100 (MQ=255)
                                                                                      catttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt             >  2:372094/3‑100 (MQ=255)
                                                                                      catttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt             >  2:1309970/3‑100 (MQ=255)
                                                                                       atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGt             >  2:1741570/2‑99 (MQ=255)
                                                                                       atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTa            >  2:1041399/2‑100 (MQ=255)
                                                                                       atttATCTCTGTCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTa            >  2:1042435/2‑100 (MQ=255)
                                                                                        tttATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAt           <  2:517102/100‑1 (MQ=255)
                                                                                          tATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAttt         >  1:2011308/1‑100 (MQ=255)
                                                                                             ctctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTAttttga      >  2:1886127/1‑98 (MQ=255)
                                                                                             ctctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCa      <  1:1503451/100‑1 (MQ=255)
                                                                                              tctGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCAATGAATTTGAAGCAGAGATCTTTGGACTGGATGACCAGTATTTTCa      >  2:2097527/1‑99 (MQ=255)
                                                                                               ntGAAAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAgg    >  2:2632417/2‑100 (MQ=255)
                                                                                               ctaaCAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACGAGTATTTTCAg     >  1:851815/4‑99 (MQ=255)
                                                                                               ctGACAGAGATAGATGAATATTTGGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGCAGCAGTGATCGGTGCACTGGATGACCAGTATTTTCAgg    <  2:395709/100‑1 (MQ=255)
                                                                                                 gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAg  >  2:2702708/1‑100 (MQ=255)
                                                                                                 gACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAg  <  1:357738/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
TAGACCCGGTGCTTAAGGAGTTGTTGTACGCCTATTCGGTAATATCCCGTGCCCGACGTTATGCTGGAATGGCTGGGGTGCCTTTGCCTTTATCTCTGACAGAGATAAATGAATATTTAGCCACTCATCCGGTATTGATTGAGCGCGATGAATTTGAAGCAGTGATCTTTGCACTGGATGACCAGTATTTTCAGGAG  >  NC_000913/1426080‑1426276

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: