Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,530,902 A→G V513A (GTC→GCC)  yhgE ← DUF4153 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,530,9020AG98.7% 294.6 / ‑4.1 76V513A (GTC→GCC) yhgEDUF4153 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (39/35);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

GATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGACTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGCATCA  >  NC_000913/3530759‑3531048
                                                                                                                                                |                                                                                                                                                  
gATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGtc                                                                                                                                             >  2:686489/1‑151 (MQ=255)
   cgcgCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGtcatc                                                                                                                                          >  2:794754/1‑151 (MQ=255)
         aCGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTaa                                                                                                                                     >  2:944772/1‑150 (MQ=255)
         aCGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTaa                                                                                                                                     >  1:672871/1‑150 (MQ=255)
          cGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTaaa                                                                                                                                    >  1:481366/1‑150 (MQ=255)
               tAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCaa                                                                                                                               >  1:30506/1‑150 (MQ=255)
                aGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAAc                                                                                                                             <  2:716895/151‑1 (MQ=255)
                  gCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAAc                                                                                                                             >  2:833472/1‑149 (MQ=255)
                    ttCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAgg                                                                                                                          >  1:46204/1‑150 (MQ=255)
                     tCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGa                                                                                                                         <  2:834972/150‑1 (MQ=255)
                     tCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGa                                                                                                                         <  2:382734/150‑1 (MQ=255)
                     tCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAATAAACTCATATCAAGCGCGTCCCATTCTTTACTGTCT‑GAGTCAAAGCCATAGACAATGGCACTGTCATCATTAAAAGCAAACCGGa                                                                                                                         <  2:815412/150‑1 (MQ=255)
                      cGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtc                                                                                                                       <  1:184905/151‑1 (MQ=255)
                      cGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAt                                                                                                                        <  2:232474/150‑1 (MQ=255)
                         ccccAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtccg                                                                                                                     <  1:34726/150‑1 (MQ=255)
                           ccAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtccgctc                                                                                                                  <  1:480747/151‑1 (MQ=255)
                            cAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtccgctcc                                                                                                                 <  2:620833/151‑1 (MQ=255)
                             aGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtccgctccg                                                                                                                <  2:442454/151‑1 (MQ=255)
                              gTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtccgctccgc                                                                                                               >  1:976151/1‑151 (MQ=255)
                              gTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGAtccgctccgc                                                                                                               >  1:835847/1‑151 (MQ=255)
                                   ccATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGAc                                                                                                          <  1:463194/151‑1 (MQ=255)
                                       ccTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAGTGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATc                                                                                                      >  2:521164/1‑151 (MQ=255)
                                       ccTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATc                                                                                                      <  1:253303/151‑1 (MQ=255)
                                        cTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATc                                                                                                      >  2:186942/1‑150 (MQ=255)
                                          tCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCAc                                                                                                    <  1:388132/150‑1 (MQ=255)
                                                    gAGTAATTTGTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCG‑GAGTCAAATCCTTAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATc                                                                                          <  2:702290/150‑1 (MQ=255)
                                                     aGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCtt                                                                                        <  2:888099/151‑1 (MQ=255)
                                                     aGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCtt                                                                                        <  1:293571/151‑1 (MQ=255)
                                                        aaTTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCt                                                                                     >  2:902089/1‑151 (MQ=255)
                                                           ttttCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCg                                                                                   >  2:521973/1‑150 (MQ=255)
                                                              tctttcGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGAcc                                                                                <  2:30502/150‑1 (MQ=255)
                                                              tctttcGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGAcc                                                                                <  2:91728/150‑1 (MQ=255)
                                                              tctttcGTTATTTCGTTCGGAAGAAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAACGCCATAGACAATGGCCCTGTCATCATTAAAAGCAAACAGGATCCGCCCCTCCTGACCATCCCTATTCAGATCTTGCTCGAcc                                                                                >  2:1072583/1‑150 (MQ=255)
                                                                 ttcGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCa                                                                               <  1:1018386/148‑1 (MQ=255)
                                                                 ttcGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAgg                                                                             >  2:736905/1‑150 (MQ=255)
                                                                   cGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACg                                                                          >  1:92623/1‑151 (MQ=255)
                                                                    gTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGT‑CCATTCTTTTTTTGCTGGGGCCAAAGCCATAAACAATGCCCTTGTCTCCTTTAAAACCAAAAAGAATCCCCTCCGCCTAACAATACCTATTCAAACTTTCCCCGACCAg                                                                              >  2:224680/1‑146 (MQ=255)
                                                                        tttCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGc                                                                     >  1:213025/1‑151 (MQ=255)
                                                                        tttCGTTCGGAAGTAAACTCATATCAAGAGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGc                                                                     >  2:493498/1‑151 (MQ=255)
                                                                         ttcgttcgGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGc                                                                     <  2:264797/150‑1 (MQ=255)
                                                                         ttcgttcgGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCg                                                                    >  1:703566/1‑151 (MQ=255)
                                                                            gttcgGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCt                                                                 >  2:95193/1‑151 (MQ=255)
                                                                             ttcgGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCtt                                                                >  1:1026819/1‑151 (MQ=255)
                                                                               cgGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCttt                                                               <  2:431690/150‑1 (MQ=255)
                                                                                ggAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCttttt                                                             <  2:337412/151‑1 (MQ=255)
                                                                                  aaGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTc                                                            >  2:559705/1‑150 (MQ=255)
                                                                                    gTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAAt                                                         >  1:785301/1‑151 (MQ=255)
                                                                                    gTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAAt                                                         >  1:227588/1‑151 (MQ=255)
                                                                                     tAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAAt                                                         >  2:722852/1‑150 (MQ=255)
                                                                                       aaCTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGACCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACa                                                      <  1:1071654/151‑1 (MQ=255)
                                                                                        aCTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAg                                                     <  2:589287/151‑1 (MQ=255)
                                                                                            atatCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGtc                                                  >  2:292603/1‑150 (MQ=255)
                                                                                                     gcgTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTtat                                        <  1:316874/151‑1 (MQ=255)
                                                                                                      cgTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTtata                                       <  2:921185/151‑1 (MQ=255)
                                                                                                        tCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGc                                     >  1:274513/1‑151 (MQ=255)
                                                                                                            aTTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATc                                 <  2:692980/151‑1 (MQ=255)
                                                                                                            aTTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATc                                 >  1:537157/1‑151 (MQ=255)
                                                                                                                ttttCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTgg                             >  2:524585/1‑151 (MQ=255)
                                                                                                                 tttCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCt                               <  1:1050962/148‑1 (MQ=255)
                                                                                                                 tttCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTgg                             >  2:411783/1‑150 (MQ=255)
                                                                                                                 tttCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTg                              <  2:981671/149‑1 (MQ=255)
                                                                                                                  ttCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGa                            <  2:795896/150‑1 (MQ=255)
                                                                                                                  ttCTGTCT‑GAGTCAAAGCCATAGACAATGGCTCTCTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGAt                           >  1:822567/1‑151 (MQ=255)
                                                                                                                    ctgtct‑gAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGAtt                          >  2:532526/1‑150 (MQ=255)
                                                                                                                      gtct‑gAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTaa                        <  1:1124746/150‑1 (MQ=255)
                                                                                                                        ct‑gAGTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGc                     >  2:572574/1‑151 (MQ=255)
                                                                                                                             gTCAAAGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGAc                 >  1:153868/1‑151 (MQ=255)
                                                                                                                                aaaGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGAcc                >  2:620862/1‑149 (MQ=255)
                                                                                                                                aaaGCCATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGAcc                >  1:211981/1‑149 (MQ=255)
                                                                                                                                    ccATAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATg          >  2:961800/1‑151 (MQ=255)
                                                                                                                                      aTAGACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTc        <  1:1048883/151‑1 (MQ=255)
                                                                                                                                         gACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGCa     <  2:703390/151‑1 (MQ=255)
                                                                                                                                         gACAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGCa     <  1:965161/151‑1 (MQ=255)
                                                                                                                                         gACAATGACTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGc      >  2:195114/1‑150 (MQ=255)
                                                                                                                                           cAATGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGCATc   >  2:681099/1‑151 (MQ=255)
                                                                                                                                            aaTGGCTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGCATCa  <  2:306220/151‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
GATCGCGCCACGCTTTAGGCTTCGTCCCCAGTTTTCCATCCTTCGCAGCTGTGAGTAATTTTTCTTTCGTTATTTCGTTCGGAAGTAAACTCATATCAAGCGCGTCCCATTCTTTTCTGTCT‑GAGTCAAAGCCATAGACAATGACTCTGTCATCATTAAAAGCAAACAGGATCCGCTCCGCCTGACCATCACTATTCAGATCTTGCTCGACCAGGACGCAGGCGTCTTTTTCAATACAGGTCATCACGTTATAGCGATCCTGGATTAAGGCCGACCAGAATGTCGCATCA  >  NC_000913/3530759‑3531048

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑