Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 465,367 G→A 5.9% L61L (CTG→TTG)  queC ← 7‑cyano‑7‑deazaguanine (preQ0) synthase; queuosine biosynthesis

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913465,3670GA 5.9% 222.2 / 2.6 68L61L (CTG→TTG) queC7‑cyano‑7‑deazaguanine (preQ0) synthase; queuosine biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (32/32);  new base A (2/2);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCCGCCAGCGTCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCT  >  NC_000913/465222‑465513
                                                                                                                                                 |                                                                                                                                                  
atTGCCGCCAGCGTCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCa                                                                                                                                                  >  8:326332/1‑148 (MQ=255)
atTGCCGCCAGCGTCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCa                                                                                                                                                  >  8:109653/1‑148 (MQ=255)
     cgccAGCGTCAGGAACAAAATATTACGCCCTGAGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTa                                                                                                                                             <  1:338124/148‑1 (MQ=255)
           cGTCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATgcgcg                                                                                                                                       <  7:79700/148‑1 (MQ=255)
           cGTCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATgcgcg                                                                                                                                       <  3:277600/148‑1 (MQ=255)
             tCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGCTACGATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCGACTGAACGCCAGCTCGTTGAGCAGGGAGACATCCAGCACCTTATGcgcgcg                                                                                                                                     >  1:282800/1‑148 (MQ=255)
                              cGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAgcg                                                                                                                    <  7:189293/148‑1 (MQ=255)
                                cccTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGcc                                                                                                                  <  1:365142/148‑1 (MQ=255)
                                   tGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGt                                                                                                               >  8:28843/1‑148 (MQ=255)
                                    gggACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCTAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGtt                                                                                                              >  8:315031/1‑148 (MQ=255)
                                      gACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCTAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTcg                                                                                                            <  7:315031/148‑1 (MQ=255)
                                            cGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGcc                                                                                                      <  6:167655/148‑1 (MQ=255)
                                            cGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGc                                                                                                       >  8:277706/1‑147 (MQ=255)
                                                ttCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGGGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACg                                                                                                   <  5:52439/147‑1 (MQ=255)
                                                        aCCATCGGCTTCAGGTTCGTAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCg                                                                                                >  1:383477/1‑142 (MQ=255)
                                                        aCCATCGGCTTCAGGTTCGTAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCg                                                                                                <  2:383477/142‑1 (MQ=255)
                                                        aCCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGtt                                                                                                                          >  5:92102/1‑116 (MQ=255)
                                                        aCCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGtt                                                                                                                          <  6:92102/116‑1 (MQ=255)
                                                           aTCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTgcg                                                                                       >  4:152918/1‑148 (MQ=255)
                                                           aTCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTgcg                                                                                       >  8:14107/1‑148 (MQ=255)
                                                            tCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTgcgc                                                                                      >  6:183945/1‑148 (MQ=255)
                                                             cGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCg                                                                                                >  5:182880/1‑137 (MQ=255)
                                                             cGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCg                                                                                                <  6:182880/137‑1 (MQ=255)
                                                              ggCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCCGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGa                                                                                    >  3:347926/1‑148 (MQ=255)
                                                                      gTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTg                                                                            <  7:170636/148‑1 (MQ=255)
                                                                           tAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATg                                                                                                                                           <  7:271280/80‑1 (MQ=255)
                                                                           tAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATg                                                                                                                                           >  8:271280/1‑80 (MQ=255)
                                                                             aTCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTaa                                                                     >  2:269637/1‑148 (MQ=255)
                                                                                 ggCACCGGAATGCTGTCACGCGTCAGGCTACTGCCCGGCAGCTCGTGGAGCAGGGTGACATCCAGCACCTTATGAGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCgc                                                                 <  2:282800/148‑2 (MQ=255)
                                                                                 ggCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGa                                                                 <  6:131986/148‑1 (MQ=255)
                                                                                 ggCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGa                                                                 <  7:326332/148‑1 (MQ=255)
                                                                                        gAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCAc                                                          >  1:107503/1‑148 (MQ=255)
                                                                                             cTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAAt                                                     >  1:288017/1‑148 (MQ=255)
                                                                                                 cACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGtt                                                                                                              >  3:197835/1‑87 (MQ=255)
                                                                                                 cACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGtt                                                                                                              <  4:197835/87‑1 (MQ=255)
                                                                                                 cACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGAc                                                 >  4:138760/1‑148 (MQ=255)
                                                                                                   cgcgTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACtt                                               >  5:33245/1‑148 (MQ=255)
                                                                                                      gTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTtcat                                            <  3:152918/148‑1 (MQ=255)
                                                                                                      gTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTtcat                                            <  7:14107/148‑1 (MQ=255)
                                                                                                         aGGCTACTGACCGCCAGCTCTTTGAGCAGTGTGAGATCCAGCACCTTATGCGCGCGTGCCCCAAGTTTCAGCGCCAGTTCGCGTGGCACGTCGATTTCTGCGCGATGCCTCTGACCGTAATCGAAAGTCACGCAATGGACTTcatcat                                         <  7:176023/148‑1 (MQ=255)
                                                                                                         aGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTtcatca                                          <  7:109653/147‑1 (MQ=255)
                                                                                                           gCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCatat                                       >  6:240991/1‑148 (MQ=255)
                                                                                                              aCTGACCGCCAGCTCGTTGAGCAGGGTGATATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATAttgt                                    >  5:73634/1‑148 (MQ=255)
                                                                                                              aCTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGGGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGAGGCCGCTGAACGTAATCGAACGTCAAGCCATGGACTTCCTCATAttgt                                    >  5:179946/1‑148 (MQ=255)
                                                                                                                  aCCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTCATCGAACGTCACGCAATGGACTTCATCATAttgttg                                  <  3:223648/146‑1 (MQ=255)
                                                                                                                  aCCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTa                                <  5:240991/148‑1 (MQ=255)
                                                                                                                  aCCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTa                                <  6:33245/148‑1 (MQ=255)
                                                                                                                  aCCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTa                                <  7:28843/148‑1 (MQ=255)
                                                                                                                  aCCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTa                                <  5:183945/148‑1 (MQ=255)
                                                                                                                   ccGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTaa                               >  1:116618/1‑148 (MQ=255)
                                                                                                                       cAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGcc                           >  6:297503/1‑148 (MQ=255)
                                                                                                                        aGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCt                          <  7:277706/148‑1 (MQ=255)
                                                                                                                        aGCTCGTTGAGCAGGGTGACATCCAACACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCt                          >  6:360373/1‑148 (MQ=255)
                                                                                                                        aGCTCGTTGAGCAGGGTGACATCCAACACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCt                          >  6:354763/1‑148 (MQ=255)
                                                                                                                             gatgcgCAGGGTGACATCCAGCCCCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGGTCGCGGGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCAc                      <  6:74668/142‑1 (MQ=255)
                                                                                                                                gAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAg                   <  2:107503/147‑1 (MQ=255)
                                                                                                                                       gTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAggtgg           >  6:56247/1‑148 (MQ=255)
                                                                                                                                        tGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCa                                                       <  6:6845/103‑1 (MQ=255)
                                                                                                                                        tGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCa                                                       >  5:6845/1‑103 (MQ=255)
                                                                                                                                        tGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAggtggt          <  4:143019/148‑1 (MQ=255)
                                                                                                                                             tCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGCCCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGAAATGCCTGCACCAGACAGGTGGTGGAAt     >  5:196921/1‑148 (MQ=255)
                                                                                                                                             tCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAAt     >  7:83670/1‑148 (MQ=255)
                                                                                                                                             tCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAAt     >  7:161353/1‑148 (MQ=255)
                                                                                                                                             tCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAAt     <  4:318774/148‑1 (MQ=255)
                                                                                                                                             tCCAACACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAAt     <  5:360373/148‑1 (MQ=255)
                                                                                                                                             tCCAACACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAAt     <  5:354763/148‑1 (MQ=255)
                                                                                                                                              caaGCACCTTATGCGCGCGTGCCCCCAGTTTCAGAGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGAAATGCCTGCACCAGACAGGTGGTGTAATc    >  5:136647/3‑148 (MQ=255)
                                                                                                                                                aGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGCCAGGTGGTGGAATCCt  >  5:4014/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATTGCCGCCAGCGTCAGGAACAAAATATTACGCCCTGGGACAAACGTATTCGGGATACCATCGGCTTCAGGTTCATAATCAGGCACCGGAATGCTGTCACGCGTCAGGCTACTGACCGCCAGCTCGTTGAGCAGGGTGACATCCAGCACCTTATGCGCGCGTGCCCCCAGTTTCAGCGCCAGTTCGCGTGCCACGTCGATTTCTGCGCGATGCCGCTGACCGTAATCGAACGTCACGCAATGGACTTCATCATATTGTTGTAATGCCTGCACCAGACAGGTGGTGGAATCCT  >  NC_000913/465222‑465513

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: ‑