Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,770,087 C→A 67.4% pseudogene (55/114 nt) yibU → pseudogene, H repeat‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,770,0870CA67.4% 205.4 / 158.8 162pseudogene (55/114 nt)yibUpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (18/35);  new base A (52/57);  total (70/92)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 166 total aligned reads displayed.

ATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGT  >  NC_000913/3769988‑3770182
                                                                                                   |                                                                                               
aTGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGc                                                                                                 <  1:1136682/100‑1 (MQ=255)
aTGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGc                                                                                                 <  1:2489562/100‑1 (MQ=255)
 tGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCa                                                                                                <  2:1601975/100‑1 (MQ=255)
   gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGa                                                                                              <  1:2118121/100‑1 (MQ=255)
   gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                              >  2:2950660/1‑100 (MQ=255)
   gAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGa                                                                                              >  2:1674621/1‑100 (MQ=255)
    aaaaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGgaagaa                                                                                             <  2:1098044/100‑1 (MQ=255)
      aaGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCGTATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACa                                                                                           >  2:671982/1‑100 (MQ=255)
        gCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAg                                                                                         >  1:1283909/1‑100 (MQ=255)
         cACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGa                                                                                        >  2:265851/1‑100 (MQ=255)
           cTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGTAAGAACAAGACg                                                                                      >  2:325134/1‑100 (MQ=255)
            tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACgg                                                                                     <  1:1286998/100‑1 (MQ=255)
            tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACgg                                                                                     <  1:736176/100‑1 (MQ=255)
            tCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACgg                                                                                     <  2:2420558/100‑1 (MQ=255)
             ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGa                                                                                    >  1:2830552/1‑100 (MQ=255)
             ccTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGa                                                                                    >  1:904992/1‑100 (MQ=255)
                gCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATaa                                                                                 >  2:1990172/1‑100 (MQ=255)
                 ctctTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATaaa                                                                                >  2:1021897/1‑100 (MQ=255)
                    ttttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATaaaa                                                                               <  2:1978515/98‑1 (MQ=255)
                     tttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAAt                                                                            <  2:1000639/100‑1 (MQ=255)
                     tttttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAAt                                                                            <  2:256393/100‑1 (MQ=255)
                       tttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATcc                                                                          <  2:2056801/100‑1 (MQ=255)
                        ttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCa                                                                         <  1:1674614/100‑1 (MQ=255)
                        ttGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCa                                                                         <  2:823234/100‑1 (MQ=255)
                         tGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCaa                                                                        >  1:267377/1‑100 (MQ=255)
                           aTTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATg                                                                      <  1:1652441/100‑1 (MQ=255)
                           aTTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATg                                                                      <  1:462568/100‑1 (MQ=255)
                            tgagagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATga                                                                     >  2:2297389/3‑100 (MQ=255)
                             tAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATgan                                                                    <  1:1063377/100‑2 (MQ=255)
                              agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATgaga                                                                   <  1:3274950/100‑1 (MQ=255)
                              agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                   >  2:1784608/1‑100 (MQ=255)
                              agagCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATgaga                                                                   >  2:404266/1‑100 (MQ=255)
                                agCCGTCGCAGGAGGAACACCTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGc                                                                 <  1:227239/100‑1 (MQ=255)
                                agCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGc                                                                 <  1:1484472/100‑1 (MQ=255)
                                agCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGc                                                                 <  2:2199504/100‑1 (MQ=255)
                                   cGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGcaca                                                              >  2:3025938/1‑100 (MQ=255)
                                   cGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGcaca                                                              >  1:1146778/1‑100 (MQ=255)
                                     tCGCAGGAGGAACAACTCATATCATCATTGCGAGCTCGACGATGTATAGTCCCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCacac                                                             <  2:1362999/99‑1 (MQ=255)
                                       gCAGAAGGAACAACTCATATCATCGTTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGAAGGATAAAAAATCCAATGAGAGCACACt                                                            <  1:946195/98‑1 (MQ=255)
                                        cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAt                                                         >  1:2795785/1‑100 (MQ=255)
                                        cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAt                                                         >  2:1895733/1‑100 (MQ=255)
                                        cAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCAAATGAGAGCACACTCAt                                                         >  1:1113297/1‑100 (MQ=255)
                                         aggaggAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATc                                                        >  2:3222186/1‑100 (MQ=255)
                                         aggaggAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAAACGGATAAAAAATCCAATTAGAGCACACTCa                                                          >  1:801378/1‑98 (MQ=255)
                                          ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATcc                                                       >  2:2425575/1‑100 (MQ=255)
                                          ggaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACGCATcc                                                       >  1:1463947/1‑100 (MQ=255)
                                           gaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCa                                                      <  2:1848464/100‑1 (MQ=255)
                                           gaggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCa                                                      <  1:2117888/100‑1 (MQ=255)
                                            aggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAg                                                     <  2:2226604/100‑1 (MQ=255)
                                             ggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                    >  2:2986301/1‑100 (MQ=255)
                                             ggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                    >  1:2512507/1‑100 (MQ=255)
                                             ggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                    >  2:2815677/1‑100 (MQ=255)
                                             ggaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGa                                                    >  2:1592161/1‑100 (MQ=255)
                                              gaACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGaa                                                   <  2:2830575/100‑1 (MQ=255)
                                                acaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACt                                                 >  2:1622742/1‑100 (MQ=255)
                                                acaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACt                                                 >  2:1457426/1‑100 (MQ=255)
                                                acaacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACt                                                 <  1:504891/100‑1 (MQ=255)
                                                 caacTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACtt                                                <  2:929547/100‑1 (MQ=255)
                                                   acTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATGCAATGAGAGCACACTCATCCAGAACttat                                              >  1:2031971/1‑100 (MQ=255)
                                                    cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                             <  1:1238958/100‑1 (MQ=255)
                                                    cTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                             <  2:1822865/100‑1 (MQ=255)
                                                     tCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatta                                            <  2:2784155/100‑1 (MQ=255)
                                                     tCATATCATCATTGAGATATCGACGATGTATGGTCTCAGCACGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACttatt                                             <  1:2374242/99‑1 (MQ=255)
                                                      cATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTaa                                           >  1:374273/1‑100 (MQ=255)
                                                        tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                         >  1:1562359/1‑100 (MQ=255)
                                                        tatCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACa                                         <  1:3304334/100‑1 (MQ=255)
                                                         atcatcATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAt                                        <  2:2251843/100‑1 (MQ=255)
                                                          tcatcaTTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATa                                       >  1:3359153/1‑100 (MQ=255)
                                                           catcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACGAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATaa                                      >  2:1107413/1‑100 (MQ=255)
                                                             tcatTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAAAAAAAAAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATg                                    <  1:2949268/100‑1 (MQ=255)
                                                              catTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGa                                   <  2:80399/100‑1 (MQ=255)
                                                              catTGCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAGCCAGAACTTATTAACATAATGa                                   <  1:2424496/100‑1 (MQ=255)
                                                                  gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                               >  2:1865103/1‑100 (MQ=255)
                                                                  gCGATCTCGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAatat                               >  2:2985346/1‑100 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTaaa                           <  1:770092/99‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  1:1695876/100‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  2:584409/100‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  2:335501/100‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  2:238967/100‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  2:2330786/100‑1 (MQ=255)
                                                                       ctcGACGATGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  1:3199072/100‑1 (MQ=255)
                                                                       ctagacaaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAACACACTCATCCAGAACTTATTAACATAATGAATATTAAAg                          <  1:665562/97‑1 (MQ=255)
                                                                           acgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaa                      <  2:2319549/100‑1 (MQ=255)
                                                                            cgaTGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaaa                     <  2:2080854/100‑1 (MQ=255)
                                                                            cgaTGTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGaaaaa                     >  2:2640858/1‑100 (MQ=255)
                                                                               tGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                  <  2:832239/100‑1 (MQ=255)
                                                                               tGTATAGTCTCATCCTGCGGAAGAACAAGACAGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatca                  <  2:693613/100‑1 (MQ=255)
                                                                                gTATAGTCTCATCCTGCGGCAGCACAAGACGGATAAAAAATCCAATGCGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcat                 >  1:3152144/1‑100 (MQ=255)
                                                                                gTATAGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATTAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcat                 <  1:252390/100‑1 (MQ=255)
                                                                                 tataGTCTCATCTTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcat                 <  1:1610740/99‑1 (MQ=255)
                                                                                 tataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc                <  2:2031664/100‑1 (MQ=255)
                                                                                 taaaGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCAACCAGAACTTATTAACATAATGAATATTAAAGAAAAatcatc                <  1:838967/97‑1 (MQ=255)
                                                                                  ataGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCg               >  1:1566669/1‑100 (MQ=255)
                                                                                    aGTCTCATCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGCGCACACCCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGaa             >  1:874063/1‑100 (MQ=255)
                                                                                          aTCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAAt       <  1:1990766/100‑1 (MQ=255)
                                                                                          aTCCTGCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACAGAATGAATATTAAAGAAAAATCATCGAAGCTAAt       >  1:101673/1‑100 (MQ=255)
                                                                                            ccTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGa     >  1:150662/1‑100 (MQ=255)
                                                                                             cTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGaa    <  2:349661/100‑1 (MQ=255)
                                                                                               gCGGAAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGt  <  2:1491499/100‑1 (MQ=255)
                                                                                                  gCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGt  <  1:2209133/97‑1 (MQ=255)
                                                                                                   |                                                                                               
ATGGAAAAGCACTCCTGCTCTTTTTTGATTAGAGCCGTCGCAGGAGGAACAACTCATATCATCATTGCGATCTCGACGATGTATAGTCTCATCCTGCGGCAGAACAAGACGGATAAAAAATCCAATGAGAGCACACTCATCCAGAACTTATTAACATAATGAATATTAAAGAAAAATCATCGAAGCTAATGAAGT  >  NC_000913/3769988‑3770182

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