Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 672,113 T→C 5.3% D25G (GAT→GGT)  lptE ← LPS assembly OM complex LptDE, lipoprotein component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913672,1130TC 5.3% 784.5 / 21.7 226D25G (GAT→GGT) lptELPS assembly OM complex LptDE, lipoprotein component
Reads supporting (aligned to +/- strand):  ref base T (104/110);  new base C (4/8);  total (108/118)
Fisher's exact test for biased strand distribution p-value = 3.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 227 total aligned reads displayed.

GATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACG  >  NC_000913/671966‑672261
                                                                                                                                                   |                                                                                                                                                    
gATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATc                                                                                                                                                     <  2:94443/149‑1 (MQ=255)
  tGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCAc                                                                                                                                                   <  4:899345/149‑1 (MQ=255)
     aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                                >  3:724078/1‑149 (MQ=255)
     aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                                >  4:448236/1‑149 (MQ=255)
     aaCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                                <  3:74145/149‑1 (MQ=255)
      aCGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                               <  2:409834/149‑1 (MQ=255)
        gTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGAt                                                                                                                                             >  1:836300/1‑149 (MQ=255)
          ccTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAATCTGCGTGATACCACGCAggttc                                                                                                                                           <  3:1438493/149‑5 (MQ=255)
                gcgTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                     <  4:132549/149‑1 (MQ=255)
                  gTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                                >  2:292541/1‑136 (MQ=255)
                  gTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCa                                                                                                                                                <  1:292541/136‑1 (MQ=255)
                    ggTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccagcca                                                                                                                                     <  2:581380/145‑1 (MQ=255)
                    ggTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACAcc                                                                                                                                 >  2:529293/1‑149 (MQ=255)
                          tttATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcgg                                                                                                                           >  1:277925/1‑149 (MQ=255)
                            tATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                         >  3:1295920/1‑149 (MQ=255)
                              tCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcg                                                                                                                            <  4:249250/144‑1 (MQ=255)
                               cAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                      <  4:1540635/149‑1 (MQ=255)
                               cAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGAtt                                                                                                                      <  2:1414234/149‑1 (MQ=255)
                                 aGCAACTTGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                    <  3:1790760/149‑1 (MQ=255)
                                 aGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTaa                                                                                                                    <  4:724078/149‑1 (MQ=255)
                                  gCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                         <  4:1107595/143‑1 (MQ=255)
                                  gCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                         >  3:1107595/1‑143 (MQ=255)
                                   cAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAaca                                                                                                                  <  4:1295920/149‑1 (MQ=255)
                                       tCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGc                                                                                                              >  2:452670/1‑149 (MQ=255)
                                         gACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACCCCCGGCGGTGATTAACACCGCCa                                                                                                            >  4:813665/1‑149 (MQ=255)
                                         gACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCa                                                                                                            >  1:1259460/1‑149 (MQ=255)
                                          acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcg                                                                                                                            <  2:495295/132‑1 (MQ=255)
                                          acacCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcg                                                                                                                            >  1:495295/1‑132 (MQ=255)
                                             ccATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAgaga                                                                                                        >  1:1374124/1‑149 (MQ=255)
                                              cATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGAt                                                                                                       >  4:308176/1‑149 (MQ=255)
                                                ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCcggcg                                                                                                                            >  2:923828/1‑126 (MQ=255)
                                                ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAaca                                                                                                                  >  1:1471394/1‑136 (MQ=255)
                                                ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaa                                                                                                     >  1:1594703/1‑149 (MQ=255)
                                                ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                               >  2:1193948/1‑107 (MQ=255)
                                                ttCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGATACCACGCAg                                                                                                                                               <  1:1193948/107‑1 (MQ=255)
                                                   aGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaa                                                                                                  <  3:448236/149‑1 (MQ=255)
                                                    gACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaac                                                                                                 <  3:1539309/149‑1 (MQ=255)
                                                          aCTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGtt                                                                                           <  1:414905/149‑1 (MQ=255)
                                                            tGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGc                                                                                         <  2:277925/149‑1 (MQ=255)
                                                               ttACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                        <  1:1145096/115‑1 (MQ=255)
                                                               ttACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGa                                                                                                                        >  2:1145096/1‑115 (MQ=255)
                                                                   gcaccgcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATgccag                                                                                                                                        <  4:1031689/95‑1 (MQ=255)
                                                                        gcacGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCAc                                                                             >  1:1512820/1‑149 (MQ=255)
                                                                             gCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                        <  1:529293/149‑1 (MQ=255)
                                                                                tAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCag                                                                                                           >  4:299224/1‑111 (MQ=255)
                                                                                tAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtccta                                                                     <  2:1594703/149‑2 (MQ=255)
                                                                                       ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTg                                                                                                                         <  3:1674278/90‑1 (MQ=255)
                                                                                       ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATaacaac                                                                                                 >  4:1881977/1‑114 (MQ=255)
                                                                                       ccGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtc                                                                        <  3:898596/139‑1 (MQ=255)
                                                                                         gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                            <  2:1796777/149‑1 (MQ=255)
                                                                                         gTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTa                                                            >  2:1551582/1‑149 (MQ=255)
                                                                                           tCGGATCGCCTGAGTCCAGGATCAGGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATCTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGc                                                          <  2:1214621/149‑1 (MQ=255)
                                                                                              gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                              >  4:1538652/1‑110 (MQ=255)
                                                                                              gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                                 <  3:1360721/123‑1 (MQ=255)
                                                                                              gATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATc                                                                                 <  1:58914/123‑1 (MQ=255)
                                                                                                tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGTTAATCACCGccgg                                                                                                           <  3:1860847/95‑3 (MQ=255)
                                                                                                tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAAt                                                                                              >  3:610104/1‑108 (MQ=255)
                                                                                                tCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACg                                                     >  1:143945/1‑149 (MQ=255)
                                                                                                  gCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCa                                                                                       <  1:9256/113‑1 (MQ=255)
                                                                                                    cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcctcc                                                                    >  2:200629/1‑130 (MQ=255)
                                                                                                    cTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcctcc                                                                    <  1:200629/130‑1 (MQ=255)
                                                                                                          ccAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTga                                           <  1:1551582/149‑1 (MQ=255)
                                                                                                           cAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          >  4:1409851/1‑149 (MQ=255)
                                                                                                            aGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagg                                         >  4:569368/1‑149 (MQ=255)
                                                                                                             ggATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                        <  2:143945/149‑1 (MQ=255)
                                                                                                              gATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgaggag                                       <  3:30453/149‑1 (MQ=255)
                                                                                                              gATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          >  2:130951/1‑146 (MQ=255)
                                                                                                                tCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                     <  3:1059704/149‑1 (MQ=255)
                                                                                                                tCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                     >  1:1327593/1‑149 (MQ=255)
                                                                                                                  aTGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          <  1:980607/142‑1 (MQ=255)
                                                                                                                  aTGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          >  2:980607/1‑142 (MQ=255)
                                                                                                                  aTGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGtt                                     <  3:254192/147‑1 (MQ=255)
                                                                                                                  aTGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTa                                   <  2:1512820/149‑1 (MQ=255)
                                                                                                                   tGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAtt                                             <  2:1708262/138‑1 (MQ=255)
                                                                                                                    gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCag                                                                                                           >  1:678345/1‑75 (MQ=255)
                                                                                                                    gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTa                                   >  4:622736/1‑147 (MQ=255)
                                                                                                                    gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                                 <  2:135977/149‑1 (MQ=255)
                                                                                                                    gACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTAcc                                 >  3:792076/1‑149 (MQ=255)
                                                                                                                       cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCa                                                                              >  3:494685/1‑101 (MQ=255)
                                                                                                                       cTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                              >  1:1778414/1‑149 (MQ=255)
                                                                                                                        ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGata                                                                                   <  1:184691/95‑1 (MQ=255)
                                                                                                                        ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGAt                                              >  4:589523/1‑132 (MQ=255)
                                                                                                                        ttCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGttt                                    >  2:651853/1‑142 (MQ=255)
                                                                                                                          cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcct                                                                      >  3:1301802/1‑106 (MQ=255)
                                                                                                                          cATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTtcct                                                                      <  4:1301802/106‑1 (MQ=255)
                                                                                                                                ggaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAAt                     >  2:756937/1‑149 (MQ=255)
                                                                                                                                 gaaggaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATc                    >  1:45518/1‑149 (MQ=255)
                                                                                                                                     gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCt                                >  3:220212/1‑133 (MQ=255)
                                                                                                                                     gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                <  4:1684646/149‑1 (MQ=255)
                                                                                                                                     gaaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACtt                >  1:1036550/1‑149 (MQ=255)
                                                                                                                                      aaCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAg                                                           <  3:1600201/105‑1 (MQ=255)
                                                                                                                                       atcTGCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgagga                                        <  3:317210/121‑1 (MQ=255)
                                                                                                                                       aCCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTgg                              >  4:162831/1‑133 (MQ=255)
                                                                                                                                           gCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTa          <  4:792076/149‑1 (MQ=255)
                                                                                                                                           gCGTGATACCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCAc                                                                             >  1:1684223/1‑82 (MQ=255)
                                                                                                                                            cGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTgag                                          <  2:147980/116‑1 (MQ=255)
                                                                                                                                                 tATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTaa    >  4:33418/1‑149 (MQ=255)
                                                                                                                                                  aTCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc   <  1:374372/149‑1 (MQ=255)
                                                                                                                                                  aTCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAAc   >  2:374372/1‑149 (MQ=255)
                                                                                                                                                   tCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACg  >  4:1800140/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATGGAACGTCCTTACGCGTGGTTTCTTTATCAAGCAACTCGACACCATTCAGACGTAACTGGTTACGCACCGCACGGCTTAATGGCCCGTTCGGATCGCCTGAGTCCAGGATCATGACCTTCATAGTGGAAGGAACCTGCGTGGTATCACGCAGATGCCAGCCACACCCGGCGGTGATTAACACCGCCAGAGATAACAACAATGTTGCCAGATATCGCACGCTTCCTCCCGCGCTTAGCCAACGACCAGATTGAGGAGTTTACCTGGTACGTAAATCACTTTACGTACAGTAACG  >  NC_000913/671966‑672261

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