Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,043,935 C→T 100% N335N (AAC→AAT)  PP_4448 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,043,9350CT94.2% 200.8 / 1.4 68N335N (AAC→AAT) PP_4448hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (34/30);  total (36/32)
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.93e-01
Rejected as polymorphism: E-value score below prediction cutoff.

TAACGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGTATAAC  >  NC_002947/5043794‑5044079
                                                                                                                                              |                                                                                                                                                
tAACGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATggagga                                                                                                                                            >  5:193993/1‑148 (MQ=255)
  aCGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGata                                                                                                                                          <  6:193993/148‑1 (MQ=255)
    gATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATgg                                                                                                                                                <  7:271547/140‑1 (MQ=255)
    gATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATgg                                                                                                                                                >  8:271547/1‑140 (MQ=255)
          aCATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATggag                                                                                                                                              >  3:298767/1‑136 (MQ=255)
          aCATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATggag                                                                                                                                              <  4:298767/136‑1 (MQ=255)
             tACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTc                                                                                                                              >  5:301425/1‑149 (MQ=255)
             tACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTc                                                                                                                              >  5:109517/1‑149 (MQ=255)
                gACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATggagga                                                                                                                                            >  6:82307/1‑132 (MQ=255)
                      ctaAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGAGCAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGAc                                                                                                                       >  1:78585/1‑148 (MQ=255)
                         aaggaaggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAAt                                                                                                                   >  4:311410/1‑148 (MQ=255)
                          aggaaggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTc                                                                                                                 <  8:294628/149‑1 (MQ=255)
                          aggaaggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTc                                                                                                                 <  3:311410/149‑1 (MQ=255)
                             aaggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTCGCTGGACAAATTCtat                                                                                                              <  6:109517/149‑1 (MQ=255)
                             aaggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTCGCTGGACAAATTCtat                                                                                                              <  6:301425/149‑1 (MQ=255)
                              aggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCtata                                                                                                             <  2:78585/149‑1 (MQ=255)
                                gTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACt                                                                                                           >  4:117475/1‑149 (MQ=255)
                                  aGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCg                                                                                                         >  4:52673/1‑149 (MQ=255)
                                    aaCAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCt                                                                                                       >  4:74361/1‑149 (MQ=255)
                                        gagCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGa                                                                                                   <  3:52673/149‑1 (MQ=255)
                                            cTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCta                                                                                                               >  6:363364/1‑133 (MQ=255)
                                             tACGACCGGTGTATCTCTAGACGA‑CAACGAACGGATATTTCGATACTATGAGAAAAAACTAGACTTCCCCCCAACCAAGCGAGAAAGCATATCAAATGCAGGATATGCAGATGGCTCGCTGGACAAATTCTAAAATAGCTTCTAGTCt                                                                                               >  4:48437/1‑148 (MQ=255)
                                             tACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGTATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCt                                                                                               >  3:257605/1‑148 (MQ=255)
                                               cGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTaaa                                                                                            <  3:117475/149‑1 (MQ=255)
                                                gACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAg                                                                                           >  1:290784/1‑149 (MQ=255)
                                                gACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAg                                                                                           >  3:293407/1‑149 (MQ=255)
                                                           ctctAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGaa                                                                                >  1:30528/1‑149 (MQ=255)
                                                           ctctAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGaa                                                                                >  7:18789/1‑149 (MQ=255)
                                                            tctAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAAc                                                                               >  2:19064/1‑149 (MQ=255)
                                                                  cga‑cAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTa                                                                          <  6:229248/148‑1 (MQ=255)
                                                                       aacgaacgGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGctt                                                                     >  4:157815/1‑149 (MQ=255)
                                                                       aacgaacgGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGct                                                                      >  4:99014/1‑148 (MQ=255)
                                                                          gaacgGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGcttctt                                                                  <  1:19064/149‑1 (MQ=255)
                                                                            acgGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGt                                                                >  5:321199/1‑149 (MQ=255)
                                                                                    ttCGATACTTTGAGAAAAAACTAGAATTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGATCGCTGGACAAATTCTATACTCGCTTCGAATCTATAGTTTCCCAGCAACGGTTAGGAATCTTCTTAGGagcg                                                        >  4:269309/1‑149 (MQ=255)
                                                                                     tCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCgg                                                       >  1:195326/1‑149 (MQ=255)
                                                                                      cGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCgg                                                       <  8:18789/148‑1 (MQ=255)
                                                                                       gATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGgc                                                     <  2:30528/149‑1 (MQ=255)
                                                                                            tttGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGt                                                                                          >  7:271568/1‑107 (MQ=255)
                                                                                            tttGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGt                                                                                          <  8:271568/107‑1 (MQ=255)
                                                                                             ttGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCa                                                 >  6:47598/1‑147 (MQ=255)
                                                                                              tGAGAAATAACTAGACTTCCACCGAAGCAAGCGAGAAAGCATATGAAATGGAGGATATGCTGATGGCTCGCTGGACAACTTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCaaa                                               <  3:48437/148‑1 (MQ=255)
                                                                                              tGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCaaa                                               <  2:195326/148‑1 (MQ=255)
                                                                                               gagaAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTa                                                                          >  1:258253/1‑120 (MQ=255)
                                                                                               gagaAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTa                                                                          <  2:258253/120‑1 (MQ=255)
                                                                                                 gaAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGaca                                           <  4:293407/149‑1 (MQ=255)
                                                                                                   aaaaaCTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGacac                                          <  3:99014/148‑1 (MQ=255)
                                                                                                    aaaaCTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGAcacca                                        <  5:47598/149‑1 (MQ=255)
                                                                                                        cTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGAT‑GCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACt                                    >  8:84194/1‑148 (MQ=255)
                                                                                                         tAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTc                                   >  3:347617/1‑149 (MQ=255)
                                                                                                         tAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATCGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTc                                   >  6:172561/1‑149 (MQ=255)
                                                                                                           gACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTc                                   <  2:290784/147‑1 (MQ=255)
                                                                                                              ttCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGacac                                          <  1:151722/137‑1 (MQ=255)
                                                                                                              ttCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGacac                                          >  2:151722/1‑137 (MQ=255)
                                                                                                              ttCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGaa                              >  5:49604/1‑149 (MQ=255)
                                                                                                              ttCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGaa                              >  6:362093/1‑149 (MQ=255)
                                                                                                              ttCCACCCAACCAAGCGAGAAAGCATATCAAATGGACGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAAc                                                                               >  7:116879/1‑100 (MQ=255)
                                                                                                              ttCCACCCAACCAAGAGAGAAAGCATATCAACTTGATGATATGCTTATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAATTTTCCCAGTAAc                                                                               <  8:116879/100‑1 (MQ=255)
                                                                                                                 cACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAacac                           <  4:347617/149‑1 (MQ=255)
                                                                                                                      aaccaaGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTcc                      <  3:74361/149‑1 (MQ=255)
                                                                                                                      aaccaaGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTcc                      <  3:157815/149‑1 (MQ=255)
                                                                                                                       accaaGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCg                     <  6:49604/149‑1 (MQ=255)
                                                                                                                               agaAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTc              >  4:336567/1‑148 (MQ=255)
                                                                                                                                  aaGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATa           <  3:336567/148‑1 (MQ=255)
                                                                                                                                   aGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAgg         <  5:172561/149‑1 (MQ=255)
                                                                                                                                     cATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGt       <  5:362093/149‑1 (MQ=255)
                                                                                                                                      atatCATATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGta      <  4:257605/149‑1 (MQ=255)
                                                                                                                                           aaaTGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGTATAAc  >  6:165978/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TAACGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGA‑CAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGTATAAC  >  NC_002947/5043794‑5044079

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