Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,633,920 T→C T51A (ACC→GCC)  guaB ← inosine 5'‑monophosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,633,9200TC98.5% 221.3 / ‑4.4 65T51A (ACC→GCC) guaBinosine 5'‑monophosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (34/30);  total (35/30)
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.22e-01

ACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGTATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTT  >  NC_000913/2633793‑2634054
                                                                                                                               |                                                                                                                                      
acaccaGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATtgctgc                                                                                                                            >  2:611425/1‑140 (MQ=255)
acaccaGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATtgctgc                                                                                                                            >  2:367205/1‑140 (MQ=255)
   ccaGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGa                                                                                                                         <  2:698732/140‑1 (MQ=255)
   ccaGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGa                                                                                                                         <  1:630699/140‑1 (MQ=255)
    caGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGaa                                                                                                                        >  2:217301/1‑140 (MQ=255)
      gATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAg                                                                                                                      >  2:140740/1‑140 (MQ=255)
         tCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGa                                                                                                                         <  2:43087/134‑1 (MQ=255)
         tCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGa                                                                                                                         >  1:43087/1‑134 (MQ=255)
            tgtTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAgg                                                                                                                >  1:393281/1‑140 (MQ=255)
            tgtTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAgg                                                                                                                >  1:544460/1‑140 (MQ=255)
                  ttCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATAtt                                                                                                          <  2:532041/140‑1 (MQ=255)
                     acacGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAg                                                                                                       <  2:757762/140‑1 (MQ=255)
                         ggCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGa                                                                                                   >  2:916825/1‑140 (MQ=255)
                         ggCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGa                                                                                                   >  1:534583/1‑140 (MQ=255)
                           cGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGaa                                                                                                  >  1:158392/1‑139 (MQ=255)
                           cGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGaa                                                                                                  <  2:158392/139‑1 (MQ=255)
                                     ctGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTc                                                                                       >  2:763054/1‑140 (MQ=255)
                                        ccTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGc                                                                                    >  2:454458/1‑140 (MQ=255)
                                         cTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCt                                                                                   >  2:919107/1‑140 (MQ=255)
                                            gCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTggg                                                                                >  1:327777/1‑140 (MQ=255)
                                            gCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTggg                                                                                >  2:497936/1‑140 (MQ=255)
                                                tCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCt                                                                            <  1:611425/140‑1 (MQ=255)
                                                tCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCt                                                                            <  1:497936/140‑1 (MQ=255)
                                                tCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCt                                                                            <  2:534583/140‑1 (MQ=255)
                                                  aaTGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGa                                                                          >  2:685002/1‑140 (MQ=255)
                                                    tGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAgg                                                                        >  1:906568/1‑140 (MQ=255)
                                                     ggACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGt                                                                       <  1:672959/140‑1 (MQ=255)
                                                        cATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTcag                                                                    <  1:339837/140‑1 (MQ=255)
                                                          tGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTcagca                                                                  <  2:76652/140‑1 (MQ=255)
                                                              tttGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTAt                                                              >  1:857838/1‑140 (MQ=255)
                                                                   ggATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGc                                                         >  2:516707/1‑140 (MQ=255)
                                                                    gATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCa                                                        <  1:516707/140‑1 (MQ=255)
                                                                      tAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGa                                                      >  2:147536/1‑140 (MQ=255)
                                                                       aaaGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGaa                                                     <  1:916825/140‑1 (MQ=255)
                                                                         aGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACg                                                   >  2:881723/1‑140 (MQ=255)
                                                                         aGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACg                                                   <  1:454458/140‑1 (MQ=255)
                                                                         aGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACg                                                   <  1:919107/140‑1 (MQ=255)
                                                                            cGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTa                                                <  2:593753/140‑1 (MQ=255)
                                                                              aTACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTaga                                              >  2:210529/1‑140 (MQ=255)
                                                                              aTACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTaga                                              <  1:192246/140‑1 (MQ=255)
                                                                               tACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGaa                                                     >  2:704389/1‑132 (MQ=255)
                                                                               tACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGaa                                                     <  1:704389/132‑1 (MQ=255)
                                                                                     cTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCa                                       <  1:685002/140‑1 (MQ=255)
                                                                                       tCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAgg                                     <  1:367205/140‑1 (MQ=255)
                                                                                         cTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGaa                                   <  1:151471/140‑1 (MQ=255)
                                                                                                 gagCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGgaga                           >  2:507042/1‑140 (MQ=255)
                                                                                                 gagCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGgaga                           >  2:482282/1‑140 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGTATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                          >  1:189119/1‑140 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGa                                    >  2:663188/1‑130 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGa                                    <  1:663188/130‑1 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                          >  1:785012/1‑140 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                          >  1:496459/1‑140 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                          >  1:403172/1‑140 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                          >  2:353714/1‑140 (MQ=255)
                                                                                                  agCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGaa                          >  2:185805/1‑140 (MQ=255)
                                                                                                   gCAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAc                         <  1:881723/140‑1 (MQ=255)
                                                                                                    cAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcg                        <  1:507042/140‑1 (MQ=255)
                                                                                                    cAATAGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcg                        <  1:482282/140‑1 (MQ=255)
                                                                                                        aGCCAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAAcgtcgt                    >  1:939826/1‑140 (MQ=255)
                                                                                                           cAGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaa                 <  2:544460/140‑1 (MQ=255)
                                                                                                            aGGCGCGCTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCaaa                >  2:830008/1‑140 (MQ=255)
                                                                                                                 cgcTTCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCa           >  2:771259/1‑140 (MQ=255)
                                                                                                                    ttCCGTTACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCagag        <  1:771259/140‑1 (MQ=255)
                                                                                                                          tACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt  <  2:857838/140‑1 (MQ=255)
                                                                                                                          tACGGCATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGcttctt  <  2:939826/140‑1 (MQ=255)
                                                                                                                               |                                                                                                                                      
ACACCAGATTCGTGTTTTTTCACACGGCGAACTTCTTCTGCCTGGCGTTCAATGGACATGTTTTTGTGGATAAAGCCGATACCGCCTTCCTGAGCCAGAGCAATAGCCAGGCGCGCTTCCGTTACGGTATCCATTGCTGCGGAAAGCATAGGGATATTCAGACGAATAGTTTTCGTCAGCTGGGTGCTGAGGTCAGCAGTATTCGGCAGAACGGTAGAGTGAGCAGGAACGAGGAGAACGTCGTCAAACGTCAGAGCTTCTT  >  NC_000913/2633793‑2634054

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