Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,570,469 T→C T7T (ACA→ACG)  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,4690TC98.4% 234.2 / NA 64T7T (ACA→ACG) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (35/28);  total (36/28)
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.98e-01

TAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTT  >  NC_000913/1570327‑1570613
                                                                                                                                              |                                                                                                                                                
tAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTAc                                                                                                                                          >  1:435197/1‑151 (MQ=255)
 aaaCCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACt                                                                                                                                         >  1:930408/1‑151 (MQ=255)
 aaaCCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACt                                                                                                                                         <  1:384971/151‑1 (MQ=255)
      aTAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATg                                                                                                                                     >  1:373716/1‑150 (MQ=255)
          aaTCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGc                                                                                                                                 >  1:599595/1‑150 (MQ=255)
             cccGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCata                                                                                                                             <  1:698332/151‑1 (MQ=255)
              gggCCTAATAGCAGGACGAAGGCTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAACCATTAACATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCatat                                                                                                                            >  2:51138/3‑151 (MQ=255)
              ccGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGCTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCatat                                                                                                                            >  2:873745/1‑151 (MQ=255)
                     aTAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATcccc                                                                                                                     >  2:394945/1‑151 (MQ=255)
                       aGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCt                                                                                                                   >  1:389249/1‑151 (MQ=255)
                         cAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTa                                                                                                                  >  1:626538/1‑150 (MQ=255)
                           ggAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTaaaa                                                                                                               <  1:68819/151‑1 (MQ=255)
                             aagaagACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACg                                                                                                             >  2:991122/1‑151 (MQ=255)
                             aagaagACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACg                                                                                                             <  1:722533/151‑1 (MQ=255)
                                   aCTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTc                                                                                                       >  2:774641/1‑151 (MQ=255)
                                   aCTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTc                                                                                                       >  2:545011/1‑151 (MQ=255)
                                     tAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAACGGGAATTctt                                                                                                     >  2:983864/1‑149 (MQ=255)
                                      aaTGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCt                                                                                                     <  1:524377/150‑1 (MQ=255)
                                      aaTGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCt                                                                                                     <  1:132816/150‑1 (MQ=255)
                                           aaaGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGa                                                                                               <  1:638301/151‑1 (MQ=255)
                                            aaGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGa                                                                                               >  2:958977/1‑150 (MQ=255)
                                            aaGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGa                                                                                               >  1:924560/1‑150 (MQ=255)
                                               cccGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAAcc                                                                                            <  2:610208/150‑1 (MQ=255)
                                               cccGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCg                                                                                           >  1:101129/1‑151 (MQ=255)
                                                  gATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACGCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCAc                                                                                        >  2:694478/1‑151 (MQ=255)
                                                  gATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCAc                                                                                        >  2:590707/1‑151 (MQ=255)
                                                   aTGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCAcc                                                                                       >  2:328/1‑151 (MQ=255)
                                                        gCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGttt                                                                                  <  1:285873/151‑1 (MQ=255)
                                                           aaGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGtttt                                                                                 <  2:648867/149‑1 (MQ=255)
                                                            aGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGtt                                                                              >  2:144598/1‑151 (MQ=255)
                                                               tAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTg                                                                           >  1:744748/1‑151 (MQ=255)
                                                                 gggTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGa                                                                          >  1:731068/1‑150 (MQ=255)
                                                                  ggTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAg                                                                        <  2:475606/151‑1 (MQ=255)
                                                                       tCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTTGGAACCGCACCTGTTTTGTTCTGAAGCGTAt                                                                   >  2:971985/1‑151 (MQ=255)
                                                                       tCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTAt                                                                   >  2:292357/1‑151 (MQ=255)
                                                                         aTAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCCCCTGTTTTGTTCTGAAGCGTAt                                                                   <  1:106734/149‑1 (MQ=255)
                                                                                 ccATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAat                                                         <  1:888012/151‑1 (MQ=255)
                                                                                 ccATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAat                                                         >  1:422635/1‑151 (MQ=255)
                                                                                  cATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGGCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAat                                                         >  2:103840/1‑150 (MQ=255)
                                                                                     tACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATAtt                                                      <  1:426229/150‑1 (MQ=255)
                                                                                      aCCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATAtttt                                                    <  1:633906/151‑1 (MQ=255)
                                                                                            gATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTg                                              <  2:766431/151‑1 (MQ=255)
                                                                                            gATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTg                                              >  1:537389/1‑151 (MQ=255)
                                                                                                     aTGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGcc                                      >  1:33527/1‑150 (MQ=255)
                                                                                                       ggCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGt                                   <  2:639232/151‑1 (MQ=255)
                                                                                                        gCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTg                                  <  2:930323/151‑1 (MQ=255)
                                                                                                           aaGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAAc                               <  1:850748/151‑1 (MQ=255)
                                                                                                           aaGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAAc                               <  2:64615/151‑1 (MQ=255)
                                                                                                              aaTCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGa                            <  1:651077/151‑1 (MQ=255)
                                                                                                               aTCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGAc                           >  2:475168/1‑151 (MQ=255)
                                                                                                               aTCCAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGAc                           <  1:677904/151‑1 (MQ=255)
                                                                                                               aTCCAAGGAATGGGAGCTGCCTAGCTTTTACCGTCTGTTCTGAAGTAACCATTTTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGAc                           <  1:983949/151‑1 (MQ=255)
                                                                                                                 ccAAGTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACt                          <  1:683207/150‑1 (MQ=255)
                                                                                                                     gTAATGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTgg                      >  1:710215/1‑150 (MQ=255)
                                                                                                                       aaTGTGAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGaa                    <  2:756363/150‑1 (MQ=255)
                                                                                                                            gAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTc              <  1:955205/151‑1 (MQ=255)
                                                                                                                            gAGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTc              <  1:185675/151‑1 (MQ=255)
                                                                                                                             aGCTGCTTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTcc             >  1:772737/1‑151 (MQ=255)
                                                                                                                               ctgctTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTc              >  1:198082/1‑148 (MQ=255)
                                                                                                                                tgctTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCg           >  1:6657/1‑150 (MQ=255)
                                                                                                                                 gctTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGa          <  2:626473/150‑1 (MQ=255)
                                                                                                                                  ctTAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGaa         >  1:479023/1‑150 (MQ=255)
                                                                                                                                    tAGCTTTACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACg      >  1:885939/1‑151 (MQ=255)
                                                                                                                                        tttACCCGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGtt  >  2:44610/1‑151 (MQ=255)
                                                                                                                                              |                                                                                                                                                
TAAACCATAAAATCCCGCCTAATAGCAGGAAGAAGACTAATGAAAAGCCCGATGTTGCGAAGGTAGGGTATTCATAAACAGCCATTACCATCGATGCCGTTATGGCAAAGAATCCAAGTAATGTGAGCTGCTTAGCTTTACCTGTCTGTACTGATGTAGCCATATTATCCCCCTAAAACGGTATTCCTGTCGGAACCGCACCTGTTTTGTTCTGAAGCGTATTCAGAACAATATTTTCCGTTGCTAATGCCAGTGAACAGACTTTGGAAATTGTCCCGAAACGGGTT  >  NC_000913/1570327‑1570613

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑