Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 568,152 G→A 47.4% pseudogene (118/213 nt) renD → pseudogene, DLP12 prophage;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913568,1520GA47.4% 5.4 / 79.0 76pseudogene (118/213 nt)renDpseudogene, DLP12 prophage;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (20/20);  new base A (18/18);  total (38/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.41e-01

GGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACA  >  NC_000913/568004‑568300
                                                                                                                                                    |                                                                                                                                                    
ggCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAgg                                                                                                                                                      <  7:123555/149‑1 (MQ=255)
 gCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGGCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGt                                                                                                                                                     <  7:60394/149‑1 (MQ=255)
         cgcgctaacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCGCGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                                                                                                                              >  4:329098‑M2/24‑148 (MQ=255)
         ccATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACg                                                                                                                                              >  4:90538/1‑148 (MQ=255)
           atatTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAAt                                                                                                                                           >  1:106493/1‑149 (MQ=255)
              taacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                                <  4:235245‑M2/123‑1 (MQ=255)
              taacaggcgcaacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                                >  3:235245‑M2/19‑141 (MQ=255)
                      gTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTaa                                                                                                                                <  7:330550/149‑1 (MQ=255)
                        aacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                                       >  7:121247/7‑140 (MQ=255)
                        aacagtaaCCAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                                       <  8:121247/134‑1 (MQ=255)
                           aTCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgt                                                                                                                           >  7:391951/1‑149 (MQ=255)
                           aTCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgt                                                                                                                           >  7:386640/1‑149 (MQ=255)
                            tCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAgtg                                                                                                                          <  3:90538/149‑1 (MQ=255)
                                ccAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                                  >  6:82422/1‑133 (MQ=255)
                                ccAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                                  <  5:82422/137‑5 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                                   >  4:236636/1‑119 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                                   <  3:236636/119‑1 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaac                                                                                                                                           <  4:8965/127‑2 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                                      <  6:134936/132‑1 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                                      >  5:134936/1‑132 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGCTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaac                                                                                                                                           >  3:211570/1‑126 (MQ=255)
                                 cAGCATAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGCTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaac                                                                                                                                           <  4:211570/127‑2 (MQ=255)
                                     aTAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCagag                                                                                                                 >  6:21128/1‑149 (MQ=255)
                                     aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtga                                                                                                                                    >  5:341955/1‑128 (MQ=255)
                                     aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtga                                                                                                                                    <  6:341955/130‑3 (MQ=255)
                                     aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                                      <  5:321363/128‑1 (MQ=255)
                                     aTAAATCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                                      >  6:321363/1‑128 (MQ=255)
                                         aTCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCAc                                                                                                                                               >  7:24685/1‑115 (MQ=255)
                                         aTCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCAc                                                                                                                                               <  8:24685/115‑1 (MQ=255)
                                         aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                                <  7:315501/114‑1 (MQ=255)
                                         aTCAGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                                >  8:315501/1‑114 (MQ=255)
                                          tCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCcg                                                                                                            <  1:324612/149‑1 (MQ=255)
                                          tCAGGCCGCGGCTAAAATGACACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                                >  6:54778/1‑113 (MQ=255)
                                          tCAGGCCGCGGCTAAAATGACACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCa                                                                                                                                                <  5:54778/113‑1 (MQ=255)
                                            aGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                                   >  7:278153/1‑108 (MQ=255)
                                            aGGCCGCGGCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGa                                                                                                                                                   <  8:278153/108‑1 (MQ=255)
                                                   ggCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtg                                                                                                                                     >  6:230930/1‑114 (MQ=255)
                                                   ggCTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtg                                                                                                                                     <  5:230930/115‑2 (MQ=255)
                                                    gCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGTGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACg                                                                                                  <  3:36282/149‑1 (MQ=255)
                                                     cTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACgg                                                                                                 <  8:203872/149‑1 (MQ=255)
                                                     cTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                                      >  1:350024/1‑112 (MQ=255)
                                                     cTAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAt                                                                                                                                      <  2:350024/112‑1 (MQ=255)
                                                      tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGACCGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                                       <  5:2046/110‑1 (MQ=255)
                                                      tAAAATGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGa                                                                                                                                       >  6:2046/1‑110 (MQ=255)
                                                        aaaTGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                                  <  8:57112/113‑5 (MQ=255)
                                                        aaaTGGCACGGGCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGAACCTGAtgaac                                                                                                                                  >  7:57112/1‑109 (MQ=255)
                                                                cgGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGt                                                                                      <  4:205542/149‑1 (MQ=255)
                                                                 gggCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                                                                                                                              >  8:240360/1‑92 (MQ=255)
                                                                 gggCAGGTCTTCTGGTTATCGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACg                                                                                                                                              <  7:240360/92‑1 (MQ=255)
                                                                    cAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGtatat                                                                                  <  5:21128/149‑1 (MQ=255)
                                                                        aaTCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATaaaa                                                                              <  6:34575/149‑1 (MQ=255)
                                                                              ggTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACAt                                                                        >  4:204262/1‑149 (MQ=255)
                                                                               gTCGTTGAAGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTTATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATc                                                                       <  3:317294/149‑1 (MQ=255)
                                                                                  gTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAAt                                                                    <  2:127901/149‑1 (MQ=255)
                                                                                      aTGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATggg                                                                 <  2:106493/148‑1 (MQ=255)
                                                                                                   gCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaacc                                                                                                                                          >  8:393197/1‑60 (MQ=255)
                                                                                                   gCGAACGGTGTATTACCGGTTTGCTACCAGGGAAGAACGGGAAGGAAAGATGAGCACGaacc                                                                                                                                          <  7:393197/62‑3 (MQ=255)
                                                                                                   gCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAatgatg                                                     <  3:204262/147‑1 (MQ=255)
                                                                                                         ggTGTATTACCTGTTCGCTACCAGAGAAGAATGGGAAGGAACTGTGAGCACGACTCTGATTTTTAAGGAGTGTCTCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAATAACATCAATGGGAACACAATGATGAAACAgtt                                             <  4:287114/149‑3 (MQ=255)
                                                                                                                     gTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTc                                  >  4:50956/1‑148 (MQ=255)
                                                                                                                           tACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGa                                                                                                                             <  5:29758/51‑1 (MQ=255)
                                                                                                                           tACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGa                                                                                                                             >  6:29758/1‑51 (MQ=255)
                                                                                                                             ccAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag                         >  6:234568/1‑149 (MQ=255)
                                                                                                                             ccAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag                         >  1:31538/1‑149 (MQ=255)
                                                                                                                              cAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCaaa                               >  8:381583/1‑142 (MQ=255)
                                                                                                                              cAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCaaa                               <  7:381583/142‑1 (MQ=255)
                                                                                                                                gagaAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAACCGGGTAGTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGCGTTCAAACTgtagta                       >  1:324151/1‑148 (MQ=255)
                                                                                                                                  gaagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag                         >  1:308802/1‑144 (MQ=255)
                                                                                                                                  gaagaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTgtag                         <  2:308802/144‑1 (MQ=255)
                                                                                                                                     gaaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtc                 >  8:158995/1‑149 (MQ=255)
                                                                                                                                      aaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGtgagttga                                      >  4:239225/1‑127 (MQ=255)
                                                                                                                                      aaTGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGtgagttga                                      <  3:239225/127‑1 (MQ=255)
                                                                                                                                          ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATtctctc             <  3:50956/148‑1 (MQ=255)
                                                                                                                                          ggaaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACAACAATGTGAACCCAATGATGAAAAAGGTGAGTAGAGTTCAAACTGTAGTACAAGtctctc             >  6:94461/1‑148 (MQ=255)
                                                                                                                                            aaggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAg          >  8:61573/1‑149 (MQ=255)
                                                                                                                                              ggaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATTAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGtt        >  7:304222/1‑149 (MQ=255)
                                                                                                                                               gaaAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGttt       >  2:72668/1‑149 (MQ=255)
                                                                                                                                                   ggTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAAc   <  5:234568/149‑1 (MQ=255)
                                                                                                                                                    gTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACa  <  1:310303/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCCTGTGTCCATATTACGTGGGTAGGATCAACCAGCATAAATCAGGCTGCGGCTAAAATGGCGCGGGCAGGAATCCTGGTCGTTGATGGTAAGGTCTGGCGAACGGTGTATTACCGGTTCGCTACCAGAGAAGAATGGGAAGGAAAGGTGAGCACGAATCTGATTTTTAAGGAGTGTCGCCAGAGTGCCGCGATGAAACGGGTATTGAGGGTATATAAAAGAACATCAATGGGAACACAATGATGAAACAGGTGAGTTGAGTTCAAACTGTAGTACAATTCTCTCCAGTTTGAACA  >  NC_000913/568004‑568300

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: