Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 736,615 T→C 7.8% pseudogene (171/240 nt) ybfQ → pseudogene, H repeat‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913736,6150TC 7.8% 182.3 / 5.5 67pseudogene (171/240 nt)ybfQpseudogene, H repeat‑associated protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base C (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCT  >  NC_000913/736469‑736764
                                                                                                                                                  |                                                                                                                                                     
cATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTg                                                                                                                                                    <  1:1293580/150‑1 (MQ=255)
   atTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGCTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAg                                                                                                                                                 >  2:1254982/1‑150 (MQ=255)
    tttCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTCTCAGATATTCTACTGTTGACTATTTGTGCCGTTAATTCTGGTGCAGAAGGTTGGGAACATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGc                                                                                                                                                <  2:13601/150‑1 (MQ=255)
         aTTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAatat                                                                                                                                           >  1:557372/1‑150 (MQ=255)
            acacCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGt                                                                                                                                        <  2:845783/150‑1 (MQ=255)
             cacCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGAAATACTACTATTGACTATTTTTGACATATTATCAGCACAAGAATGATCGCAAGAAAAAGAGGATTTCGGGTAAACAACTCACGATTTATTTAATCAATATGGTg                                                                                                                                       >  1:997696/1‑150 (MQ=9)
              acCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGa                                                                                                                                      <  2:23751/150‑1 (MQ=255)
               cccGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGAt                                                                                                                                     >  1:429415/1‑150 (MQ=255)
               cccGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTATTGAAGCAATATGGTGAt                                                                                                                                     >  1:170655/1‑150 (MQ=255)
               cccGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGAGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGAt                                                                                                                                     <  2:655582/150‑1 (MQ=255)
                  gATTACAGACAAGCCTGGAATGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGAtttt                                                                                                                                  <  2:1016297/150‑1 (MQ=255)
                   aTTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTg                                                                                                                                 <  2:549036/150‑1 (MQ=255)
                     tagaGACAAGCCTTGAAAGTGGAGCATAAATTGTCAGTTATTCTACTTTTGATTATTTGTGCCGTTATTTCTGGTGAAGAAAGTTGGTAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTCTTGAAGCAATATGGTGATTTTGaa                                                                                                                               <  1:208710/147‑1 (MQ=11)
                        aGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTTTGCCGTTATTTCTGGTGCAGAATGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAAt                                                                                                                            <  1:164322/150‑1 (MQ=255)
                         gACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATTGTGATTTTGAAAATg                                                                                                                           >  2:715967/1‑150 (MQ=255)
                         caCAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATg                                                                                                                           <  2:578766/149‑1 (MQ=255)
                           aaaGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTTTTGACTATTTTTGCCGTTATTTCTGGTGCAGAAAGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTCTTGAAGCAATATGGTGATTTTGAAAATGGt                                                                                                                         <  2:937077/149‑1 (MQ=255)
                             aGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAt                                                                                                                       >  1:530351/1‑150 (MQ=255)
                              gCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAtt                                                                                                                      >  2:540410/1‑150 (MQ=255)
                              gCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTAtt                                                                                                                      >  1:489024/1‑150 (MQ=255)
                               ccTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTc                                                                                                                     >  2:1121431/1‑150 (MQ=255)
                                cTGGACAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTTTTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTcc                                                                                                                    <  1:1232372/150‑1 (MQ=255)
                                cTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAAGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTCTTGAAGc                                                                                                                                                >  2:960519/1‑122 (MQ=255)
                                cTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTCTTGACTATTTGTGCCGTTATTTCTGGTGCAGAATGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTCTTGAAGc                                                                                                                                                <  1:960519/122‑1 (MQ=255)
                                 tGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCt                                                                                                                   <  2:699672/150‑1 (MQ=255)
                                 tGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCt                                                                                                                   >  1:1042847/1‑150 (MQ=255)
                                         ggtgCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATATCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGa                                                                                                           >  1:1145108/1‑150 (MQ=255)
                                           tgCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATa                                                                                                         <  2:668326/150‑1 (MQ=38)
                                           tgCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTCTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATa                                                                                                         >  2:41972/1‑150 (MQ=38)
                                                 aaTTGTCAGATATTCTACTGTTGACAATTTGTGCCGTAATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTg                                                                                                   <  1:541657/150‑1 (MQ=255)
                                                   ttGTCAGAAATTCTTATGATCACTATTTGTGCCGTTATTTCTGGTGCAGACCGTTGGGACGTTAAAGAGGATTTCGGGGATACTCATCTCGATTTTTTGAAGCAATATGGTTATTTTGAAAATTGTATTCCTGTTCACGATACCATTGcc                                                                                                 <  2:1176611/150‑1 (MQ=11)
                                                            atTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTa                                                                                        <  1:482981/150‑1 (MQ=32)
                                                                 aCTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTGTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACTATTGCCAGAGTTGTATcctg                                                                                   >  1:820375/1‑148 (MQ=2)
                                                                  cTGTTGACTATTTGTGCGGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAg                                                                                  <  2:666265/150‑1 (MQ=21)
                                                                    gTTGACAATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGTTATAGACTATTTCGTGTAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCCGAGTTGTCTCCCAggg                                                                                <  2:1228315/150‑1 (MQ=11)
                                                                      tGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGaa                                                                              >  2:1108963/1‑150 (MQ=38)
                                                                      tGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATACCTGTTCACGATACCATTGCCAGAGTTGTATCCAAGGGaa                                                                              >  1:826705/1‑150 (MQ=25)
                                                                            tttGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCa                                                                        <  1:792687/150‑1 (MQ=255)
                                                                               gtgCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGt                                                                     <  1:715005/150‑1 (MQ=255)
                                                                                    gTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGc                                                                <  2:435367/150‑1 (MQ=255)
                                                                                          tCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACttt                                                          >  1:613673/1‑150 (MQ=255)
                                                                                            gggTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACttttt                                                        <  2:539550/149‑1 (MQ=255)
                                                                                               tGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGt                                                     <  2:648861/150‑1 (MQ=255)
                                                                                               ggCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGt                                                     <  1:715967/149‑1 (MQ=255)
                                                                                                gCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTACACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTa                                                    <  1:1185508/150‑1 (MQ=255)
                                                                                                    aaGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTTTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTaa                                                >  2:675246/1‑150 (MQ=255)
                                                                                                      ggTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAaat                                              >  2:838215/1‑150 (MQ=255)
                                                                                                         tGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAaataat                                           >  1:354912/1‑150 (MQ=255)
                                                                                                           ggAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTc                                         >  1:862342/1‑150 (MQ=255)
                                                                                                                 tataGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCAttt                                   >  2:161083/1‑150 (MQ=255)
                                                                                                                 tataGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCAttt                                   >  2:699492/1‑150 (MQ=255)
                                                                                                                   taGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAGAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTa                                 >  2:261152/1‑150 (MQ=255)
                                                                                                                     gagGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATg                               >  1:327999/1‑150 (MQ=255)
                                                                                                                       ggATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGtt                             >  2:1288893/1‑150 (MQ=255)
                                                                                                                        gATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGttt                            >  2:1231219/1‑150 (MQ=255)
                                                                                                                              ggtgAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAAt                      <  2:860450/147‑1 (MQ=255)
                                                                                                                              ggggAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAAt                      <  1:1254982/150‑1 (MQ=255)
                                                                                                                              ggggAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAAt                      >  1:151744/1‑150 (MQ=255)
                                                                                                                                 gAAACACATCTCGATTTTTTGAAGCAATATGGTTATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTga                   <  2:182306/150‑1 (MQ=255)
                                                                                                                                 gAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTGCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTga                   <  2:540297/150‑1 (MQ=255)
                                                                                                                                  aaaCACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTgag                  <  2:429415/150‑1 (MQ=255)
                                                                                                                                       caTCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATAtt             >  1:389585/1‑150 (MQ=255)
                                                                                                                                       caTCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATAtt             >  2:8726/1‑150 (MQ=255)
                                                                                                                                       caTCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCAAGTAACTCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATAtt             >  1:144954/1‑150 (MQ=255)
                                                                                                                                         tctcGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTca           >  1:313896/1‑149 (MQ=255)
                                                                                                                                           acGATTATTTGAAGCAATATGGTGTTGTTGAAACTGGTATTCCTGTTCACGAAACCATTGCCAGAGTTGTATCCCATGGAAAGATCACGTAACGCTACTATTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCtt         <  2:176958/149‑1 (MQ=37)
                                                                                                                                            cGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACTCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCtta        >  1:262688/1‑150 (MQ=255)
                                                                                                                                               ttttttGAAGCATTATGGTGATTTTGAAAATGGTATTCCTGTTCACGATAGCATTGCCAGTTTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCttatta     <  1:284128/150‑1 (MQ=255)
                                                                                                                                               ttttttGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATAACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTTTGTTTAAAAATTGAGATATTCCttatta     <  1:551566/150‑1 (MQ=255)
                                                                                                                                                  tttGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCt  >  1:1040599/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CATATTTCTATTACACCCGATTACAGACAAGCCTGGAAAGTGGTGCATAAATTGTCAGATATTCTACTGTTGACTATTTGTGCCGTTATTTCTGGTGCAGAAGGTTGGGAAGATATAGAGGATTTCGGGGAAACACATCTCGATTTTTTGAAGCAATATGGTGATTTTGAAAATGGTATTCCTGTTCACGATACCATTGCCAGAGTTGTATCCCAGGGAAAGATCACGTAACGCTACTTTTTTGTACTAAATAATTCGCATTTTATGTTTAAAAATTGAGATATTCCTTATTACCT  >  NC_000913/736469‑736764

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