Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,570,252 A→T 6.5% S80T (TCA→ACA)  gadC ← glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,570,2520AT 6.5% 188.6 / 2.2 61S80T (TCA→ACA) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base A (27/30);  new base T (2/2);  total (29/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.80e-01

ATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTC  >  NC_000913/1570107‑1570399
                                                                                                                                                 |                                                                                                                                                   
aTGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATAc                                                                                                                                                  >  2:338053/1‑149 (MQ=255)
   gggTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATAccc                                                                                                                                                <  1:338053/148‑1 (MQ=255)
      tCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGCGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGc                                                                                                                                            <  2:246670/149‑1 (MQ=255)
           attcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGa                                                                                                                                       <  8:326817/149‑1 (MQ=255)
            ttcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCTTCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGac                                                                                                                                      <  3:408142/149‑1 (MQ=255)
            ttcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGac                                                                                                                                      <  1:426958/149‑1 (MQ=255)
            ttcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGac                                                                                                                                      <  3:401787/149‑1 (MQ=255)
            ttcaGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGa                                                                                                                                       >  3:415364/1‑148 (MQ=255)
               aGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGacac                                                                                                                                    <  6:170004/148‑1 (MQ=255)
                  gcTGGCCATTTCAGGATGTAGGAGAGTGACCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTTCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACaccacc                                                                                                                                <  1:331918/149‑1 (MQ=255)
                         aTTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTcc                                                                                                                         >  5:273275/1‑149 (MQ=255)
                          tttCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGTTACCCAGGCGAAGACACCACCTTCTTcc                                                                                                                         >  2:269618/1‑148 (MQ=255)
                           ttCAGGATGTAGGATATTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTTCATTTAGCCAAATGATCTCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGa                                                                                                                                       <  5:82635/133‑1 (MQ=255)
                                     aGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACg                                                                                                             >  6:153702/1‑149 (MQ=255)
                                      ggAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACgg                                                                                                            >  2:398758/1‑149 (MQ=255)
                                        agagTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTg                                                                                                          <  6:372247/149‑1 (MQ=255)
                                            tGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCa                                                                                                      >  3:430886/1‑149 (MQ=255)
                                               cccTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCAttt                                                                                                   >  6:346026/1‑149 (MQ=255)
                                                 cTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGAt                                                                                                                                                    <  6:224106/98‑1 (MQ=255)
                                                 cTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTcc                                                                                                 <  6:273275/149‑1 (MQ=255)
                                                    acacGAAATAGAGCATCGGATTAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCg                                                                                                <  7:171280/147‑1 (MQ=255)
                                                           aTAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACaaa                                                                                       <  2:348385/149‑1 (MQ=255)
                                                              gagCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGTTACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGt                                                                                     <  1:269618/148‑1 (MQ=255)
                                                              gagCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTc                                                                                    >  4:222657/1‑149 (MQ=255)
                                                                          aaaaCCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATa                                                                         >  3:385819/1‑148 (MQ=255)
                                                                            aaCCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAAc                                                                      >  7:68339/1‑149 (MQ=255)
                                                                              ccAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCa                                                                    >  7:195670/1‑149 (MQ=255)
                                                                               cAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCATCCTTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCAt                                                                   <  3:426970/149‑1 (MQ=255)
                                                                               cAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCa                                                                    >  5:399931/1‑148 (MQ=255)
                                                                                    gCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATaaaa                                                               <  4:430886/148‑1 (MQ=255)
                                                                                       aTTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATccc                                                           <  5:153702/149‑1 (MQ=255)
                                                                                        tttGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCg                                                          <  3:305623/149‑1 (MQ=255)
                                                                                        tttGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTCATACCCAGGCGAAGACACCGCCTTCTTCCCAGCCGACGCCGGTCGGCTTTACCGCAGCACAATCTCCCCCGGGAAAAACCCAAAAAGTCACg                                                          >  7:431649/1‑149 (MQ=255)
                                                                                          tGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGc                                                         <  1:398758/148‑1 (MQ=255)
                                                                                            cAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCt                                                       >  6:343762/1‑148 (MQ=255)
                                                                                                  gCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAg                                                  >  2:290906/1‑147 (MQ=255)
                                                                                                  gCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCa                                                <  6:399931/149‑1 (MQ=255)
                                                                                                     aaaTGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACCCCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCCCAAAGTCCCACGGGAATAACCCATAAAATCCCGCCTAATAGCAgg                                              >  4:47446/1‑148 (MQ=255)
                                                                                                     aaaTGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAgg                                              >  3:432288/1‑148 (MQ=255)
                                                                                                          agaTCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGgaagaa                                         <  5:343762/148‑1 (MQ=255)
                                                                                                           gaTCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAagaaga                                       <  8:195670/149‑1 (MQ=255)
                                                                                                           gaTCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAagaaga                                       <  5:346026/149‑1 (MQ=255)
                                                                                                            aTCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGTTACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAagaaga                                       >  1:113080/1‑148 (MQ=255)
                                                                                                            aTCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACATTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCACCCGGGAATACACCATAAAATCCAGCCGAATAGCAGGAAGAacac                                      >  6:113832/1‑149 (MQ=255)
                                                                                                               gCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTaa                                   <  3:222657/149‑1 (MQ=255)
                                                                                                                cTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAAt                                  >  2:297051/1‑149 (MQ=255)
                                                                                                                   cAAATCGCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAAt                                  <  4:332092/146‑1 (MQ=255)
                                                                                                                   cAAATCCCCATCTCGGCCCCAGAGTATTTGTTACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATg                                 <  2:113080/147‑1 (MQ=255)
                                                                                                                    aaaTCACCATCTCGGCCACAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTACGCAGCAAAAAGTCCCACGGGAATAAAAAATAAAATCCCGACTAAAAGCAGGAAGAAGACTAATGaa                               >  8:79032/1‑148 (MQ=255)
                                                                                                                        ccccATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACCCCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATTCACCCTAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGcc                          >  2:61406/1‑149 (MQ=255)
                                                                                                                         cccATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGTAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGaaaa                             <  5:17075/145‑1 (MQ=255)
                                                                                                                         cccATCTCGGCCCCAGAGTATTTGATAACCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGCGAAAAAACCATAAAATCCCGCCAAATAGCAGGAAGAAGACTAATGaaaa                             >  6:17075/1‑145 (MQ=255)
                                                                                                                          ccATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGACGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGccc                         >  7:450786/1‑148 (MQ=255)
                                                                                                                                    cccAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATcc                                                            <  4:286527/103‑1 (MQ=255)
                                                                                                                                    cccAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATcc                                                            >  3:286527/1‑103 (MQ=255)
                                                                                                                                      cAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAgg            >  6:110830/1‑149 (MQ=255)
                                                                                                                                        gagTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTa          <  5:110830/149‑1 (MQ=255)
                                                                                                                                           tATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAggg       <  8:68339/149‑1 (MQ=255)
                                                                                                                                             tttGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTa     >  3:353471/1‑149 (MQ=255)
                                                                                                                                                gATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGaa              >  6:327025/1‑137 (MQ=255)
                                                                                                                                                gATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGaa              <  5:327025/137‑1 (MQ=255)
                                                                                                                                                 aTACCCAGGCGAAGACACCACCTTCTTCCCAGCCGGCGACGGTGGCCATTTCCGCAGCTCAAAGTCCCACTGGAATAAACCATAACATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTAtgc  <  1:61406/148‑3 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
ATGGGGTCTTCATTCAGCGCTGGCCATTTCAGGATGTAGGAGAGTGCCCCTAACACGAAATAGAGCATCGGAATAAAACCAATGGCGATTTGCAGATAGCCAAATGAGATCGCTGCAAATCCCCATCTCGGCCCCAGAGTATTTGATACCCAGGCGAAGACACCACCTTCTTCCCAGCCGTCGACGGTGGCCATTTCCGCAGCACAAAGTCCCACGGGAATAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTC  >  NC_000913/1570107‑1570399

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