Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,488,051 G→A 6.2% intergenic (+16/‑181) ydcF → / → aldA DUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,488,0510GA 6.2% 406.9 / 9.5 129intergenic (+16/‑181)ydcF/aldADUF218 superfamily protein, SAM‑binding/aldehyde dehydrogenase A, NAD‑linked
Reads supporting (aligned to +/- strand):  ref base G (54/66);  new base A (3/5);  total (57/72)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.95e-01

Only 100 of 134 total aligned reads displayed.

GCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCG  >  NC_000913/1487902‑1488200
                                                                                                                                                     |                                                                                                                                                     
gCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTtgtg                                                                                                                                                       >  1:118298/1‑150 (MQ=255)
gCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTtgtg                                                                                                                                                       <  1:1370499/150‑1 (MQ=255)
 cGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGa                                                                                                                                                      >  2:1574436/1‑150 (MQ=255)
  gATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGaa                                                                                                                                                     <  2:1662414/150‑1 (MQ=255)
  gATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGaa                                                                                                                                                     <  1:1003726/150‑1 (MQ=255)
   aTGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAAc                                                                                                                                                    >  1:1719058/1‑150 (MQ=255)
   aTGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATACATGCCTGGCAAACGCTGACACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAAc                                                                                                                                                    <  1:870870/150‑1 (MQ=255)
     ggCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCa                                                                                                                                                  >  2:1280262/1‑150 (MQ=255)
       cTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCGCGTTACGTTAAAAATTGCCCGTTTGTGAACCACt                                                                                                                                                <  2:208884/150‑1 (MQ=255)
           ggTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                            >  1:1217518/1‑150 (MQ=255)
           ggTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGtt                                                                                                                                            <  2:1194361/150‑1 (MQ=255)
            gTCCCCGCGGTCGAGATTTTATCGTTCACGTTGAGTTTCCGGCAGAAGGCATCCATGAATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGttt                                                                                                                                           <  1:266600/150‑1 (MQ=255)
               cccGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCa                                                                                                                                        <  1:168688/150‑1 (MQ=255)
                ccGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCaa                                                                                                                                       <  2:1624497/150‑1 (MQ=255)
                   cgGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACg                                                                                                                                    <  1:1410787/150‑1 (MQ=255)
                     ggCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACggg                                                                                                                                  >  1:1452779/1‑150 (MQ=255)
                      gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCTTCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGc                                                                                                                                 <  2:1164471/150‑1 (MQ=255)
                      gCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGc                                                                                                                                 <  2:1044616/150‑1 (MQ=255)
                       cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                >  1:866806/1‑150 (MQ=255)
                       cGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCAGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCa                                                                                                                                >  1:612843/1‑150 (MQ=255)
                        gagaTTTTATCGTTCACGTTGATTTTCCGCCAGAAGTCATCCATGCATGGCAACCGCTGAAACATGATGCGGAGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCAt                                                                                                                               <  2:387647/150‑1 (MQ=255)
                         agaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATg                                                                                                                              >  2:311433/1‑150 (MQ=255)
                          gaTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCaatt                                                                                                                             <  1:1290969/150‑4 (MQ=38)
                           aTTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGAGGCGGTGCTCATCGAGGCGCTGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGAc                                                                                                                            <  1:327038/150‑1 (MQ=255)
                            ttttATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTACAAACGGGCATGACt                                                                                                                           >  1:114912/1‑150 (MQ=255)
                                    ttCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACtt                                                                                                                   <  1:157341/150‑1 (MQ=255)
                                      cACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACtttt                                                                                                                 >  2:689271/1‑150 (MQ=255)
                                      cACGTTGATTTTCCGGCAGAAGTCATCCATGCATCGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACtttt                                                                                                                 <  1:478736/150‑1 (MQ=255)
                                       aCGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTa                                                                                                                <  2:280399/150‑1 (MQ=255)
                                       aCGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTa                                                                                                                >  2:382107/1‑150 (MQ=255)
                                       aCGTTGAGTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAAGATGATGCGGAGCTCATCGAGGCGATTGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTa                                                                                                                <  1:305353/150‑1 (MQ=255)
                                        cGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAt                                                                                                               >  2:1276631/1‑150 (MQ=255)
                                        cGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAt                                                                                                               >  1:1505852/1‑150 (MQ=255)
                                        cGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTAt                                                                                                               >  2:829688/1‑150 (MQ=255)
                                          ttGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTTCGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTCCAAACGGGCATGACTCCTGACTTTTAttt                                                                                                             >  2:1081346/1‑150 (MQ=255)
                                             aTTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTg                                                                                                          <  2:425731/150‑1 (MQ=255)
                                              ttttCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGTAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGc                                                                                                         <  2:1310416/150‑1 (MQ=255)
                                                   cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                    >  1:1017549/1‑150 (MQ=255)
                                                   cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                    >  2:27325/1‑150 (MQ=255)
                                                   cGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCtttt                                                                                                    >  2:317071/1‑150 (MQ=255)
                                                    ggCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTa                                                                                                   <  2:1183325/150‑1 (MQ=255)
                                                       aGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTc                                                                                                <  1:311433/150‑1 (MQ=255)
                                                        taaGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTcc                                                                                               <  2:1500746/149‑1 (MQ=255)
                                                        gAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTcc                                                                                               <  2:389189/150‑1 (MQ=255)
                                                         aaGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCt                                                                                              <  2:787306/150‑1 (MQ=255)
                                                          aGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtt                                                                                             >  2:177105/1‑150 (MQ=255)
                                                          aGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTAATCCtt                                                                                             >  2:1144592/1‑150 (MQ=255)
                                                           gTCATCCATGCATGGCAAACGCTGAAACATGCTGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCttt                                                                                            <  2:1159084/150‑1 (MQ=255)
                                                           gTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCttt                                                                                            <  2:866806/150‑1 (MQ=255)
                                                            tCATCCATGCATGTCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                           <  2:412401/150‑1 (MQ=255)
                                                            gcATCCATGCATGGCAACCGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAATATTGCCCGTTTGTGAACCAATTGTTTGCAAAAGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCtttt                                                                                           <  2:277217/149‑1 (MQ=255)
                                                             cATCCATGCATGGCACAAGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGACCCTGTTGTGTGCAAACGGGGAGTACTCATGACTTTTATTTCAGCATTTTATTCCTGTTa                                                                                          >  2:111084/1‑150 (MQ=255)
                                                             cATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTa                                                                                          >  2:1426397/1‑150 (MQ=255)
                                                             cATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTAATTCCTTTTa                                                                                          >  1:69469/1‑150 (MQ=255)
                                                                  atgcatgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCATGACTTTTATTTCTGCCTTTTATTCCTTTTACACtt                                                                                     >  1:280521/1‑150 (MQ=255)
                                                                   tgcatgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTg                                                                                    <  1:12579/150‑1 (MQ=255)
                                                                   tgcatgGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTg                                                                                    <  1:1710320/150‑1 (MQ=255)
                                                                          cAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTAt                                                                             >  1:1413594/1‑150 (MQ=255)
                                                                          cAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGACATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTAt                                                                             <  2:998712/150‑1 (MQ=255)
                                                                            accGCTGAAACATGATGCGGTGCTCATCTAGGCGATGGAAAGTCCCTCGTTAAGTTAAAAATTGCCCGTTTGTGAACCACTTGGTTGAAAACTGTCATGACTCATTCCATTTATTTCTGGCTTTTATTCATTTTAGACTTGTTTTTAtaa                                                                           >  1:142107/3‑148 (MQ=255)
                                                                                cTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGcc                                                                       >  2:891174/1‑150 (MQ=255)
                                                                                    aaCATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGATACGTTAAAAATTGCCCGTTTGTGAACCACTTTTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTc                                                                   >  2:1161768/1‑150 (MQ=255)
                                                                                     aCATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTca                                                                  <  2:87750/150‑1 (MQ=255)
                                                                                         tattcgttGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCATTTGTTTGCAAACGGGCATGACTCCTGTCTTTTATTTCTGCCTTTTATTCCTTTTTCACTTGTTTTTATGAAGCCATTAACAGa                                                              <  2:677094/143‑1 (MQ=255)
                                                                                                tGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTcc                                                       <  2:1017549/150‑1 (MQ=255)
                                                                                                tGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTcc                                                       <  2:983139/150‑1 (MQ=255)
                                                                                                   tCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCttt                                                    >  2:936/1‑150 (MQ=255)
                                                                                                    cATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGATTTGTCCTTTc                                                   >  1:282198/1‑150 (MQ=255)
                                                                                                     aTCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCa                                                  <  2:1206423/150‑1 (MQ=255)
                                                                                                     aTCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCa                                                  <  2:1719058/150‑1 (MQ=255)
                                                                                                      tCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCAc                                                 <  1:1082914/150‑1 (MQ=255)
                                                                                                       cGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACg                                                >  2:1648203/1‑150 (MQ=255)
                                                                                                       cGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTAAACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACg                                                >  2:1146742/1‑150 (MQ=255)
                                                                                                        gAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGCCTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGa                                               >  1:1224354/1‑150 (MQ=255)
                                                                                                             gATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCg                                          <  1:546104/150‑1 (MQ=255)
                                                                                                               tGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                        <  1:664156/150‑1 (MQ=255)
                                                                                                               gggAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtc                                        <  1:1103041/149‑1 (MQ=255)
                                                                                                                 gAAAGTCTCTCGTTACGTTAAAAATTGCCCGTTTGTGATCCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                      <  2:52717/150‑1 (MQ=255)
                                                                                                                 gAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGtctc                                      <  1:526408/150‑1 (MQ=255)
                                                                                                                     gTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtg                                  >  1:1072297/1‑150 (MQ=255)
                                                                                                                     gTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtg                                  >  1:1037258/1‑150 (MQ=255)
                                                                                                                      tCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTGGGCAAACGTGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTAGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtga                                 <  2:85831/150‑1 (MQ=255)
                                                                                                                       cGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgat                                >  2:572808/1‑150 (MQ=255)
                                                                                                                        gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                               >  2:718083/1‑150 (MQ=255)
                                                                                                                        gCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCtgatg                               <  1:1276631/150‑1 (MQ=255)
                                                                                                                             ttacgttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGAtgattga                          >  2:1401851/1‑150 (MQ=255)
                                                                                                                                 gttaAAAATTTCCCGTTTGTTAACAACTTGTTTGCAAACGGGCATGACTCCTGAATTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtt                      <  2:550251/150‑1 (MQ=255)
                                                                                                                                 gttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtt                      <  2:1413594/150‑1 (MQ=255)
                                                                                                                                 gttaAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtt                      >  2:408191/1‑150 (MQ=255)
                                                                                                                                 gttaAAAATTGCCAGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGtt                      >  1:738975/1‑150 (MQ=255)
                                                                                                                                      gtcatGCCCGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGttaatta                 >  1:778606/5‑150 (MQ=38)
                                                                                                                                        aTTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATTTTAATTAAc               <  1:605386/150‑1 (MQ=255)
                                                                                                                                          tGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGACGATTGATGTTAATTAACaa             <  1:942013/150‑1 (MQ=255)
                                                                                                                                            cccGTTTGCAAACAAGTGGTTCACAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATg           >  2:706138/1‑150 (MQ=255)
                                                                                                                                              cGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTGCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGTTTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTa         <  2:1592959/150‑1 (MQ=255)
                                                                                                                                                tttGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTAtt       <  2:1037258/150‑1 (MQ=255)
                                                                                                                                                tttGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTAtt       <  2:1072297/150‑1 (MQ=255)
                                                                                                                                                   gtgAACCACTTGTTTTCAAACGGGAATGACTCATTACTTTTATTTCTGACTTTTATTCCTTTTACACTTGTTTTTATTAAGCCCTTCTCAGAATTGTCCTTTCACGATTCCGTCTCTCTGATTATTGATGTTAATTAACAATGTATTCAc    <  2:935537/150‑1 (MQ=255)
                                                                                                                                                     gAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTTATGAACAATGTATTCACCg  >  2:346060/1‑150 (MQ=255)
                                                                                                                                                     gAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCg  >  1:1149577/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGATGGCTACGGTCCCCGCGGGCGAGATTTTATCGTTCACGTTGATTTTCCGGCAGAAGTCATCCATGCATGGCAAACGCTGAAACATGATGCGGTGCTCATCGAGGCGATGGAAAGTCGCTCGTTACGTTAAAAATTGCCCGTTTGTGAACCACTTGTTTGCAAACGGGCATGACTCCTGACTTTTATTTCTGCCTTTTATTCCTTTTACACTTGTTTTTATGAAGCCCTTCACAGAATTGTCCTTTCACGATTCCGTCTCTCTGATGATTGATGTTAATTAACAATGTATTCACCG  >  NC_000913/1487902‑1488200

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: