Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,365 A→C 8.0% intergenic (+129/+103) rhlE → / ← ybiA ATP‑dependent RNA helicase RhlE/N‑glycosidase YbiA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,3650AC 8.0% 346.0 / 17.4 111intergenic (+129/+103)rhlE/ybiAATP‑dependent RNA helicase RhlE/N‑glycosidase YbiA
Reads supporting (aligned to +/- strand):  ref base A (54/48);  new base C (4/5);  total (58/53)
Fisher's exact test for biased strand distribution p-value = 7.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

AAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAG  >  NC_000913/832216‑832513
                                                                                                                                                     |                                                                                                                                                    
aaaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATaa                                                                                                                                                      <  1:128353/150‑1 (MQ=255)
aaaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATaa                                                                                                                                                      >  2:750686/1‑150 (MQ=255)
 aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAg                                                                                                                                                     <  2:616257/150‑1 (MQ=255)
 aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAg                                                                                                                                                     >  1:1315547/1‑150 (MQ=255)
  aCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAgg                                                                                                                                                    <  2:1070901/150‑1 (MQ=255)
    cTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCg                                                                                                                                                  <  2:321612/150‑1 (MQ=255)
         gCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAc                                                                                                                                             <  1:1313155/150‑1 (MQ=255)
          cTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcg                                                                                                                                            >  1:426895/1‑150 (MQ=255)
             cgcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                         <  1:1003523/150‑1 (MQ=255)
             cgcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                         <  1:720745/150‑1 (MQ=255)
               cAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                       <  1:947355/150‑1 (MQ=255)
                 gTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                     <  2:386192/150‑1 (MQ=255)
                  tAATCTTTTATGCCGTGCTATGCCCGGCATCAGGGTGATGAACAACCGCAATACTGCCTGATGCGCTTCGCTTATCAGGCCTACGTGAACTCTTCAATTTATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    <  1:1479899/150‑1 (MQ=21)
                  tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                    <  2:1475415/150‑1 (MQ=25)
                   aaTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                                                                                                                                   >  2:1038740/1‑150 (MQ=255)
                    aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                  <  2:1454293/150‑1 (MQ=255)
                    aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCgg                                                                                                                                  <  2:1435945/150‑1 (MQ=255)
                      gcctgaTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                >  2:431997/6‑150 (MQ=255)
                      cTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcc                                                                                                                                <  2:1232003/150‑2 (MQ=21)
                      cTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCa                                                                                                                                <  1:289014/150‑1 (MQ=255)
                      cTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGAT‑‑GCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                              >  2:1457764/1‑150 (MQ=255)
                       cctgaTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                               >  1:1522750/5‑150 (MQ=255)
                       cctgaTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                               >  1:76583/5‑150 (MQ=255)
                       cctgaTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAt                                                                                                                               >  1:1182116/5‑150 (MQ=255)
                        tttATGCCGGGCTATGCCCGGCATCAGGCTGATTAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAtt                                                                                                                              <  2:670765/150‑1 (MQ=255)
                         ttATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCCGATGCGCTACGCTTATCAGGCCTACGTTAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctac                                                                                                                             <  1:158727/150‑5 (MQ=21)
                          tATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATGCGGCGTAAACGCCTTATCCGGcctaca                                                                                                                            <  2:176246/150‑6 (MQ=21)
                           aTGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTca                                                                                                                           <  2:1138702/150‑1 (MQ=255)
                            tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTcac                                                                                                                          <  1:1064142/150‑1 (MQ=255)
                                 ggCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCAcaacaa                                                                                                                     <  2:528512/150‑1 (MQ=255)
                                      tGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACt                                                                                                                >  1:850337/1‑150 (MQ=255)
                                       gCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACtt                                                                                                               >  2:126840/1‑150 (MQ=255)
                                         ccGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGt                                                                                                             >  1:491288/1‑150 (MQ=255)
                                            gCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAg                                                                                                          <  2:1505773/150‑1 (MQ=255)
                                            gCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCACTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAg                                                                                                          <  1:220250/150‑1 (MQ=255)
                                             cATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGc                                                                                                         >  2:202841/1‑150 (MQ=255)
                                                cAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAAt                                                                                                      >  1:765196/1‑150 (MQ=255)
                                                 aGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAAtt                                                                                                     >  1:1161464/1‑150 (MQ=255)
                                                   gCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTg                                                                                                   <  2:902786/150‑1 (MQ=255)
                                                      gatgaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTAAGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTgaga                                                                                                >  2:746303/1‑150 (MQ=255)
                                                         gaACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAaca                                                                                             >  2:447659/1‑150 (MQ=255)
                                                          aaCAAACGCAAACCTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAacac                                                                                            >  1:1334102/1‑150 (MQ=255)
                                                            cacaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACgg                                                                                          <  2:604075/147‑1 (MQ=255)
                                                            cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTc                                                                                                            <  1:1041249/132‑1 (MQ=255)
                                                            cAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTc                                                                                                            >  2:1041249/1‑132 (MQ=255)
                                                              aaCGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGa                                                                                        >  1:1374540/1‑150 (MQ=255)
                                                                cGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGaaa                                                                                      >  2:985606/1‑150 (MQ=255)
                                                                  cAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATTTTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGaaaaa                                                                                    >  2:1429505/1‑150 (MQ=255)
                                                                       cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTGAGCAATTTGAGAACACGGGAAAAAATTTa                                                                               <  2:633605/150‑1 (MQ=255)
                                                                       cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTa                                                                               <  1:840225/150‑1 (MQ=255)
                                                                       cTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTa                                                                               <  2:1405130/150‑1 (MQ=255)
                                                                         gCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATc                                                                             <  2:1310957/150‑1 (MQ=255)
                                                                         gCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATc                                                                             <  2:15756/150‑1 (MQ=255)
                                                                         gCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATc                                                                             <  1:187971/150‑1 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTTTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTg                                                                           >  2:403558/1‑150 (MQ=255)
                                                                           cTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTg                                                                           >  2:893091/1‑150 (MQ=255)
                                                                              aTGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGttt                                                                        >  1:1019016/1‑150 (MQ=255)
                                                                              aTGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGttt                                                                        <  1:837477/150‑1 (MQ=255)
                                                                                gcgcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTa                                                                      <  1:397876/150‑1 (MQ=255)
                                                                                 cgcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAc                                                                     <  2:377640/150‑1 (MQ=255)
                                                                                 cgcTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAc                                                                     <  2:170860/150‑1 (MQ=255)
                                                                                      cGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGccc                                                                <  1:746303/150‑1 (MQ=255)
                                                                                           aTCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCCCTTGTCCGCAATTTGAGAACACGGGAAAAAATTTATCAGTTTTACCGCCCATAAg                                                           >  2:1184167/1‑150 (MQ=255)
                                                                                            tCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGc                                                          >  2:920022/1‑150 (MQ=255)
                                                                                              aGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCa                                                        >  2:506815/1‑150 (MQ=255)
                                                                                                    aCGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCgaga                                                  >  1:1166976/1‑150 (MQ=255)
                                                                                                       tGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGt                                               >  1:418333/1‑150 (MQ=255)
                                                                                                         aaCTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTa                                             >  1:282436/1‑150 (MQ=255)
                                                                                                               tcAATTTATTGAATTTGATTGCTTTTGTAGGCCTGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTTTTTTACGGCCCATAAGCCACCGAGACGTTACTtctc                                       <  2:1315547/149‑1 (MQ=38)
                                                                                                                cAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTctct                                      <  2:626478/150‑1 (MQ=255)
                                                                                                                 aaTATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCta                                     <  1:1038740/150‑1 (MQ=255)
                                                                                                                    atatTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAg                                  <  1:191180/150‑1 (MQ=255)
                                                                                                                    atatTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAg                                  <  2:1019016/150‑1 (MQ=255)
                                                                                                                     tatTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGc                                 >  2:727527/1‑150 (MQ=255)
                                                                                                                     aatTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGc                                 <  2:1295020/149‑1 (MQ=255)
                                                                                                                      atTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGcc                                >  2:1327289/1‑150 (MQ=255)
                                                                                                                           aTTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTg                           <  2:426895/150‑1 (MQ=255)
                                                                                                                                cATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTcgc                      <  1:920022/150‑1 (MQ=255)
                                                                                                                                  tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCgcgc                    >  2:122048/1‑150 (MQ=255)
                                                                                                                                  tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCgcgc                    <  2:850337/150‑1 (MQ=255)
                                                                                                                                  tGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCgcgc                    >  2:1336365/1‑150 (MQ=255)
                                                                                                                                    cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCaa                  >  1:45076/1‑150 (MQ=255)
                                                                                                                                    cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCaa                  >  2:206883/1‑150 (MQ=255)
                                                                                                                                    cTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTAACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCaa                  >  1:873724/1‑150 (MQ=255)
                                                                                                                                      tttGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACTATCACTTGTCAGCAATTTGAGAACACGTGAAAAAATTTATGTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACt                <  2:807123/150‑1 (MQ=255)
                                                                                                                                      gtgaaatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGAGAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACt                >  1:480802/8‑150 (MQ=38)
                                                                                                                                        tGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTcc              >  1:1514218/1‑150 (MQ=255)
                                                                                                                                         gTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCa             >  2:1036447/1‑150 (MQ=255)
                                                                                                                                          aatgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAt            >  1:1073556/4‑150 (MQ=255)
                                                                                                                                           atgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAtt           >  1:1164259/3‑150 (MQ=255)
                                                                                                                                           aGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCAtt           >  1:1056464/1‑150 (MQ=255)
                                                                                                                                            tgCCGGATGCGGCGTAAACGCCTTATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa          >  1:490922/2‑150 (MQ=255)
                                                                                                                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCGGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGATGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTa          >  2:577380/1‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTaaa        >  1:1127418/1‑150 (MQ=255)
                                                                                                                                              ccGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTaaa        >  2:1287804/1‑150 (MQ=255)
                                                                                                                                                ggATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAg      <  2:884699/150‑1 (MQ=255)
                                                                                                                                                 gATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAgg     >  1:1303749/1‑150 (MQ=255)
                                                                                                                                                   tAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTa   >  2:1065124/1‑150 (MQ=255)
                                                                                                                                                    aaGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAg  <  1:1429505/150‑1 (MQ=255)
                                                                                                                                                    aaGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAg  <  1:202841/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAG  >  NC_000913/832216‑832513

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