Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,412,645 T→A 9.7% intergenic (‑51/+32) 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,6450TA 9.7% 353.2 / 21.2 120intergenic (‑51/+32)yfbU/yfbVUPF0304 family protein YfbU/UPF0208 membrane protein YfbV
Reads supporting (aligned to +/- strand):  ref base T (43/65);  new base A (6/6);  total (49/71)
Fisher's exact test for biased strand distribution p-value = 5.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 120 total aligned reads displayed.

GCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTG  >  NC_000913/2412498‑2412793
                                                                                                                                                     |                                                                                                                                                    
gCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTCCTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGc                                                                                                                                                    >  2:1948352/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: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:922609/150‑1 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                     <  1:1996354/150‑1 (MQ=255)
               gTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCg                                                                                                                                     >  1:1988364/1‑150 (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)
                     cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                             >  2:1665393/1‑136 (MQ=255)
                     cATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCttt                                                                                                                                             <  1:1665393/136‑1 (MQ=255)
                         tGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                           <  1:78993/150‑1 (MQ=255)
                          gCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGt                                                                                                                          >  1:684141/1‑150 (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)
                                     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                                                                                                   >  1:714098/1‑148 (MQ=255)
                                                   tGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAg                                                                                                   <  2:714098/148‑1 (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:65479/150‑1 (MQ=255)
                                                    gTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGttt                                                                                                <  2:684393/150‑1 (MQ=255)
                                                      aCTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTa                                                                                              <  2:1447744/150‑1 (MQ=255)
                                                       cTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTAt                                                                                             <  2:217980/150‑1 (MQ=255)
                                                       cTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTAt                                                                                             <  2:1180881/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)
                                                                           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:1120992/1‑150 (MQ=255)
                                                                                      gTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTc                                                              >  1:584446/1‑150 (MQ=255)
                                                                                        cATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                           <  1:1064732/87‑1 (MQ=255)
                                                                                        cATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAgggg                                                                                                                           >  2:1064732/1‑87 (MQ=255)
                                                                                         aTTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGc                                                           <  1:1205062/150‑1 (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)
                                                                                                            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)
                                                                                                                  ctctACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCt                                  <  1:114172/150‑1 (MQ=255)
                                                                                                                  ctctACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGATTACCAGGcattagc                                  <  1:1386303/150‑7 (MQ=255)
                                                                                                                        tacGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACg                            <  2:418243/150‑1 (MQ=255)
                                                                                                                             tG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGc                       <  2:1988364/150‑1 (MQ=255)
                                                                                                                                catataaaaGAAGGCATACTCAGCCTTCTTTTATATGCGTAATCAGGGGTCGATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGcc                      >  2:436681/6‑150 (MQ=255)
                                                                                                                                 aaaaaaaaGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGACAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCa                     >  1:1561939/1‑150 (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                >  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)
                                                                                                                                          aagGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTc            <  1:1129014/150‑1 (MQ=255)
                                                                                                                                           agGCTGAGTATGCCTTCTTTTATATGCGTATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCa          >  1:1112180/1‑150 (MQ=255)
                                                                                                                                            cgatGATTATGCCTTCTTTTATAAGAGTAATCAGGGATCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCa          >  2:1767831/4‑150 (MQ=255)
                                                                                                                                             gCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAg         >  1:1763838/1‑150 (MQ=255)
                                                                                                                                                  gTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTcc    <  1:910343/150‑1 (MQ=255)
                                                                                                                                                   tATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCt   >  2:1314674/1‑150 (MQ=255)
                                                                                                                                                    aTGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCATTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTg  <  1:872439/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCAGGCGACGGTAACGTTCAGCATTTGCCGGATCGAGCATAGTCATCATCTTGTACTGATTAGACAAAATAAGACGTTGCGCGTTGGTCATTTCCATTGTTGACTCCTGTATCACTCTACTACGGTG‑‑AAAAAAAAGAAGGCTGAGTATGCCTTCTTTTATATGCGTAATCAGGGGTCAATTACAAATCATCAAGGAAAGTTTTATCCAGTTGTTTGAAGGCGCGCTTAAGCGTGTCAGCTAATGCCTGGTAATCAGGCTTGCCTTCAACGGGTGCCAACACCTGTCCAGACTCCTG  >  NC_000913/2412498‑2412793

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: