Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 344,704 A→T 5.9% K177* (AAA→TAA) ‡ yahL → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913344,7040AT 5.9% 429.4 / 12.1 118K177* (AAA→TAA) ‡yahLuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (51/60);  new base T (3/4);  total (54/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 122 total aligned reads displayed.

ATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTG  >  NC_000913/344558‑344850
                                                                                                                                                  |                                                                                                                                                  
atgatgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                                                                  <  3:848514/149‑1 (MQ=255)
atgatgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                                                                  <  1:1489024/149‑1 (MQ=255)
 tgatgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAg                                                                                                                                                 <  1:1203802/149‑1 (MQ=255)
 tgatgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAg                                                                                                                                                 <  1:77309/149‑1 (MQ=255)
  gatgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAgg                                                                                                                                                <  1:333432/149‑1 (MQ=255)
    tgCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGa                                                                                                                                              <  2:732852/149‑1 (MQ=255)
      cGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                            <  4:120721/149‑1 (MQ=255)
       gTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTTTAGGGaaaa                                                                                                                                           <  1:231169/149‑2 (MQ=255)
        tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                                                                  >  2:1821107/1‑141 (MQ=255)
        tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTaaa                                                                                                                                                  <  1:1821107/141‑1 (MQ=255)
        tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                                                                          <  1:1602107/149‑1 (MQ=255)
        tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                                                                           <  4:823058/148‑1 (MQ=255)
        tttCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAag                                                                                                                                           >  3:823058/1‑148 (MQ=255)
            ttGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgtgt                                                                                                                                      >  4:1858303/1‑149 (MQ=255)
               atatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTAt                                                                                                                                    <  2:1285328/148‑1 (MQ=255)
               atatTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTAt                                                                                                                                    >  1:1285328/1‑148 (MQ=255)
                 atTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAgg                                                                                                                                                <  1:569734/134‑1 (MQ=255)
                 atTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAgg                                                                                                                                                >  2:569734/1‑134 (MQ=255)
                  ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                            >  1:426022/1‑137 (MQ=255)
                  ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                            <  2:426022/137‑1 (MQ=255)
                  ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                                                                          <  3:527441/139‑1 (MQ=255)
                  ttAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAaga                                                                                                                                          >  4:527441/1‑139 (MQ=255)
                     aaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgtg                                                                                                                             <  2:653857/149‑1 (MQ=255)
                     aaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgtg                                                                                                                             >  2:434298/1‑149 (MQ=255)
                     aaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                                                              >  1:1797006/1‑148 (MQ=255)
                     aaGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                                                              <  2:1797006/148‑1 (MQ=255)
                       gCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgt                                                                                                                                        <  2:988272/136‑1 (MQ=255)
                       gCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgt                                                                                                                                        >  1:988272/1‑136 (MQ=255)
                         aaaCGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGaaa                                                                                                                         >  2:1816874/1‑149 (MQ=255)
                            cGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAAt                                                                                                                      <  3:367784/149‑1 (MQ=255)
                             gTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATg                                                                                                                     >  4:899521/1‑149 (MQ=255)
                                  ggCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                            >  2:892751/1‑121 (MQ=255)
                                  ggCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGaaa                                                                                                                                            <  1:892751/121‑1 (MQ=255)
                                   gCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgtgt                                                                                                                                      <  3:773250/126‑1 (MQ=255)
                                   gCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAgtgt                                                                                                                                      >  4:773250/1‑126 (MQ=255)
                                     tGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTcgcg                                                                                                             >  2:484322/1‑149 (MQ=255)
                                       gACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaa                                                                                                           >  3:1885501/1‑149 (MQ=255)
                                        aCGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaaa                                                                                                          >  2:1226183/1‑149 (MQ=255)
                                             ttACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcg                                                                                                     <  3:899521/149‑1 (MQ=255)
                                               aCCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgtc                                                                                                   <  3:313282/149‑1 (MQ=255)
                                                      ttAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                                                            <  1:853649/149‑1 (MQ=255)
                                                      ttAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCg                                                                                            >  4:329398/1‑149 (MQ=255)
                                                        aaGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                          <  2:298235/149‑1 (MQ=255)
                                                        aaGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                          <  1:484322/149‑1 (MQ=255)
                                                        aaGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                          <  4:199920/149‑1 (MQ=255)
                                                          gAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGgaga                                                                                        <  1:1226183/149‑1 (MQ=255)
                                                           aaaTTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagac                                                                                       >  2:1402125/1‑149 (MQ=255)
                                                            aaTTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtc                                                                                                      >  4:1454341/1‑133 (MQ=255)
                                                            aaTTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtc                                                                                                      <  3:1454341/133‑1 (MQ=255)
                                                             aTTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacag                                                                                     >  3:1413995/1‑149 (MQ=255)
                                                              ttAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                                                              >  2:621391/1‑107 (MQ=255)
                                                              ttAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCAtgt                                                                                                                              <  1:621391/107‑1 (MQ=255)
                                                               tAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                                   <  3:1858303/149‑1 (MQ=255)
                                                                    cAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaaa                                                                                                          <  1:137260/121‑1 (MQ=255)
                                                                    cAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGaaa                                                                                                          >  2:137260/1‑121 (MQ=255)
                                                                           tCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                                   >  2:785887/1‑137 (MQ=255)
                                                                           tCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                                   <  1:785887/137‑1 (MQ=255)
                                                                             gATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCa                                                                     <  3:329398/149‑1 (MQ=255)
                                                                               ttCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAgg                                                                   <  4:1885501/149‑1 (MQ=255)
                                                                                tCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgc                                                                  <  1:434298/149‑1 (MQ=255)
                                                                                tCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgc                                                                  <  1:22245/149‑1 (MQ=255)
                                                                                   gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                                                                         <  1:9952/107‑1 (MQ=255)
                                                                                   gTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAAt                                                                                                         >  2:9952/1‑107 (MQ=255)
                                                                                          cAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCAt                                                        <  4:1413995/149‑1 (MQ=255)
                                                                                                  aaaTGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCACTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgcg                                                                 <  3:1367437/132‑1 (MQ=255)
                                                                                                        aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgtc                                                                                                   >  4:741152/1‑92 (MQ=255)
                                                                                                        aTTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATgtcgtc                                                                                                   <  3:741152/92‑1 (MQ=255)
                                                                                                          tcttctTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTAtt                                                                                               >  4:1406036/1‑94 (MQ=255)
                                                                                                            ttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                >  4:374754/1‑139 (MQ=255)
                                                                                                            ttcttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                <  3:374754/139‑1 (MQ=255)
                                                                                                               ttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCAt                                                                           >  4:340710/1‑109 (MQ=255)
                                                                                                               ttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                <  1:223792/136‑1 (MQ=255)
                                                                                                               ttGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                >  2:223792/1‑136 (MQ=255)
                                                                                                                       cATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGgc                                                                  <  1:1778173/110‑1 (MQ=255)
                                                                                                                       cATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTa                           >  2:312233/1‑149 (MQ=255)
                                                                                                                               tttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGttt                   >  4:1115653/1‑149 (MQ=255)
                                                                                                                               tttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGttt                   >  4:123894/1‑149 (MQ=255)
                                                                                                                               tttGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGAGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                >  2:1530010/1‑120 (MQ=255)
                                                                                                                               tttGCTGAGAATTGTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGACATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                <  1:1530010/120‑1 (MQ=255)
                                                                                                                                   cTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTcc               <  3:1115653/149‑1 (MQ=255)
                                                                                                                                    tGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                          >  4:1579919/1‑73 (MQ=255)
                                                                                                                                    tGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTaa                                                                        >  2:73399/1‑91 (MQ=255)
                                                                                                                                    tGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTAc                          >  4:1707275/1‑137 (MQ=255)
                                                                                                                                    tGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTAc                          <  3:1707275/137‑1 (MQ=255)
                                                                                                                                    tGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGtt                    >  4:1209509/1‑143 (MQ=255)
                                                                                                                                    tGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGtt                    <  3:1209509/143‑1 (MQ=255)
                                                                                                                                    tGAGAATCTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGga                                                                                          <  3:1579919/73‑1 (MQ=255)
                                                                                                                                     gagaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGa                                                     >  1:558803/1‑109 (MQ=255)
                                                                                                                                     gagaATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGGCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGa                                                     <  2:558803/109‑1 (MQ=255)
                                                                                                                                         tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatta                                               <  3:1795680/109‑1 (MQ=255)
                                                                                                                                         tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatta                                               >  4:1795680/3‑111 (MQ=255)
                                                                                                                                         tctttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattatt                                                <  1:689963/108‑1 (MQ=255)
                                                                                                                                            ttCCCTTTAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACt      <  1:1167862/149‑1 (MQ=255)
                                                                                                                                              cccTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTgg    <  1:312233/149‑1 (MQ=255)
                                                                                                                                                cTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTg  >  2:1149469/1‑149 (MQ=255)
                                                                                                                                                cTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTg  >  2:1098143/1‑149 (MQ=255)
                                                                                                                                                 tAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                                   >  4:517208/1‑67 (MQ=255)
                                                                                                                                                 tAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGagacagac                                                                                   <  3:517208/67‑1 (MQ=255)
                                                                                                                                                 tAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattat                                                 <  1:914734/101‑1 (MQ=255)
                                                                                                                                                 tAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGattat                                                 >  2:914734/1‑101 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATGATGCGTTTCTTGATATTAAAGCAAACGTGCAGGCTGGACGCTTTACCGGGCTTAAGAAATTAGTGCAGTTTATCGATTCAGTGCGAACAGATTGCAAATGTATTCTTCTTGAATACCATATGCCTTTGCTGAGAATTTTCCCTAAAGGGAAAGAGTGTATGCATGTGGAAAAATGGCTTCGCGAAATGTCGTCTATTCCGGAGACAGACGCGCCCATTAAGCAGGCGCTTGCGCATGGATTATTATTACACTTAAAAAATATCTACCCTGTTTTTCCTGAATCACTGGTG  >  NC_000913/344558‑344850

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: