Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 752,065 G→C 6.3% D35E (GAC→GAG) ‡ ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0650GC 6.3% 279.2 / 7.7 79D35E (GAC→GAG) ‡ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (33/41);  new base C (2/3);  total (35/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

TTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGAT  >  NC_000913/751919‑752209
                                                                                                                                                  |                                                                                                                                                
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:417252/148‑1 (MQ=255)
                                              tttGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAt                                                                                                   <  1:439583/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                                                                                            >  4:25449/1‑139 (MQ=255)
                                                              tttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCaa                                                                                            <  3:25449/139‑1 (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:569338/1‑148 (MQ=255)
                                                                  ttgttATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGa                                                                               <  1:664055/148‑1 (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                                               <  4:804820/140‑1 (MQ=255)
                                                                                                          cAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAata                                               >  3:804820/1‑140 (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                              <  2:65841/140‑1 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATAcc                              >  1:65841/1‑140 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAata                     >  3:277642/1‑149 (MQ=255)
                                                                                                                           tttAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAata                     >  4:527493/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)
                                                                                                                                                  |                                                                                                                                                
TTGATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGCGAT  >  NC_000913/751919‑752209

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: ‑