Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,464,224 T→C 5.1% intergenic (+154/+28) yfdF → / ← mlaA uncharacterized protein/ABC transporter maintaining OM lipid asymmetry, OM lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,464,2240TC 5.1% 1091.8 / 30.0 316intergenic (+154/+28)yfdF/mlaAuncharacterized protein/ABC transporter maintaining OM lipid asymmetry, OM lipoprotein component
Reads supporting (aligned to +/- strand):  ref base T (160/139);  new base C (3/13);  total (164/152)
Fisher's exact test for biased strand distribution p-value = 8.73e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 324 total aligned reads displayed.

TTTTTCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCACC  >  NC_000913/2464076‑2464372
                                                                                                                                                    |                                                                                                                                                    
tttttCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGttt                                                                                                                                                     >  1:350619/1‑150 (MQ=255)
 ttttCCAGTTTTTTTATTTTTCAACAGCTGAATATCGGAAATAGTAAGTGACGAACATCAGTATTTTTAAATAATTAAGTCTAGGTCATTTTCGGTTTATCTCCATTTAATTTTTGTAATTGAGGGATAACAAAAAAAAATGTTATtttt                                                                                                                                                    >  1:618232/1‑150 (MQ=255)
     ccAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGa                                                                                                                                                <  2:1823192/150‑1 (MQ=255)
     ccAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGa                                                                                                                                                >  2:4221666/1‑150 (MQ=255)
     ccAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGa                                                                                                                                                >  2:4091615/1‑150 (MQ=255)
      cAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGAc                                                                                                                                               >  2:6639430/1‑150 (MQ=255)
      cAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGAc                                                                                                                                               >  2:4649699/1‑150 (MQ=255)
       tgTTTTTTTATTTTTCAACAGCTGAATATCGTAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACt                                                                                                                                              <  1:5043761/149‑1 (MQ=255)
        gTTTTTTGATTTTTCAACAGCTGAATATCGTAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTc                                                                                                                                             <  1:2484972/150‑1 (MQ=255)
         ttttttGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCa                                                                                                                                            <  2:2098429/150‑1 (MQ=255)
          tttttGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCAc                                                                                                                                           <  1:706135/150‑1 (MQ=255)
              tGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCttt                                                                                                                                       <  2:3413640/150‑1 (MQ=255)
               gATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCtttt                                                                                                                                      >  2:6073108/1‑150 (MQ=255)
               gATTTTTCAACAGCTGAATATCGGAAATAGTCAGGGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCtttt                                                                                                                                      <  2:5329174/150‑1 (MQ=255)
                aTTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCttttt                                                                                                                                     >  1:4620979/1‑150 (MQ=255)
                aTTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCttttt                                                                                                                                     <  1:4222091/150‑1 (MQ=255)
                aTTTGTCAACAGCTGACTATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCttttt                                                                                                                                     <  2:5616126/150‑1 (MQ=255)
                 tttttCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCtttttt                                                                                                                                    <  2:2975677/150‑1 (MQ=255)
                     tCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTAttt                                                                                                                                <  2:3231759/150‑1 (MQ=255)
                     tCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTAttt                                                                                                                                >  2:511589/1‑150 (MQ=255)
                     tCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCTTTTTTGttg                                                                                                                                <  2:7171935/150‑2 (MQ=39)
                     tCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCTTTTTTGttg                                                                                                                                <  1:4461711/150‑2 (MQ=39)
                      cAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTCGCAAAACTCACCTTTTTTGttgt                                                                                                                               <  1:1376355/150‑3 (MQ=38)
                       caCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGt                                                                                                                              <  2:4752209/149‑1 (MQ=255)
                       aaCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGt                                                                                                                              <  1:5633143/150‑1 (MQ=255)
                       aaCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGt                                                                                                                              >  1:5286768/1‑150 (MQ=255)
                        aCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGtt                                                                                                                             <  2:3925508/150‑1 (MQ=255)
                        aCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGtt                                                                                                                             <  1:2663219/150‑1 (MQ=255)
                        aCAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGtt                                                                                                                             <  1:585146/150‑1 (MQ=255)
                         cAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGttt                                                                                                                            <  2:1270050/150‑1 (MQ=255)
                         cAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGttt                                                                                                                            <  1:2729577/150‑1 (MQ=255)
                          aGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTc                                                                                                                           >  1:3694267/1‑150 (MQ=255)
                          aGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTc                                                                                                                           >  1:911722/1‑150 (MQ=255)
                            cTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCtt                                                                                                                         <  1:3338820/150‑1 (MQ=255)
                             tGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTa                                                                                                                        <  1:2847964/150‑1 (MQ=255)
                             tGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTa                                                                                                                        >  2:3286627/1‑150 (MQ=255)
                                  atCGGAAATAGTCAGTGACGAACATCCGTATTTTTACATACTTCCGTCTAGGTCATTGTCGGTTTATCTCCATTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAg                                                                                                                   <  2:2519233/150‑1 (MQ=255)
                                     ggAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAAt                                                                                                                >  1:5210738/1‑150 (MQ=255)
                                           aGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAatat                                                                                                          >  2:6203540/1‑150 (MQ=255)
                                             tCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCt                                                                                                        >  2:5713877/1‑150 (MQ=255)
                                               aGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCttt                                                                                                      >  2:3594131/1‑150 (MQ=255)
                                               aGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCttt                                                                                                      <  2:7045407/150‑1 (MQ=255)
                                               aGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCttt                                                                                                      >  2:6678150/1‑150 (MQ=255)
                                                gTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCtttt                                                                                                     <  1:1932841/150‑1 (MQ=255)
                                                 tGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTa                                                                                                    <  2:1712492/150‑1 (MQ=255)
                                                 tGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTa                                                                                                    <  2:2757383/150‑1 (MQ=255)
                                                  gACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTaa                                                                                                   <  2:6521957/150‑1 (MQ=255)
                                                     gAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAatc                                                                                                <  1:4338476/150‑1 (MQ=255)
                                                     gAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAata                                                                                                >  1:6692608/1‑149 (MQ=255)
                                                        aaTCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAatcatc                                                                                             <  2:5165212/149‑1 (MQ=255)
                                                         aTCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTGTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATcc                                                                                            >  1:2137982/1‑150 (MQ=255)
                                                         aTCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGAGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATcc                                                                                            >  2:1024660/1‑150 (MQ=255)
                                                             gTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGaa                                                                                        >  2:6876840/1‑150 (MQ=255)
                                                              tATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTTTTTCTTATTCAGAATAAATATATTTTAAATCATCCTGAAt                                                                                       >  1:4160696/1‑150 (MQ=255)
                                                                tttttACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCg                                                                                     <  2:5286768/150‑1 (MQ=255)
                                                                  tttACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTAACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCt                                                                                   >  1:1006664/1‑150 (MQ=255)
                                                                       aTAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGc                                                                              <  1:4056802/150‑1 (MQ=255)
                                                                           ttACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTTGGAATTGAGGGATAACAACAAAAAAGGTGAGGTTTGCGACTCACATTTTTTATTTGTTTCTGATGCATAATCAATATATTTTAAATCATCCTGAATAGATTGTGCGTTc                                                                          >  2:232538/1‑150 (MQ=255)
                                                                           ttACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTc                                                                          <  1:4921082/150‑1 (MQ=255)
                                                                            tACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCg                                                                         <  2:4515548/150‑1 (MQ=255)
                                                                                tCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGAtt                                                                     >  2:4042881/1‑150 (MQ=255)
                                                                                   aGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTc                                                                  <  1:4046021/150‑1 (MQ=255)
                                                                                   aGGTCATTGTCGGTTTATATCCTTTTACGTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTATTATTCAGAATCACTCTCTTTTAAATAATCCTGAATCGCTTGTGCGTTCGGAttttg                                                                  >  1:6943999/1‑149 (MQ=255)
                                                                                     gTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTAGTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCt                                                                >  2:6443144/1‑150 (MQ=255)
                                                                                         ttGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCgg                                                            <  1:832253/150‑1 (MQ=255)
                                                                                            tCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGttt                                                         >  2:4951351/1‑150 (MQ=255)
                                                                                              ggTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGa                                                       <  1:6409433/150‑1 (MQ=255)
                                                                                              ggTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGa                                                       <  1:3228335/150‑1 (MQ=255)
                                                                                              ggTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGa                                                       <  1:1160874/150‑1 (MQ=255)
                                                                                                 ttATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTGTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGtt                                                    <  1:7076394/150‑1 (MQ=255)
                                                                                                  tATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTc                                                   <  2:3694267/150‑1 (MQ=255)
                                                                                                   aTCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcg                                                  >  1:1556700/1‑150 (MQ=255)
                                                                                                    tctcCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgc                                                 >  2:2914473/1‑150 (MQ=255)
                                                                                                     ctcCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgcc                                                >  1:4406183/1‑150 (MQ=255)
                                                                                                      tcCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGAGACTCACCTTTTTTATTTGTTTCTTATTCAGAACCAATATCTTTTAAATCATAATGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTagccg                                               >  2:3771801/1‑150 (MQ=255)
                                                                                                       ccTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTcgccgc                                              <  1:4318500/150‑1 (MQ=255)
                                                                                                          tttACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCAt                                           >  1:2574541/1‑150 (MQ=255)
                                                                                                              cTTGTGGTAATTGAGGGATAACAACAAACAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGc                                       <  2:5850713/150‑1 (MQ=255)
                                                                                                               ttGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCg                                      <  1:3194226/150‑1 (MQ=255)
                                                                                                                   ggTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGa                                  >  1:3074465/1‑150 (MQ=255)
                                                                                                                   ggTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGa                                  >  1:6659087/1‑150 (MQ=255)
                                                                                                                   ggTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGATCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGa                                  >  2:4286110/1‑150 (MQ=255)
                                                                                                                   ggTAATTGAGGGAAAACAACAAAAAAGGTGAGTTTTGCGACTCACCTATTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAAACATACTGAATCGCTAGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGa                                  >  1:1987607/1‑150 (MQ=255)
                                                                                                                       aTTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATc                              <  1:715102/150‑1 (MQ=255)
                                                                                                                       aTTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATc                              >  1:1629751/1‑150 (MQ=255)
                                                                                                                        ttGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCa                             <  1:5891533/150‑1 (MQ=255)
                                                                                                                        ttGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCa                             <  1:2773698/150‑1 (MQ=255)
                                                                                                                         tGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCAt                            >  1:97531/1‑150 (MQ=255)
                                                                                                                           aGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGa                          >  2:173250/1‑150 (MQ=255)
                                                                                                                               aTAACAACAAAAAAGGTGAGTTTTGCGAATCACCTTTTTTATTTGTTTCTTATTCAGATTCTATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGTTTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCt                      >  2:596514/1‑150 (MQ=255)
                                                                                                                                tAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTg                     <  2:1858061/150‑1 (MQ=255)
                                                                                                                                 acaaataaaaaaGGTGAGTCGCAAAACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCGTGTGCGTTCGGATTTTCCTGCGGTGTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTgg                    >  2:5554487/7‑150 (MQ=38)
                                                                                                                                 aacaacAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTgg                    >  2:3479078/1‑150 (MQ=255)
                                                                                                                                    aacaaAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAg                 >  2:1559994/1‑150 (MQ=255)
                                                                                                                                    aacaaAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAg                 >  1:6486449/1‑150 (MQ=255)
                                                                                                                                       aaaaaaGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTAc              >  2:2566709/1‑150 (MQ=255)
                                                                                                                                         aaaaGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGc            >  2:5518635/1‑150 (MQ=255)
                                                                                                                                          aaaGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCt           <  1:1854280/150‑1 (MQ=255)
                                                                                                                                                   ttttGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGTAATCATGACGCTGGAAGTACGCTTCGCGCAcc  >  2:4383203/1‑150 (MQ=255)
                                                                                                                                                   ttttGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCAcc  >  2:6395674/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTTTTCCAGTTTTTTGATTTTTCAACAGCTGAATATCGGAAATAGTCAGTGACGAACATCAGTATTTTTACATAATTACGTCTAGGTCATTGTCGGTTTATCTCCTTTTACTTGTGGTAATTGAGGGATAACAACAAAAAAGGTGAGTTTTGCGACTCACCTTTTTTATTTGTTTCTTATTCAGAATCAATATCTTTTAAATCATCCTGAATCGCTTGTGCGTTCGGATTTTCCTGCGGTTTGAGTTCGCCGCCATTAGCGATGAAATCATGACGCTGGAAGTACGCTTCGCGCACC  >  NC_000913/2464076‑2464372

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: