Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 752,066 T→G 6.4% D35A (GAC→GCC) ‡ ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0660TG 6.4% 275.3 / 9.2 78D35A (GAC→GCC) ‡ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (32/41);  new base G (2/3);  total (34/44)
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.93e-01

TTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGC  >  NC_000913/751919‑752214
                                                                                                                                                   |                                                                                                                                                    
ttGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGt                                                                                                                                                      >  1:444921/1‑148 (MQ=255)
       aTGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAtt                                                                                                                                              <  3:454391/149‑1 (MQ=255)
       aTGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAt                                                                                                                                               <  2:149655/148‑1 (MQ=255)
       aTGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAt                                                                                                                                               <  2:724315/148‑1 (MQ=255)
        tGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAt                                                                                                                                               <  3:134571/147‑1 (MQ=255)
        tGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAt                                                                                                                                               >  4:134571/1‑147 (MQ=255)
         ggAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAtt                                                                                                                                              <  4:606620/147‑1 (MQ=255)
         ggAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCAtt                                                                                                                                              >  3:606620/1‑147 (MQ=255)
          gAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGa                                                                                                                                            <  2:2522/148‑1 (MQ=255)
          gAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGa                                                                                                                                            <  1:309971/148‑1 (MQ=255)
           aaTGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAta                                                                                                                                          <  4:259838/149‑2 (MQ=255)
              ggTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTa                                                                                                                                        <  2:444921/148‑1 (MQ=255)
               gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATc                                                                                                                                      <  3:723785/149‑1 (MQ=255)
               gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATc                                                                                                                                      >  4:672983/1‑149 (MQ=255)
                ttttGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCg                                                                                                                                     >  3:17063/1‑149 (MQ=255)
                        cccTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACAt                                                                                                                              >  1:397588/1‑148 (MQ=255)
                        cccTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACAt                                                                                                                              >  2:417252/1‑148 (MQ=255)
                          cTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTc                                                                                                                            >  1:55791/1‑148 (MQ=255)
                          cTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCa                                                                                                                           >  4:505083/1‑149 (MQ=255)
                           tttttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAg                                                                                                                          >  4:641134/1‑149 (MQ=255)
                           tttttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGa                                                                                                                                            >  1:672970/1‑131 (MQ=255)
                            ttttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAta                                                                                                                                          >  2:278602/1‑131 (MQ=255)
                            ttttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAta                                                                                                                                          <  1:278602/132‑2 (MQ=255)
                             tttttGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGa                                                                                                                         <  2:290162/148‑1 (MQ=255)
                                  gCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAAt                                                                                                                   <  4:892497/149‑1 (MQ=255)
                                  gCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAAt                                                                                                                   <  4:178779/149‑1 (MQ=255)
                                   ccAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAAt                                                                                                                   >  2:768513/1‑148 (MQ=255)
                                    cAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTc                                                                                                                 <  4:17063/149‑1 (MQ=255)
                                     aGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTc                                                                                                                 <  2:397588/148‑1 (MQ=255)
                                       ttAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCga                                                                                                                                    >  4:540200/1‑127 (MQ=255)
                                       ttAATTTTGTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCga                                                                                                                                    <  3:540200/127‑1 (MQ=255)
                                              tttGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAt                                                                                                        <  1:439583/148‑1 (MQ=255)
                                              tttGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAt                                                                                                        <  1:417252/148‑1 (MQ=255)
                                                      ttACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAgg                                                                                               >  3:105386/1‑149 (MQ=255)
                                                          tGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCCTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCAt                                                                                            >  1:81360/1‑148 (MQ=255)
                                                           gATTTTGTGGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATACAGATTGAATTCGACCGCCATTGCGCAAGGCATc                                                                                           >  1:783329/1‑148 (MQ=255)
                                                              tttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCaa                                                                                                 <  3:25449/139‑1 (MQ=255)
                                                              tttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCaa                                                                                                 >  4:25449/1‑139 (MQ=255)
                                                              tttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGcc                                                                                        >  1:879458/1‑148 (MQ=255)
                                                              tttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGcc                                                                                        <  1:764922/148‑1 (MQ=255)
                                                               ttgttgTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCa                                                                                       <  1:768513/148‑1 (MQ=255)
                                                                 gttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATg                                                                                     <  2:857550/148‑1 (MQ=255)
                                                                  ttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                                    <  1:664055/148‑1 (MQ=255)
                                                                  ttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                                    >  1:569338/1‑148 (MQ=255)
                                                                         tCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggca                                                                             >  1:375098/1‑148 (MQ=255)
                                                                          cACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACcaggcag                                                                            <  2:879458/148‑1 (MQ=255)
                                                                            ccGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGa                                                                          >  1:76290/1‑148 (MQ=255)
                                                                                aaCGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACaa                                                                     >  4:699703/1‑149 (MQ=255)
                                                                                  cGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACaaaa                                                                   <  3:505083/149‑1 (MQ=255)
                                                                                      aaCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                                >  2:310654/1‑148 (MQ=255)
                                                                                      aaCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagag                                                                >  2:499551/1‑148 (MQ=255)
                                                                                       aCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAagaga                                                               <  2:81360/148‑1 (MQ=255)
                                                                                        cAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGt                                                             <  3:789551/149‑1 (MQ=255)
                                                                                         aaaaTATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGt                                                             <  2:76290/148‑1 (MQ=255)
                                                                                         aaaaTATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGt                                                             <  2:569338/148‑1 (MQ=255)
                                                                                               ttCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATaaa                                                      >  3:844791/1‑149 (MQ=255)
                                                                                                     gggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTc                                                 <  1:465269/148‑1 (MQ=255)
                                                                                                      ggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCa                                                <  1:310654/148‑1 (MQ=255)
                                                                                                      ggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCAc                                               <  3:672983/149‑1 (MQ=255)
                                                                                                          cAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAata                                                    >  3:804820/1‑140 (MQ=255)
                                                                                                          cAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAata                                                    <  4:804820/140‑1 (MQ=255)
                                                                                                              gACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTcc                                        >  1:420289/1‑148 (MQ=255)
                                                                                                               aCTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACAGTGTCCAt                                      <  3:835042/149‑1 (MQ=255)
                                                                                                                     ttttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTg                                 <  1:499551/148‑1 (MQ=255)
                                                                                                                     ttttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGa                                <  3:641134/149‑1 (MQ=255)
                                                                                                                      tttCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGa                                <  2:375098/148‑1 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATAcc                                   >  1:65841/1‑140 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATAcc                                   <  2:65841/140‑1 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAata                          >  4:527493/1‑149 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAata                          >  3:277642/1‑149 (MQ=255)
                                                                                                                              aataatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATAttt                       <  3:699703/149‑1 (MQ=255)
                                                                                                                                 aatGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTAt                    >  3:594537/1‑149 (MQ=255)
                                                                                                                                   tGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATCGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGa                  >  4:617833/1‑149 (MQ=255)
                                                                                                                                       atatCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATaaa                                                      <  2:146374/109‑1 (MQ=255)
                                                                                                                                       atatCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATaaa                                                      >  1:146374/1‑109 (MQ=255)
                                                                                                                                       atatCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAgg              <  4:844791/149‑1 (MQ=255)
                                                                                                                                        tatCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAATAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAgg              <  2:783329/148‑1 (MQ=255)
                                                                                                                                            aaTGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGgcggcg         >  3:151398/1‑149 (MQ=255)
                                                                                                                                              tGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCgat       <  4:753577/149‑1 (MQ=255)
                                                                                                                                                   tCGCGCATTGATTTATCGAGAACATTCAGATTGAAt                                                                                                                   >  4:306763/1‑36 (MQ=255)
                                                                                                                                                   tCGCGCATTGATTTATCGAGAACATTCAGATTGAAt                                                                                                                   <  3:306763/36‑1 (MQ=255)
                                                                                                                                                   tCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATgc  <  4:105386/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGATGATGC  >  NC_000913/751919‑752214

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: