Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,437,019 A→G 6.2% K89E (AAA→GAA)  cysQ → 3'(2'),5'‑bisphosphate nucleotidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,437,0190AG 6.2% 188.0 / 3.0 64K89E (AAA→GAA) cysQ3'(2'),5'‑bisphosphate nucleotidase
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base G (2/2);  total (29/35)
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.26e-01

GATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTG  >  NC_000913/4436875‑4437166
                                                                                                                                                |                                                                                                                                                   
gATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAaga                                                                                                                                                 <  1:211819/149‑1 (MQ=255)
       ccGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGTAAGTCCGTCAGCACTGGCAGCGTTACTTGCTGGTAGACCCGCTGGATGGTACTAAAGAGTttat                                                                                                                                           <  7:166301/148‑1 (MQ=255)
          cTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTaaa                                                                                                                                       >  8:12421/1‑149 (MQ=255)
          cTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTaaa                                                                                                                                       <  5:12119/149‑1 (MQ=255)
             acacCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAAc                                                                                                                                      <  3:148546/147‑1 (MQ=255)
               acCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTaa                                                                                                                                  <  7:239955/149‑1 (MQ=255)
                     aTCATGGACGGTTTAAGTACGCTGACACCGGATGTTGCGGTCCTTTATGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACAGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCg                                                                                                                             >  8:58210/1‑148 (MQ=255)
                      tCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGAACCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGaa                                                                                                                           <  1:127754/149‑1 (MQ=255)
                        aTGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAAtt                                                                                                                         <  2:39953/149‑1 (MQ=255)
                              ggTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGc                                                                                                                              <  3:172504/138‑1 (MQ=255)
                              ggTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCACCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGc                                                                                                                              >  4:172504/1‑138 (MQ=255)
                                     gtacgCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACAtt                                                                                                            <  7:37472/149‑1 (MQ=255)
                                     gtacgCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACAtt                                                                                                            >  3:92497/1‑149 (MQ=255)
                                           tGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGTCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCt                                                                                                       >  2:193658/1‑148 (MQ=255)
                                           tGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCt                                                                                                       >  3:102173/1‑148 (MQ=255)
                                             acacCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCACGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGCCCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATTGCGAATTCACCGTTAACATTGCGCTCAt                                                                                                    <  4:145960/149‑1 (MQ=255)
                                                 cGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGAc                                                                                                >  2:160064/1‑149 (MQ=255)
                                                 cGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGAc                                                                                                <  3:202043/149‑1 (MQ=255)
                                                       ttCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTa                                                                                                                                                     <  8:127410/90‑1 (MQ=255)
                                                       ttCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTa                                                                                                                                                     >  7:127410/1‑90 (MQ=255)
                                                          cGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCaaa                                                                                       >  3:72819/1‑149 (MQ=255)
                                                           ggTCCTTTCTGAAGAAGATCCTCCCGGGTGGGAAGTCCGGCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCGCATTGACCATGGCAAAc                                                                                      <  6:71813/149‑1 (MQ=255)
                                                                tttCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAtt                                                                                 >  7:33675/1‑149 (MQ=255)
                                                                   cTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGAtttt                                                                               <  1:160064/148‑1 (MQ=255)
                                                                     gaagaaGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTGAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAgg                                                                            >  3:87193/1‑149 (MQ=255)
                                                                     gaagaaGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTGAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAgg                                                                            >  3:83757/1‑149 (MQ=255)
                                                                     gaagaaGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAgg                                                                            >  2:193729/1‑149 (MQ=255)
                                                                        gaagaTCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgt                                                                         >  6:114985/1‑149 (MQ=255)
                                                                        gaagaTCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTGACTGGCTGGTAGACCCGGTGGATGGTACTAAAGAGTTTATTAAACGTATTGGCGAATTCACCGTTAACCTTGGGCTCATTGACCATGGCAAACCGAGTTTAGGCgt                                                                         >  6:64030/1‑149 (MQ=255)
                                                                         aagaTCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTGAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgtg                                                                        <  4:87193/149‑1 (MQ=255)
                                                                          agaTCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTGAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgtg                                                                        <  4:83757/148‑1 (MQ=255)
                                                                          agaTCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCgtgg                                                                       >  5:195336/1‑149 (MQ=255)
                                                                                tccCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCa                                                                                              >  3:179182/1‑120 (MQ=255)
                                                                                tccCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCa                                                                                              <  4:179182/120‑1 (MQ=255)
                                                                                         ggAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTaa                                                        >  4:34057/1‑149 (MQ=255)
                                                                                          gAAGTCCGTCAGCACTGGCAGCGGTACTGGCTGGTAGACCCGCTGGAGGGGACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCGCATTGACCATGGCAAACCGATTTTAGGCGTGGGGTATGCGCCGGTAAt                                                       <  3:201449/149‑1 (MQ=255)
                                                                                            aGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGa                                                     <  4:193803/149‑1 (MQ=255)
                                                                                                cGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGc                                                                           >  2:217465/1‑123 (MQ=255)
                                                                                                cGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGc                                                                           <  1:217465/123‑1 (MQ=255)
                                                                                                 gTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGt                                                 >  5:127570/1‑148 (MQ=255)
                                                                                                  tCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTaa                                               >  7:37069/1‑149 (MQ=255)
                                                                                                  tCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTaa                                               <  8:37069/149‑1 (MQ=255)
                                                                                                    aGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAAc                                                                                      >  1:38701/1‑108 (MQ=255)
                                                                                                    aGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAAc                                                                                      <  2:38701/108‑1 (MQ=255)
                                                                                                            cAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAgcgc                                     <  7:58210/149‑1 (MQ=255)
                                                                                                              gCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGATTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCgg                                   <  2:85084/149‑1 (MQ=255)
                                                                                                                gTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTgc                                                                                                          >  8:211123/1‑76 (MQ=255)
                                                                                                                gTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCgg                                   >  4:87604/1‑147 (MQ=255)
                                                                                                                     tggctggTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAgg                            >  3:38462/1‑149 (MQ=255)
                                                                                                                       gctggTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGc                           >  4:11516/1‑148 (MQ=255)
                                                                                                                       gctggTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGc                           <  5:114985/148‑1 (MQ=255)
                                                                                                                         tggTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCa                          >  5:130958/1‑147 (MQ=255)
                                                                                                                                 ccGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACa                                         >  1:183478/1‑124 (MQ=255)
                                                                                                                                 ccGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACa                                         <  2:183478/124‑1 (MQ=255)
                                                                                                                                   gCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaag              <  8:33675/149‑1 (MQ=255)
                                                                                                                                       gATGGTACTAAAGAGTTTATTCAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaagaa            <  6:195336/147‑1 (MQ=255)
                                                                                                                                       gATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaagaa            <  1:193729/147‑1 (MQ=255)
                                                                                                                                       gATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaagaa            <  3:87604/147‑1 (MQ=255)
                                                                                                                                       gATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAaagaa            <  4:92497/147‑1 (MQ=255)
                                                                                                                                       gATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAagaaga          <  1:193658/149‑1 (MQ=255)
                                                                                                                                          ggTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTg       >  1:75517/1‑149 (MQ=255)
                                                                                                                                            tACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCg     <  3:34057/149‑1 (MQ=255)
                                                                                                                                              cTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGgt   <  7:199211/149‑1 (MQ=255)
                                                                                                                                              cTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGgt   <  6:127570/149‑1 (MQ=255)
                                                                                                                                              cTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGgt   <  4:102173/149‑1 (MQ=255)
                                                                                                                                               tAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGgtg  >  8:168553/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
GATATTGCCGCTCACACCGTTATCATGGACGGTTTACGTACGCTGACACCGGATGTTCCGGTCCTTTCTGAAGAAGATCCTCCCGGTTGGGAAGTCCGTCAGCACTGGCAGCGTTACTGGCTGGTAGACCCGCTGGATGGTACTAAAGAGTTTATTAAACGTAATGGCGAATTCACCGTTAACATTGCGCTCATTGACCATGGCAAACCGATTTTAGGCGTGGTGTATGCGCCGGTAATGAACGTAATGTACAGCGCGGCAGAAGGCAAAGCGTGGAAAGAAGAGTGCGGTG  >  NC_000913/4436875‑4437166

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