Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,584,541 G→A 7.1% I397I (ATC→ATT) ‡ parC ← DNA topoisomerase IV subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,584,5410GA 7.1% 207.5 / 3.8 72I397I (ATC→ATT) ‡parCDNA topoisomerase IV subunit A
Reads supporting (aligned to +/- strand):  ref base G (24/41);  new base A (2/3);  total (26/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

ACTCATCCTGCTCACCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCACTCCACCAATAGGGTTCGCAGGTTT  >  NC_002947/5584393‑5584688
                                                                                                                                                       |                                                                                                                                                   
aCTCATCCTGCTCACCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTgg                                                                                                                                                     <  8:143762/149‑1 (MQ=255)
     tCCTGCTCACCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCAc                                                                                                                                                <  8:34521/149‑1 (MQ=255)
       cTGCTCACCGCGGATTTTCTTTTCTTCCAGTCGCGCCAATTGCCGATGGTGGGTTTTGTGTATTTAGCCTCNGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCAc                                                                                                                                                <  1:131064/149‑1 (MQ=255)
       cTGCTCACCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGGCGATATTCCTCGGTACGGATGATGTGGATCACt                                                                                                                                               <  1:257192/149‑1 (MQ=255)
            cACCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑GCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGAAGTGATCCACATCATcc                                                                                                                                         <  1:152285/149‑1 (MQ=255)
             aCCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGACAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATc                                                                                                                                          >  1:75901/1‑147 (MQ=255)
              ccGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATACTCCTCGGTACGGATGAAGTGATCCACATCATccgt                                                                                                                                       <  8:2798/149‑3 (MQ=255)
                     tCTTCATCTCTTCCAGGCGCGCAAACTGCCGGAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAg                                                                                                                                 <  8:278547/148‑1 (MQ=255)
                       ttCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTc                                                                                                                                             <  6:157867/134‑1 (MQ=255)
                       ttCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTc                                                                                                                                             >  5:157867/1‑134 (MQ=255)
                        tCATCTCTTCCAGGCGCGCCAACTCCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTAAGGATGATGTGGAGCACTTCATCCAGGTGGAGGAAt                                                                                                                             >  8:309963/1‑149 (MQ=255)
                         cATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAgg                                                                                                                                <  4:154770/145‑1 (MQ=255)
                         cATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAg                                                                                                                                 <  3:58780/144‑1 (MQ=255)
                         cATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATg                                                                                                                            >  3:387480/1‑149 (MQ=255)
                         cATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑GCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAAt                                                                                                                             <  8:261943/148‑1 (MQ=255)
                          atatctTCCAGGCGCGCCAACTGCCGCAGGCGGGGCTCGAGGATTTAG‑‑TCTGACTGTATCTCGGTCAGCCCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATACCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATg                                                                                                                            <  4:76728/145‑1 (MQ=255)
                            ctctTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGAAGTGATCCACATCATccgtac                                                                                                                                     <  1:301471/137‑5 (MQ=255)
                            ctctTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGAAGTGATCCACATCATccgtac                                                                                                                                     >  2:301471/1‑133 (MQ=255)
                             tctTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGt                                                                                                                        >  7:76008/1‑149 (MQ=255)
                               ttgcAGGCGCGGCTACTGCCGCAGGCGGGTCTCGAGGCTGTAG‑‑GCGTCTTGTTTGTCTGTCATCTCGATACGTTCGATGAGGGCCTGCTTTG‑GATATTGCTCGGTACGGATGATGTGGTTCACTTGATCCAGGTTGAGGCATGCAg                                                                                                                         <  8:316808/143‑1 (MQ=255)
                               ttCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCa                                                                                                                          <  1:2051/145‑1 (MQ=255)
                               ttCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGt                                                                                                                        <  8:76008/147‑1 (MQ=255)
                               ttCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGt                                                                                                                        >  3:170341/1‑147 (MQ=255)
                                tCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAg                                                                                                                     >  4:131382/1‑149 (MQ=255)
                                 ccAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGACAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAg                                                                                                                     <  2:75901/148‑1 (MQ=255)
                                      cgcgcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAAcc                                                                                                               <  3:385751/149‑1 (MQ=255)
                                       gcgcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATc                                                                                                                                          >  7:105471/1‑121 (MQ=255)
                                       gcgcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATc                                                                                                                                          <  8:105471/121‑1 (MQ=255)
                                       gcgcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCg                                                                                                              <  2:219230/149‑1 (MQ=255)
                                       gcgcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAAAGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATACAGGTTGAGGAATGCAGTGAGCAAAcct                                                                                                              >  8:354045/1‑148 (MQ=255)
                                         gcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCAc                                                                                                                                                >  7:53134/1‑113 (MQ=255)
                                         gcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCAc                                                                                                                                                <  8:53134/113‑1 (MQ=255)
                                         gcCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTc                                                                                                            <  1:33729/149‑1 (MQ=255)
                                              cTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACa                                                                                                       <  2:320616/149‑1 (MQ=255)
                                                  cGCAGGCGGGTCTCGAGGATGTTG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGTTGTGGATCACTTGATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTg                                                                                                   <  3:56420/149‑1 (MQ=255)
                                                  cGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTg                                                                                                   >  6:393806/1‑149 (MQ=255)
                                                   gCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGc                                                                                                  <  2:238117/149‑1 (MQ=255)
                                                                gAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCa                                                                                                                          <  6:173164/112‑1 (MQ=255)
                                                                 aGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGAc                                                                                    <  2:229388/149‑1 (MQ=255)
                                                                  ggATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGAcc                                                                                   <  4:170341/149‑1 (MQ=255)
                                                                     tGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTa                                                                                <  5:393806/149‑1 (MQ=255)
                                                                             cGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGa                                                                          <  3:131382/149‑1 (MQ=255)
                                                                              ggCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGa                                                                          <  2:193374/148‑1 (MQ=255)
                                                                                      tctcGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTc                                                                                                            >  4:24226/1‑106 (MQ=255)
                                                                                      tctcGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTc                                                                                                            <  3:24226/106‑1 (MQ=255)
                                                                                       ctcGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGAtgctg                                                                  >  4:383778/1‑147 (MQ=255)
                                                                                        tcGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCaaa                                                                                                                 >  5:349453/1‑99 (MQ=255)
                                                                                        tcGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGAGGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCaaa                                                                                                                 <  6:349453/99‑1 (MQ=255)
                                                                                          ggTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCa                                                                <  5:55903/146‑1 (MQ=255)
                                                                                           gTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCAAGACGATGCTGCAggcg                                                            >  3:49452/1‑149 (MQ=255)
                                                                                                 tggAAACGTTCGATCTGGGCCTGCTTGG‑GATATTCCTCGGTACGGATTATGTGGATCACTTCATCCAGGTTGAGGAATGCAGt                                                                                                                        <  2:367057/81‑1 (MQ=255)
                                                                                                   gAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACGGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGcgacg                                                      <  7:309963/147‑1 (MQ=255)
                                                                                                   gAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGaa                                                    >  4:51398/1‑149 (MQ=255)
                                                                                                       cGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGt                                                <  3:383778/149‑1 (MQ=255)
                                                                                                         ttCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGTATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTAt                                              <  5:132496/149‑1 (MQ=255)
                                                                                                            gATCAGGGCCTTCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATTCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTg                                            <  7:354045/148‑1 (MQ=255)
                                                                                                             aTCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGt                                          <  3:51398/149‑1 (MQ=255)
                                                                                                                 gggCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGAtgctg                                                                  <  8:291532/121‑1 (MQ=255)
                                                                                                                 gggCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGAtgctg                                                                  >  7:291532/1‑121 (MQ=255)
                                                                                                                 gggCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCg                                       >  1:57974/1‑148 (MQ=255)
                                                                                                                      tGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTc                                 >  4:188548/1‑149 (MQ=255)
                                                                                                                       gCTTTG‑GATATTGCTCTGTACTGATGATGTTTTTCACTTCATCCAGGTTGTGGAATGCATTTTGCAAACCGTCCAACTGGTGCAGACGCTTCTCGACCTTATCCAGACTATGCTGCAGGCGGCGACGAACGGTATTGTTGCGGAACt                                  <  2:368539/147‑1 (MQ=255)
                                                                                                                        cTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTc                                 <  4:387480/147‑1 (MQ=255)
                                                                                                                           gg‑gATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGTGACGAACGGTATTGGTGCGGAACTCCAGc                             >  2:52707/1‑148 (MQ=255)
                                                                                                                              gATATTCCTCGGTACGGATGATGTTGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACGGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACt                                  <  5:252326/141‑1 (MQ=255)
                                                                                                                              gATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGTGACGAACGGTATTGGTGCGGAACTCCAGCCAc                          <  1:52707/149‑1 (MQ=255)
                                                                                                                              gATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGAAGCTGAAGGCGGCGACGAACGGTATTCGTGCGGAACt                                  >  6:252326/1‑141 (MQ=255)
                                                                                                                                  ttCCTCGGTACGGATGAAGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTATCGCCCTTATCCAGACGATGCTGCAGGCGACGACCAAAGGTATTGCTGAGGAACTCCAGCCACTcca                      >  5:827/1‑149 (MQ=255)
                                                                                                                                   tCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAAc                                                                                                        <  4:2696/66‑1 (MQ=255)
                                                                                                                                   tCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAAc                                                                                                        >  3:2696/1‑66 (MQ=255)
                                                                                                                                          tACGTATTATGAGTATCACTTCATCCAGTTTGAGGAATGCAGTGAGCAAACCGTCGAAAAGTTGCAGACGCTTCTGGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTTGTGCGGAACTCCAGCCACTGCACCAATAgg              <  6:827/149‑1 (MQ=255)
                                                                                                                                           ctggATGAAGTGATCCACATCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCACTCCACCAATAggg             >  2:187681/3‑149 (MQ=255)
                                                                                                                                                    gtgGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCACTCCACCAATAGGGTTCGCAGGt    >  3:156357/1‑149 (MQ=255)
                                                                                                                                                      ggATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGAGCGGAACTCCAGCCACTCCACCAATAGGGTTCGCAGGttt  >  5:255008/1‑149 (MQ=255)
                                                                                                                                                      ggATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGAGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCACTCCACCAATAGGGTTCGCAGGttt  >  1:308911/1‑149 (MQ=255)
                                                                                                                                                       gATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCAc                          >  4:196941/1‑124 (MQ=255)
                                                                                                                                                       gATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCAc                          <  3:196941/124‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
ACTCATCCTGCTCACCGCGGATCTTCATCTCTTCCAGGCGCGCCAACTGCCGCAGGCGGGTCTCGAGGATGTAG‑‑TCGGCCTGTATCTCGGTCAGCTCGAAACGTTCGATCAGGGCCTGCTTGG‑GATATTCCTCGGTACGGATGATGTGGATCACTTCATCCAGGTTGAGGAATGCAGTGAGCAAACCGTCCAACAGGTGCAGACGCTTCTCGACCTTATCCAGACGATGCTGCAGGCGGCGACGAACGGTATTGGTGCGGAACTCCAGCCACTCCACCAATAGGGTTCGCAGGTTT  >  NC_002947/5584393‑5584688

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