Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,677 A→T 5.3% intergenic (+31/‑335) fadL → / → yfdF long‑chain fatty acid outer membrane transporter/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,462,6770AT 5.3% 337.3 / 6.4 98intergenic (+31/‑335)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (44/46);  new base T (2/3);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.21e-01

GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTT  >  NC_000913/2462530‑2462823
                                                                                                                                                   |                                                                                                                                                   
gACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCAt                                                                                                                                                    >  4:614365/1‑149 (MQ=255)
 aCGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATa                                                                                                                                                   <  3:459947/149‑1 (MQ=255)
    ttGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                               <  1:197314/149‑1 (MQ=255)
     tGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                              <  1:395560/149‑1 (MQ=255)
     tGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAg                                                                                                                                              <  3:614365/149‑1 (MQ=255)
       gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGt                                                                                                                                            <  1:110694/149‑1 (MQ=255)
        tgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTg                                                                                                                                           >  2:19378/1‑149 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGa                                                                                                                                          <  1:878354/149‑1 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGa                                                                                                                                          <  1:669461/149‑1 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                          <  1:19378/149‑1 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                          <  4:181394/149‑1 (MQ=255)
         gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGC‑TATGCAGGTGa                                                                                                                                          <  2:687574/149‑1 (MQ=255)
          tttCTTATATGCACGGTCAGAGCGTGTGAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTAGGTACTAACTTTAACTACGCGTTCTGATAACTCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGAc                                                                                                                                         >  1:142654/1‑149 (MQ=255)
           ttCCTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGACt                                                                                                                                        <  4:484341/149‑1 (MQ=255)
            tCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACtt                                                                                                                                       <  3:476726/149‑1 (MQ=255)
             cTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACttt                                                                                                                                      >  3:447026/1‑149 (MQ=255)
              ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                               >  3:367886/1‑139 (MQ=255)
              ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCa                                                                                                                                               <  4:367886/139‑1 (MQ=255)
              ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTa                                                                                                                                     >  4:249483/1‑149 (MQ=255)
               tataTGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTATTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTaa                                                                                                                                    <  2:178538/149‑1 (MQ=255)
                atatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTAT‑GCAGGTGACTTTatc                                                                                                                                   <  4:498061/149‑3 (MQ=255)
                atatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAAc                                                                                                                                   <  1:458819/149‑1 (MQ=255)
                     cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcccc                                                                                                                              <  1:414307/149‑1 (MQ=255)
                     cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTcccc                                                                                                                              <  2:424900/149‑1 (MQ=255)
                                 gTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                  >  3:38155/1‑149 (MQ=255)
                                  tGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGaaa                                                                                                                 <  1:581082/149‑1 (MQ=255)
                                   gAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAag                                                                                                                >  2:378314/1‑149 (MQ=255)
                                   gAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAag                                                                                                                >  3:797841/1‑149 (MQ=255)
                                       aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                            <  3:653361/149‑1 (MQ=255)
                                       aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                            >  1:8654/1‑149 (MQ=255)
                                       aTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAagtagt                                                                                                            <  1:658266/149‑1 (MQ=255)
                                           aCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTcc                                                                                                        >  4:136977/1‑149 (MQ=255)
                                           aCGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTcc                                                                                                        >  4:441371/1‑149 (MQ=255)
                                                   ccATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCt                                                                                                >  3:297428/1‑149 (MQ=255)
                                                    cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCtt                                                                                               >  3:168503/1‑149 (MQ=255)
                                                    cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCtt                                                                                               <  4:11750/149‑1 (MQ=255)
                                                      tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                   <  2:125375/143‑1 (MQ=255)
                                                      tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAtt                                                                                                   >  1:125375/1‑143 (MQ=255)
                                                      tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAc                                                                                             <  3:442737/149‑1 (MQ=255)
                                                      tACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGCCACCTGCATA‑GCAGGTGACTTTAACTCCACCACTTCACCGAAAGTAGTGTCCTCATTGCTTAc                                                                                             <  2:613998/149‑1 (MQ=255)
                                                       aCCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCTTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTAcc                                                                                            <  2:227943/149‑1 (MQ=255)
                                                           gTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATTACGTGTGCGCCGGGATAAAGTCACCAGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCTAAAGTAGTGTCCTCATTGCTTAcctcct                                                                                        <  2:638658/149‑1 (MQ=255)
                                                             tCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                      <  3:615048/149‑1 (MQ=255)
                                                             tCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGa                                                                                      <  3:819310/149‑1 (MQ=255)
                                                               gAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGt                                                                                    >  1:16099/1‑149 (MQ=255)
                                                                       aGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                            >  1:504083/1‑149 (MQ=255)
                                                                          tAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAg                                                                         >  2:303870/1‑149 (MQ=255)
                                                                           aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                        <  1:869403/149‑1 (MQ=255)
                                                                           aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                        <  1:340778/149‑1 (MQ=255)
                                                                             aGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCt                                                                      <  4:447026/149‑1 (MQ=255)
                                                                                cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                   <  1:378314/149‑1 (MQ=255)
                                                                                  ggCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                                 >  3:593102/1‑149 (MQ=255)
                                                                                  ggCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                                 >  2:367291/1‑149 (MQ=255)
                                                                                        tCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCtt                                                           <  1:252474/149‑1 (MQ=255)
                                                                                           gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                        <  1:110110/149‑1 (MQ=255)
                                                                                           gTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGcc                                                        >  4:524101/1‑149 (MQ=255)
                                                                                              cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                            >  4:18168/1‑126 (MQ=255)
                                                                                              cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                            <  3:18168/126‑1 (MQ=255)
                                                                                              cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                     >  1:510754/1‑149 (MQ=255)
                                                                                              cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                     >  3:595337/1‑149 (MQ=255)
                                                                                              cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                     >  1:87542/1‑149 (MQ=255)
                                                                                               tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                    >  3:842864/1‑149 (MQ=255)
                                                                                                  cTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                                 >  3:267211/1‑149 (MQ=255)
                                                                                                   tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                >  3:397935/1‑149 (MQ=255)
                                                                                                   tttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                                >  1:730418/1‑149 (MQ=255)
                                                                                                       aCTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCaa                                            <  3:249483/149‑1 (MQ=255)
                                                                                                        cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAAt                                           <  4:297428/149‑1 (MQ=255)
                                                                                                        cTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAAt                                           <  2:470190/149‑1 (MQ=255)
                                                                                                         tACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCGGCTTGCCCTTTCAAATCAATa                                          >  1:844497/1‑149 (MQ=255)
                                                                                                              gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCgg                                     >  1:362429/1‑149 (MQ=255)
                                                                                                              gTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCgg                                     >  4:61736/1‑149 (MQ=255)
                                                                                                                  tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATc                                              <  2:209377/136‑1 (MQ=255)
                                                                                                                  tGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATc                                              >  1:209377/1‑136 (MQ=255)
                                                                                                                      aaCGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcaca                             >  1:297366/1‑149 (MQ=255)
                                                                                                                             tCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAAc                      <  3:461006/149‑1 (MQ=255)
                                                                                                                                 cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGGAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                            >  2:767310/1‑91 (MQ=255)
                                                                                                                                 cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGa                  >  3:252624/1‑149 (MQ=255)
                                                                                                                                 cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGa                  >  2:549722/1‑149 (MQ=255)
                                                                                                                                 cTGGATAAAGTCACCTGCATA‑GCAGGTGACTTGAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCaaa                                                                            <  1:767310/91‑1 (MQ=255)
                                                                                                                                  tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                 >  1:525345/1‑149 (MQ=255)
                                                                                                                                  tGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa                 <  2:586108/149‑1 (MQ=255)
                                                                                                                                   ggATAAAGTCACCTGCATA‑GCAGGTGCCTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg                <  2:142654/149‑1 (MQ=255)
                                                                                                                                   ggATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg                <  1:303870/149‑1 (MQ=255)
                                                                                                                                    gATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                                 <  4:596072/131‑1 (MQ=255)
                                                                                                                                    gATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATc                                 >  3:596072/1‑131 (MQ=255)
                                                                                                                                     aTAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAggg              >  1:175226/1‑149 (MQ=255)
                                                                                                                                      tAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGa             <  4:397935/149‑1 (MQ=255)
                                                                                                                                       aaaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa            >  3:225062/1‑149 (MQ=255)
                                                                                                                                        aaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa           >  4:701763/1‑149 (MQ=255)
                                                                                                                                        aaGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa           <  1:175229/149‑1 (MQ=255)
                                                                                                                                         aGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg          <  2:643126/149‑1 (MQ=255)
                                                                                                                                         aGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg          <  1:232745/149‑1 (MQ=255)
                                                                                                                                         aGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAg          <  3:45234/149‑1 (MQ=255)
                                                                                                                                           tCACCTGCTAT‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGc        >  4:736711/1‑149 (MQ=255)
                                                                                                                                             aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATcc                                <  1:269385/123‑1 (MQ=255)
                                                                                                                                             aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATcc                                >  2:269385/1‑123 (MQ=255)
                                                                                                                                             aCCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat      >  4:499962/1‑149 (MQ=255)
                                                                                                                                                 gCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTACCCTTTCAAATGAATAAGCGGTATCCACAGAGAAACCTGAAGGGAATGGCATAttt  >  2:117598/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATA‑GCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATATTT  >  NC_000913/2462530‑2462823

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