Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,462,677 A→T 5.4% 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.4% 321.6 / 6.3 97intergenic (+31/‑335)fadL/yfdFlong‑chain fatty acid outer membrane transporter/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (37/51);  new base T (2/3);  total (41/56)
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.95e-01

ACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAT  >  NC_000913/2462531‑2462821
                                                                                                                                                  |                                                                                                                                                 
aCGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCAt                                                                                                                                                  >  1:178637/1‑148 (MQ=255)
     ggTGTTACTTATATGTACGGTCAGAGCGTGAAAATTCACGAAGGCCCATTCCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAg                                                                                                                                            <  1:752444/148‑1 (MQ=255)
      gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          <  3:371539/149‑1 (MQ=255)
      gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          >  3:473865/1‑149 (MQ=255)
      gtgtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          <  4:131005/149‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTg                                                                                                                                         <  1:121739/148‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTg                                                                                                                                         <  1:606174/148‑1 (MQ=255)
        gtTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCAGACGAGTTCGAGTCTGAAGGTAAAGTCTGGCTGGTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGTCTGGATAAAGTCACCTGCATAGCAGGTg                                                                                                                                         >  2:53652/1‑148 (MQ=255)
         tttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGa                                                                                                                                        <  2:33840/148‑1 (MQ=255)
         tttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGAc                                                                                                                                       <  3:275508/149‑1 (MQ=255)
          ttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGAc                                                                                                                                       <  1:109998/148‑1 (MQ=255)
          ttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGAc                                                                                                                                       <  2:500550/148‑1 (MQ=255)
          ttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGAc                                                                                                                                       >  2:97796/1‑148 (MQ=255)
          ttCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGAc                                                                                                                                       >  1:782708/1‑148 (MQ=255)
             ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCTATGCAGGTGACTTTa                                                                                                                                   <  3:747648/149‑1 (MQ=255)
             ttATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTa                                                                                                                                   >  3:492497/1‑149 (MQ=255)
              tataTGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTaa                                                                                                                                  <  4:48685/149‑1 (MQ=255)
              tataTGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTaa                                                                                                                                  <  4:739434/149‑1 (MQ=255)
              tataTGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTa                                                                                                                                   <  2:789660/148‑1 (MQ=255)
                tatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAAc                                                                                                                                 >  2:393968/1‑148 (MQ=255)
                tatGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAAc                                                                                                                                 >  1:657730/1‑148 (MQ=255)
                  tGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTc                                                                                                                               <  1:899422/148‑1 (MQ=255)
                    cACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTccc                                                                                                                             <  1:732140/148‑1 (MQ=255)
                     aCGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTccccc                                                                                                                           >  3:696557/1‑149 (MQ=255)
                           agagCGTGAAATTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGTCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCa                                                                                                                     >  4:88974/1‑149 (MQ=255)
                               cGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCg                                                                                                                  >  2:695330/1‑148 (MQ=255)
                                 tGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTTGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaaa                                                                                                               <  3:786069/149‑1 (MQ=255)
                                 tGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGaa                                                                                                                <  1:842923/148‑1 (MQ=255)
                                    aaaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagta                                                                                                            <  4:696557/149‑1 (MQ=255)
                                    aaaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagt                                                                                                             <  1:676756/148‑1 (MQ=255)
                                    aaaTTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAagt                                                                                                             <  1:285203/148‑1 (MQ=255)
                                           cGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTcc                                                                                                      <  2:771886/148‑1 (MQ=255)
                                                   cATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCt                                                                                              <  2:731808/148‑1 (MQ=255)
                                                        cAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACttt                                                                                                                                    <  4:937313/105‑1 (MQ=255)
                                                        cAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACttt                                                                                                                                    >  3:937313/1‑105 (MQ=255)
                                                             cGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAg                                                                                   >  3:875795/1‑149 (MQ=255)
                                                               aGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGt                                                                                  >  2:535880/1‑148 (MQ=255)
                                                                gTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGttt                                                                                <  3:213456/149‑1 (MQ=255)
                                                                 tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGtttt                                                                               <  4:875795/149‑1 (MQ=255)
                                                                 tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGttt                                                                                <  1:894530/148‑1 (MQ=255)
                                                                 tCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGttt                                                                                <  2:76353/148‑1 (MQ=255)
                                                                   tGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGc                                                                             >  4:378626/1‑149 (MQ=255)
                                                                    gAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCa                                                                            <  4:492497/149‑1 (MQ=255)
                                                                        gTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACa                                                                        >  4:845676/1‑149 (MQ=255)
                                                                          aaaGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGc                                                                      >  3:96552/1‑149 (MQ=255)
                                                                            aGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGcc                                                                     <  1:99759/148‑1 (MQ=255)
                                                                              ccTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGt                                                                  <  3:904085/149‑1 (MQ=255)
                                                                               cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                 >  4:115221/1‑149 (MQ=255)
                                                                               cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGtt                                                                 <  3:318254/149‑1 (MQ=255)
                                                                               cTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGt                                                                  <  1:535880/148‑1 (MQ=255)
                                                                                tGGCTGTTCGGTACTAACTTTAACTACTCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACGTGCATAGCAGGTGACTGTGACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGta                                                                 >  2:692221/1‑147 (MQ=255)
                                                                                tGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGTTAAAGTCACCTGCTATGCAGGTGACTTTatcca                                                                                                                               <  2:135824/86‑5 (MQ=25)
                                                                                tGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGTTAAAGTCACCTGCTATGCAGGTGACTTTatcca                                                                                                                               >  1:135824/1‑82 (MQ=25)
                                                                                 ggCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTg                                                                <  1:554132/148‑1 (MQ=255)
                                                                                  gCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTgg                                                               <  2:573719/148‑1 (MQ=255)
                                                                                    tGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          >  3:649655/1‑71 (MQ=255)
                                                                                    tGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          <  4:649655/71‑1 (MQ=255)
                                                                                     gTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAg                                                            >  2:439684/1‑148 (MQ=255)
                                                                                     gTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAg                                                            >  1:611289/1‑148 (MQ=255)
                                                                                       tCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCt                                                          <  2:420715/148‑1 (MQ=255)
                                                                                             cTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                   >  3:869159/1‑149 (MQ=255)
                                                                                              tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCttt                                                  <  4:538507/149‑1 (MQ=255)
                                                                                              tAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCtt                                                   >  2:210647/1‑148 (MQ=255)
                                                                                               aaCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTc                                                 >  3:813656/1‑149 (MQ=255)
                                                                                                aCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAt                                                                                                  >  2:654769/1‑99 (MQ=255)
                                                                                                aCTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCAt                                                                                                  <  1:654769/99‑1 (MQ=255)
                                                                                                 cTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                               >  3:730883/1‑149 (MQ=255)
                                                                                                 cTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaa                                               >  4:337855/1‑149 (MQ=255)
                                                                                                   ttAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCaaa                                              >  2:946167/1‑148 (MQ=255)
                                                                                                         aCGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          <  3:832385/50‑1 (MQ=255)
                                                                                                         aCGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGt                                                                                                                                          >  4:832385/1‑50 (MQ=255)
                                                                                                         aCGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATa                                        <  2:381192/148‑1 (MQ=255)
                                                                                                           gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAg                                      <  1:45299/148‑1 (MQ=255)
                                                                                                           gcgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGc                                     <  3:213707/149‑1 (MQ=255)
                                                                                                            cgTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCg                                    <  4:730883/149‑1 (MQ=255)
                                                                                                               tCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGt                                  <  2:782708/148‑1 (MQ=255)
                                                                                                               tCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGt                                  >  2:625288/1‑148 (MQ=255)
                                                                                                                     aaCGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCcac                            <  1:593021/148‑1 (MQ=255)
                                                                                                                      aCGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAg                          >  3:688735/1‑149 (MQ=255)
                                                                                                                       cgcgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTa                                                                                            >  2:276781/1‑82 (MQ=255)
                                                                                                                       cgcgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTa                                                                                            <  1:276781/82‑1 (MQ=255)
                                                                                                                         cgTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGAc                        >  1:323922/1‑148 (MQ=255)
                                                                                                                           ttCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaa                      <  1:562461/148‑1 (MQ=255)
                                                                                                                            tCGCCTGGATAAAGTCACCTGGATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCACCTGAGTTTTGCAAACAGGCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaaa                     <  2:240936/148‑1 (MQ=255)
                                                                                                                            tCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACaaa                     <  1:393968/148‑1 (MQ=255)
                                                                                                                                cTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTg                 <  2:73539/148‑1 (MQ=255)
                                                                                                                                 tGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGaa               <  4:473865/149‑1 (MQ=255)
                                                                                                                                   gATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAg              >  2:771257/1‑148 (MQ=255)
                                                                                                                                    aTAAAGTCACCTGCGTAGCGGGTGACTTTAACGCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGAATCCACAGACAAACCTGAAggg            >  4:213543/1‑149 (MQ=255)
                                                                                                                                      aaaGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaa          <  3:911698/149‑1 (MQ=255)
                                                                                                                                        aGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGaaa         <  2:178637/148‑1 (MQ=255)
                                                                                                                                          tCACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGc      >  3:815795/1‑149 (MQ=255)
                                                                                                                                           cACCTGCTATGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGGGTCCTCATTGCTTATCTCCTGAGTTTTGCAAATAGCCTGTTGGCAGCTTGCCCTTCCAAATCAATAAGCGGTGTCCACAGACAAACCTGAAGGGAAACGc      >  1:495914/1‑148 (MQ=255)
                                                                                                                                            aCCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCat    >  4:586560/1‑149 (MQ=255)
                                                                                                                                            aCCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCa     <  1:226491/148‑1 (MQ=255)
                                                                                                                                              cTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCatat  >  4:102183/1‑149 (MQ=255)
                                                                                                                                              cTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCata   <  1:771257/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACGTTGGTGTTTCTTATATGCACGGTCAGAGCGTGAAAATTAACGAAGGCCCATACCAGTTCGAGTCTGAAGGTAAAGCCTGGCTGTTCGGTACTAACTTTAACTACGCGTTCTGATAACGCGTTCGCCTGGATAAAGTCACCTGCATAGCAGGTGACTTTAACTCCCCCACTTCACCGAAAGTAGTGTCCTCATTGCTTACCTCCTGAGTTTTGCAAACAGCCTGTTGGCAGCTTGCCCTTTCAAATCAATAAGCGGTATCCACAGACAAACCTGAAGGGAAAGGCATAT  >  NC_000913/2462531‑2462821

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: