Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,639 T→A 10.4% intergenic (‑45/+38) 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,6390TA10.4% 319.5 / 23.7 116intergenic (‑45/+38)yfbU/yfbVUPF0304 family protein YfbU/UPF0208 membrane protein YfbV
Reads supporting (aligned to +/- strand):  ref base T (41/63);  new base A (6/6);  total (47/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 = 2.10e-01

Only 100 of 117 total aligned reads displayed.

GCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAG  >  NC_000913/2412498‑2412786
                                                                                                                                               |                                                                                                                                                   
   ggCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCctt                                                                                                                                          >  1:65479/1‑150 (MQ=255)
   ggCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCctt                                                                                                                                          >  1:1887208/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)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                                                     <  2:1680496/150‑1 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCtttt                                                                                                                                     >  1:1165222/1‑150 (MQ=255)
        cGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCATACTCAGCCTTCtttt                                                                                                                                     <  1:1179546/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)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                              <  1:1996354/150‑1 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                              >  1:1988364/1‑150 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                              <  1:530730/150‑1 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                              <  1:922609/150‑1 (MQ=255)
                  cAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTaa                                                                                                                           <  2:1138710/150‑1 (MQ=255)
                    gCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATc                                                                                                                         >  1:1187662/1‑150 (MQ=255)
                     cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                      >  2:1665393/1‑136 (MQ=255)
                     cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                      <  1:1665393/136‑1 (MQ=255)
                         tGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                    <  2:118204/150‑1 (MQ=255)
                          gCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGt                                                                                                                   >  1:684141/1‑150 (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:642275/150‑1 (MQ=255)
                                    gCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                                                         <  1:174060/150‑1 (MQ=255)
                                    gCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACaaa                                                                                                         <  2:136298/150‑1 (MQ=255)
                                     cATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAat                                                                                                        <  2:1187662/150‑1 (MQ=255)
                                         gTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAatcatc                                                                                                    >  1:1393124/1‑150 (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:1393124/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)
                                                         gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                    <  2:1887208/150‑1 (MQ=255)
                                                         gATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATcc                                                                                    >  1:149061/1‑150 (MQ=255)
                                                           ttAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAg                                                                                  >  2:1124645/1‑150 (MQ=255)
                                                              gACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCGTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttt                                                                               >  1:1729659/1‑149 (MQ=255)
                                                               aCAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgt                                                                              >  2:1129014/1‑150 (MQ=255)
                                                                cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                             <  1:1644166/150‑1 (MQ=255)
                                                                cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                             <  2:767554/150‑1 (MQ=255)
                                                                cAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAgttgtt                                                                             >  2:767994/1‑150 (MQ=255)
                                                                          cGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                                                       <  2:505593/130‑1 (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                                                           >  2:1205062/1‑150 (MQ=255)
                                                                                  gTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCg                                                           <  2:1226511/150‑1 (MQ=255)
                                                                                  gTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCg                                                           <  1:1157156/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)
                                                                                                  gttgACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTg                                           >  2:1551870/1‑150 (MQ=255)
                                                                                                  gttgACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTg                                           <  1:1944689/150‑1 (MQ=255)
                                                                                                            gTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc                                 >  2:1344102/1‑150 (MQ=255)
                                                                                                            gTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGc                                 >  2:1937384/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)
                                                                                                                      aatacGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGGGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCaa                       <  1:883329/148‑1 (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)
                                                                                                                             tG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGc                <  2:1988364/150‑1 (MQ=255)
                                                                                                                                 aaaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGACAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCa              >  1:1561939/1‑150 (MQ=255)
                                                                                                                                   aaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAac            <  2:509792/150‑1 (MQ=255)
                                                                                                                                      aaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc         >  1:507155/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCGATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc         <  1:436681/150‑1 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc               >  2:254719/1‑144 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc               <  1:254719/144‑1 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc         >  1:1375876/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc         >  2:355649/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc         >  2:1844327/1‑150 (MQ=255)
                                                                                                                                      aaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACAcc         >  1:116304/1‑150 (MQ=255)
                                                                                                                                           agGCTGAGTATGCCTTCTTTTATATGCGTATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCa   >  1:1112180/1‑150 (MQ=255)
                                                                                                                                            cgatGATTATGCCTTCTTTTATAAGAGTAATCAGGGATCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCa   >  2:1767831/4‑150 (MQ=255)
                                                                                                                                             gCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAg  >  1:1763838/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAG  >  NC_000913/2412498‑2412786

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: