Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,385,650 A→G 5.1% E47G (GAA→GGA)  ycjF → UPF0283 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,385,6500AG 5.1% 273.7 / 3.2 78E47G (GAA→GGA) ycjFUPF0283 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base G (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.22e-01

ATTGCGATGACCGAACCGTTAAAACCACGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTGGACAATGAATGCCT  >  NC_000913/1385505‑1385790
                                                                                                                                                 |                                                                                                                                            
aTTGCGATGACCGAACCGTTAAAACCACGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGgaag                                                                                                                                            <  1:1010114/148‑1 (MQ=255)
            gAACCGTTAAAACCACGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAg                                                                                                                                >  3:506132/1‑148 (MQ=255)
                 gTTAAAACCACGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGaa                                                                                                                           >  3:929823/1‑148 (MQ=255)
                       aCCACGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTa                                                                                                                     <  4:506132/148‑1 (MQ=255)
                           cGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATgg                                                                                                                 >  4:216063/1‑148 (MQ=255)
                           cGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATg                                                                                                                  >  2:710100/1‑147 (MQ=255)
                           cGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATg                                                                                                                  <  1:710100/147‑1 (MQ=255)
                            gTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGa                                                                                                                <  4:929823/148‑1 (MQ=255)
                               ttGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGc                                                                                                             >  4:1028141/1‑148 (MQ=255)
                                   tttCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAagaaga                                                                                                                                        <  4:1043183/117‑1 (MQ=255)
                                   tttCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAagaaga                                                                                                                                        >  3:1043183/1‑117 (MQ=255)
                                        aCGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGaa                                                                                                                           <  4:207381/125‑1 (MQ=255)
                                        aCGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGaa                                                                                                                           >  3:207381/1‑125 (MQ=255)
                                                ctGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAAc                                                                                            >  2:987786/1‑148 (MQ=255)
                                                    aGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAg                                                                                        <  3:1028141/148‑1 (MQ=255)
                                                     ggTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTa                                                                                                                     <  4:918174/118‑1 (MQ=255)
                                                     ggTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTa                                                                                                                     >  3:918174/1‑118 (MQ=255)
                                                       tCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCt                                                                                     >  3:226930/1‑148 (MQ=255)
                                                        cGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCtg                                                                                    <  1:987786/148‑1 (MQ=255)
                                                          aTCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATg                                                                                                                  <  1:744487/116‑1 (MQ=255)
                                                          aTCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATg                                                                                                                  >  2:744487/1‑116 (MQ=255)
                                                           tCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACg                                                                                           <  3:1096948/138‑1 (MQ=255)
                                                           tCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACg                                                                                           >  4:1096948/1‑138 (MQ=255)
                                                               aaTCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTc                                                                                                  <  1:5314/127‑1 (MQ=255)
                                                               aaTCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTc                                                                                                  >  2:5314/1‑127 (MQ=255)
                                                                    tAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATgg                                                                                                                 <  1:1016609/107‑1 (MQ=255)
                                                                    tAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATgg                                                                                                                 >  2:1016609/1‑107 (MQ=255)
                                                                    tAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATg                                                                        >  4:26334/1‑148 (MQ=255)
                                                                     aaaTTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGa                                                                                                                            >  3:688607/1‑95 (MQ=255)
                                                                     aaaTTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGa                                                                                                                            <  4:688607/95‑1 (MQ=255)
                                                                     aaaTTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGgcg                                                                               <  4:654828/140‑1 (MQ=255)
                                                                     aaaTTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGgcg                                                                               >  3:654828/1‑140 (MQ=255)
                                                                       aTTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCg                                                                                                         <  4:23626/112‑1 (MQ=255)
                                                                       aTTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCg                                                                                                         >  3:23626/1‑112 (MQ=255)
                                                                          cAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATg                                                                  <  3:216063/148‑1 (MQ=255)
                                                                            gggCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGaagaag                                                                                                                                         <  4:675834/75‑1 (MQ=255)
                                                                            gggCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGaagaag                                                                                                                                         >  3:675834/1‑75 (MQ=255)
                                                                            gggCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATggg                                                                <  3:26334/148‑1 (MQ=255)
                                                                                      gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAggg                                                                                                                                     <  2:446836/69‑1 (MQ=255)
                                                                                      gACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAggg                                                                                                                                     >  1:446836/1‑69 (MQ=255)
                                                                                       aCCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATgg                                                                       >  4:833712/1‑130 (MQ=255)
                                                                                       aCCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATgg                                                                       <  3:833712/130‑1 (MQ=255)
                                                                                           ttGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGc                                                                                                             <  3:167516/88‑1 (MQ=255)
                                                                                           ttGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGc                                                                                                             >  4:167516/1‑88 (MQ=255)
                                                                                              aCGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTgg                                              <  4:226930/148‑1 (MQ=255)
                                                                                               cGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGgc                                             >  2:1092637/1‑148 (MQ=255)
                                                                                                 aaaaTCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGt                                                                      >  3:420851/1‑121 (MQ=255)
                                                                                                 aaaaTCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGt                                                                      <  4:420851/121‑1 (MQ=255)
                                                                                                 aaaaTCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGgcgc                                           >  4:19680/1‑148 (MQ=255)
                                                                                                    aTCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAg                                        <  1:1092637/148‑1 (MQ=255)
                                                                                                      cAGGTGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGttt                                                <  2:707752/138‑1 (MQ=255)
                                                                                                      cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCg                                                                                                         >  2:199781/1‑81 (MQ=255)
                                                                                                      cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCg                                                                                                         <  1:199781/81‑1 (MQ=255)
                                                                                                      cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGttt                                                >  1:707752/1‑138 (MQ=255)
                                                                                                      cAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCg                                      >  1:52974/1‑148 (MQ=255)
                                                                                                         gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGGAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGctc                                                     <  4:603173/130‑1 (MQ=255)
                                                                                                         gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGGAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGctc                                                     >  1:375885/1‑130 (MQ=255)
                                                                                                         gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGGAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGctc                                                     <  2:375885/130‑1 (MQ=255)
                                                                                                         gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGGAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGctc                                                     >  3:603173/1‑130 (MQ=255)
                                                                                                         gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTgg                                              <  2:945323/137‑1 (MQ=255)
                                                                                                         gcgcAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTgg                                              >  1:945323/1‑137 (MQ=255)
                                                                                                         gcgcAAAATTTTGCACCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTg                                   >  4:1093071/1‑148 (MQ=255)
                                                                                                            cAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCg                                <  3:1093071/148‑1 (MQ=255)
                                                                                                              aaaTTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCggg                              >  1:118831/1‑148 (MQ=255)
                                                                                                                aTTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAg                                        >  1:60561/1‑136 (MQ=255)
                                                                                                                aTTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAg                                        <  2:60561/136‑1 (MQ=255)
                                                                                                                 ttttGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATg                                                                  <  2:313769/109‑1 (MQ=255)
                                                                                                                 ttttGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATg                                                                  >  1:313769/1‑109 (MQ=255)
                                                                                                                            ccACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTg                >  1:68358/1‑148 (MQ=255)
                                                                                                                             cACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTgg               <  2:52974/148‑1 (MQ=255)
                                                                                                                             cACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTgg               <  3:19680/148‑1 (MQ=255)
                                                                                                                                 cTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAAt                                                                         <  4:1133579/86‑1 (MQ=255)
                                                                                                                                 cTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAAt                                                                         >  3:1133579/1‑86 (MQ=255)
                                                                                                                                     acgaAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAg                  <  1:303135/137‑1 (MQ=255)
                                                                                                                                     acgaAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAg                  >  2:303135/1‑137 (MQ=255)
                                                                                                                                      cgaAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGt                                  >  3:69870/1‑120 (MQ=255)
                                                                                                                                      cgaAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGt                                  <  4:69870/120‑1 (MQ=255)
                                                                                                                                       gaAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTGGACaatgaatg     >  3:680999/1‑148 (MQ=255)
                                                                                                                                        aaGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTGGACaat          >  1:206708/1‑142 (MQ=255)
                                                                                                                                        aaGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTGGACaat          <  2:206708/142‑1 (MQ=255)
                                                                                                                                          gcgcAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTGGACAATGAATGCCt  <  4:680999/148‑1 (MQ=255)
                                                                                                                                                 aagaagAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAg                  <  3:814702/125‑1 (MQ=255)
                                                                                                                                                 aagaagAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAg                  >  4:814702/1‑125 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
ATTGCGATGACCGAACCGTTAAAACCACGTATTGATTTCGACGGTCCTCTGGAGGTCGATCAGAATCCTAAATTCAGGGCGCAGCAGACCTTTGACGAAAATCAGGCGCAAAATTTTGCCCCGGCCACGCTCGACGAAGCGCAGGAAGAAGAGGGGCAAGTCGAAGCGGTAATGGACGCAGCGTTACGTCCGAAACGCAGCCTGTGGCGCAAAATGGTGATGGGAGGGCTGGCTCTGTTTGGCGCAAGCGTTGTCGGGCAGGGTGTACAGTGGACAATGAATGCCT  >  NC_000913/1385505‑1385790

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