Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,642 G→C 10.3% intergenic (‑48/+35) yfbU ← / ← yfbV UPF0304 family protein YfbU/UPF0208 membrane protein YfbV

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,412,6420GC10.3% 367.4 / 27.5 117intergenic (‑48/+35)yfbU/yfbVUPF0304 family protein YfbU/UPF0208 membrane protein YfbV
Reads supporting (aligned to +/- strand):  ref base G (41/63);  new base C (6/6);  total (48/69)
Fisher's exact test for biased strand distribution p-value = 5.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.88e-01

Only 100 of 118 total aligned reads displayed.

GCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCC  >  NC_000913/2412498‑2412791
                                                                                                                                                  |                                                                                                                                                     
gCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTCCTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGc                                                                                                                                                  >  2:1948352/1‑150 (MQ=255)
   ggCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCctt                                                                                                                                               >  1:1887208/1‑150 (MQ=255)
   ggCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCctt                                                                                                                                               >  1:65479/1‑150 (MQ=255)
     cGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCcttct                                                                                                                                             <  2:960776/150‑1 (MQ=255)
     cGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCcttct                                                                                                                                             >  1:1447744/1‑150 (MQ=255)
      gACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                                                            <  1:580872/150‑1 (MQ=255)
      gACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCcttctt                                                                                                                                            >  2:1670419/1‑150 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                                                          <  2:1680496/150‑1 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                                                          >  1:1165222/1‑150 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCATACTCAGCCTTCtttt                                                                                                                                          <  1:1179546/150‑1 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCATACTCAGCCTTCtttt                                                                                                                                          <  1:186297/150‑1 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCATACTCAGCCTTCtttt                                                                                                                                          <  1:935046/150‑1 (MQ=255)
          gTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCATACTCAGCCTTCtttttt                                                                                                                                        <  2:627747/150‑3 (MQ=255)
           tAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTATACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTtata                                                                                                                                       <  1:1116012/150‑1 (MQ=255)
              cGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGc                                                                                                                                    >  2:1386303/1‑150 (MQ=255)
              cGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGc                                                                                                                                    >  2:1379352/1‑150 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                   <  1:1996354/150‑1 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                   <  1:922609/150‑1 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                   <  1:530730/150‑1 (MQ=255)
                ttCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGt                                                                                                                                  >  1:1664730/1‑150 (MQ=255)
                 tCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTa                                                                                                                                 <  1:1435033/150‑1 (MQ=255)
                  cAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTaa                                                                                                                                <  2:1138710/150‑1 (MQ=255)
                    gCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATc                                                                                                                              >  1:1187662/1‑150 (MQ=255)
                    gCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATc                                                                                                                              >  1:217980/1‑150 (MQ=255)
                     cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                           <  1:1665393/136‑1 (MQ=255)
                     cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                           >  2:1665393/1‑136 (MQ=255)
                         tGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                         <  2:118204/150‑1 (MQ=255)
                         tGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                         <  1:78993/150‑1 (MQ=255)
                           ccGGTTCTAGCATATTAATCATCTTGTACTTATTAGACAAAATAAGACGTTTCGCTTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggggc                                                                                                                       <  1:697742/150‑3 (MQ=255)
                           ccGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTc                                                                                                                       <  1:1511741/150‑1 (MQ=255)
                              gATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAAt                                                                                                                    <  2:91674/150‑1 (MQ=255)
                                 cGAGCATAGTGATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTTCTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTAc                                                                                                                 <  1:1703372/150‑1 (MQ=255)
                                 cGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTAc                                                                                                                 >  1:1226511/1‑150 (MQ=255)
                                    gCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                                                              <  2:136298/150‑1 (MQ=255)
                                    gCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                                                              <  2:642275/150‑1 (MQ=255)
                                     cATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAat                                                                                                             <  2:1187662/150‑1 (MQ=255)
                                         gTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGAGGTCAATTACAAatcatc                                                                                                         <  1:796876/150‑1 (MQ=255)
                                          tcatcaTCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTatat                                                                                                                                      >  1:91950/1‑120 (MQ=255)
                                          tcatcaTCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTatat                                                                                                                                      <  2:91950/120‑1 (MQ=255)
                                           catcatCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCaa                                                                                                       <  2:1458230/150‑1 (MQ=255)
                                             tcatcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAgg                                                                                                     <  1:1085868/150‑1 (MQ=255)
                                               atcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGaa                                                                                                   >  1:509792/1‑150 (MQ=255)
                                                tcTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGaaa                                                                                                  <  2:684141/150‑1 (MQ=255)
                                                 cTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAg                                                                                                 <  2:1165222/150‑1 (MQ=255)
                                                   tGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAg                                                                                                 <  2:714098/148‑1 (MQ=255)
                                                   tGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAg                                                                                                 >  1:714098/1‑148 (MQ=255)
                                                    ttACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGttt                                                                                              <  2:538265/149‑1 (MQ=255)
                                                    gTACTGATTAGACATAATAAGACCTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTCCTCCGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACATATCATGAAGGAAAGttt                                                                                              <  2:888795/150‑1 (MQ=255)
                                                    gTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGttt                                                                                              <  2:684393/150‑1 (MQ=255)
                                                     tACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTAAATTACAAATCATCAAGGAAAGtttt                                                                                             >  2:224885/1‑150 (MQ=255)
                                                      aCTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                                                            <  2:1447744/150‑1 (MQ=255)
                                                      aCTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                                                            <  2:667148/150‑1 (MQ=255)
                                                       cTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTAt                                                                                           <  2:217980/150‑1 (MQ=255)
                                                       cTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTAt                                                                                           <  2:1180881/150‑1 (MQ=255)
                                                         gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                         <  2:1887208/150‑1 (MQ=255)
                                                         gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                         >  1:149061/1‑150 (MQ=255)
                                                              gACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCGTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttt                                                                                    >  1:1729659/1‑149 (MQ=255)
                                                               aCAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgt                                                                                   >  2:1129014/1‑150 (MQ=255)
                                                                cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                                  <  2:767554/150‑1 (MQ=255)
                                                                cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                                  >  2:767994/1‑150 (MQ=255)
                                                                          cGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                                                            >  1:505593/1‑130 (MQ=255)
                                                                          cGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGgcgcg                                                                        <  1:1670419/150‑1 (MQ=255)
                                                                           gTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGGCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGgcgcgc                                                                       <  2:562968/150‑1 (MQ=255)
                                                                                gcgTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTGAAAAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATACAGTTGTTTGAAGGCGCGCTTa                                                                    >  1:1331407/1‑150 (MQ=255)
                                                                                  gTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCg                                                                <  1:1157156/150‑1 (MQ=255)
                                                                                  gTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCg                                                                <  2:1226511/150‑1 (MQ=255)
                                                                                      gTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTc                                                            >  1:584446/1‑150 (MQ=255)
                                                                                      gTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTc                                                            >  1:1120992/1‑150 (MQ=255)
                                                                                        cATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                         <  1:1064732/87‑1 (MQ=255)
                                                                                        cATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                         >  2:1064732/1‑87 (MQ=255)
                                                                                          tttCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCt                                                        >  1:809991/1‑150 (MQ=255)
                                                                                               aTTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGc                                                   >  1:1898310/1‑150 (MQ=255)
                                                                                                  tttgACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTg                                                <  1:1948352/149‑1 (MQ=255)
                                                                                                  gttgACTCCTGTTTCACTCTACTACGGTG‑‑AAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTgg                                               <  1:550354/150‑1 (MQ=255)
                                                                                                  gttgACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTg                                                <  1:1944689/150‑1 (MQ=255)
                                                                                                  gttgACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTg                                                >  2:1551870/1‑150 (MQ=255)
                                                                                                            gTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc                                      >  2:1937384/1‑150 (MQ=255)
                                                                                                            gTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc                                      >  2:1344102/1‑150 (MQ=255)
                                                                                                               tCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCCGATTACCAGGcatt                                   <  2:1898310/150‑4 (MQ=255)
                                                                                                                cACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGc                                  >  2:931951/1‑150 (MQ=255)
                                                                                                                 aCTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGcc                                 <  1:1491462/150‑1 (MQ=255)
                                                                                                                  ctctACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCt                                <  1:114172/150‑1 (MQ=255)
                                                                                                                  ctctACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGATTACCAGGcattagc                                <  1:1386303/150‑7 (MQ=255)
                                                                                                                  ctctACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGATTACCAGGcattagc                                <  1:1379352/150‑7 (MQ=255)
                                                                                                                        tacGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACg                          <  2:418243/150‑1 (MQ=255)
                                                                                                                        tacGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACg                          <  1:1124645/150‑1 (MQ=255)
                                                                                                                            gTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTg                      <  2:1120992/150‑1 (MQ=255)
                                                                                                                                catataaaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCGATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc                    >  2:436681/6‑150 (MQ=255)
                                                                                                                                 aaaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGACAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCa                   >  1:1561939/1‑150 (MQ=255)
                                                                                                                                   aaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAac                 <  2:509792/150‑1 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCGATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc              <  1:436681/150‑1 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc                    >  2:254719/1‑144 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc                    <  1:254719/144‑1 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc              >  2:355649/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc              >  2:1844327/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc              >  1:116304/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc              >  1:1375876/1‑150 (MQ=255)
                                                                                                                                           agGCTGAGTATGCCTTCTTTTATATGCGTATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCa        >  1:1112180/1‑150 (MQ=255)
                                                                                                                                             gCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAg       >  1:1763838/1‑150 (MQ=255)
                                                                                                                                                  gTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTcc  <  1:910343/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCC  >  NC_000913/2412498‑2412791

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