Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 459,834 C→A A316E (GCG→GAG)  lon → Lon protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,8340CA98.5% 185.0 / NA 70A316E (GCG→GAG) lonLon protease
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/34);  total (36/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

GAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTA  >  NC_000913/459701‑459971
                                                                                                                                      |                                                                                                                                         
gAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAAt                                                                                                                                  >  1:890125/1‑143 (MQ=255)
    cAGCGTTGCAGAAGCTGAAACTGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCtt                                                                                                                              <  1:121934/143‑1 (MQ=255)
     agagTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTg                                                                                                                             <  2:55971/143‑1 (MQ=255)
      gagTTTCAGAAGCTGAAAAGGATGTCTCCGATGTCGGAAGAAGCGACCGTAGGGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGa                                                                                                                            <  1:1284684/143‑1 (MQ=255)
       agTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTGGTGCGTGGTTATATCGACTGGATGGGACAGGTGCCGTGGAATGCGCGTAGCAAGGGCAAAAAAGACCTGCGTCAGGAGCCGGCAAGCCTTGAt                                                                                                                           >  1:96019/1‑143 (MQ=255)
         ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATAc                                                                                                                         >  2:545462/1‑143 (MQ=255)
            cAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGa                                                                                                                      <  1:631935/143‑1 (MQ=255)
             aGAAGCCGAAAATGATGTCCCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGAc                                                                                                                     <  1:556922/143‑1 (MQ=255)
                   tGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTAt                                                                                                               >  1:652139/1‑143 (MQ=255)
                    gAAGATGCTGTCTGCGATGTCGGCTGAAGGGACGGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATg                                                                                                              <  2:821165/143‑1 (MQ=255)
                    gAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTATAGGTGAAGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGAACATCCTTGATACCGACCATTATg                                                                                                              >  2:174678/1‑143 (MQ=255)
                            tgTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGGCTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGAGCATTATGGTCTGGAg                                                                                                      >  1:1029908/1‑143 (MQ=255)
                             gTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAACGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATCGTCTGGAgc                                                                                                     >  2:246416/1‑143 (MQ=255)
                              tctcCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGGGGAATGCGCGTAGCAAGGTCGAAAACGACCTGCGCCAGGAGCAGGAAATCCTTGATACCGGCCATTAGGGTCggggggg                                                                                                    >  1:880016/1‑139 (MQ=255)
                              tctcCGACGTCCGCAGAAGCGACCGTAGTGCGCGGTTATATCGACTTGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAgcg                                                                                                    <  1:1468444/143‑1 (MQ=255)
                                tctgATGTCGGCAGAAGCGACCGGAGTGCGTGGTTGTATCGACTGGATGGTACAGGTGCCGGGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGcgcg                                                                                                  <  2:1253854/140‑1 (MQ=255)
                                tcCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGGGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGGGCAGGAAATCCTTGATACCGCCCACTATGGTCCGGAGcgcg                                                                                                  >  2:716294/1‑143 (MQ=255)
                                 ccGGTGTCGGCAGAAGCGACCGTAGTGCGTGGTGAGAGCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGt                                                                                                 <  2:343933/143‑1 (MQ=255)
                                  cGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCTGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTg                                                                                                >  2:776129/1‑143 (MQ=255)
                                    aTGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGaa                                                                                              >  1:834161/1‑143 (MQ=255)
                                     tGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGaaa                                                                                             <  2:216751/143‑1 (MQ=255)
                                       tCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATACTTGATACCGACCATCATGGTCTGGAGCGCGTGAAAGa                                                                                           >  2:853810/1‑143 (MQ=255)
                                        cGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGAt                                                                                          <  2:890125/143‑1 (MQ=255)
                                         tgaaGACGCGACTGTAATGAGTGGTTATATTGACTGGATGCTACAGGTGCCGTGGAATGCGCGTAGCAAGAGCAAAGAATACTTGACTTGGCAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                         <  2:518552/140‑1 (MQ=255)
                                         tccAGAAGCGCCCGTAGTGCGTGGTTATATCGACTGGCTGGCGCAGGTGCCGTGGAATTCCTTCATCTTAGTCAACAACGACCCGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                         <  1:1458548/141‑1 (MQ=255)
                                         ggaggAAGCGATCGTGGTGCGTGGTTATATCGAGAGAATGGTAGAGGTGGCGTGGAATGCGCGTAGCAAGGTCAAAAAAGATCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                         <  2:884242/139‑1 (MQ=255)
                                         ggCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTGCAGGGCGTAAAATCCTTGATACCGACCGTTATGGTCTGGCGCGCGTGAAAGATc                                                                                         >  2:1039022/1‑143 (MQ=255)
                                         ggCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                         <  2:457670/143‑1 (MQ=255)
                                          gCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCg                                                                                        >  2:1200143/1‑143 (MQ=255)
                                                gCGACCGTAGTGCGTGGGTATATCGACGGGATGGTACAGGTGCCGTGGAATGCGCGGGGCAAGGGCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGCGAAAGATCGAATCCt                                                                                  >  1:1197962/1‑143 (MQ=255)
                                                   tgcggaGTGCGTGGTTGTATCGACTGGGTGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGGCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGGGTGGAGCGCGTGAAAGATCGAATCCTTGa                                                                               <  2:681652/138‑1 (MQ=255)
                                                   gaggtggTGCGTGGTTATATCGACTGGAGGGTACAGGTGCCTGGGAATGCGCGTAGCAAGGTCATAAAAGCCCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGa                                                                               <  2:173630/137‑1 (MQ=255)
                                                   aCCGTAGTGGGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGa                                                                               <  2:32297/143‑1 (MQ=255)
                                                       tAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCGTGCGTAt                                                                           >  1:197245/1‑143 (MQ=255)
                                                          tGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGGCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGGTCGAACGCTCGAGTATCtt                                                                        >  2:198061/1‑143 (MQ=255)
                                                           gggtagTTATCTCGACTGAATAGTACAGGTGCCATGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTg                                                                       <  2:660243/138‑1 (MQ=255)
                                                            cGTGGTTATATCGACTGGAGGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCGGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGc                                                                      <  2:608045/143‑1 (MQ=255)
                                                               ggTCATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGt                                                                   >  2:815509/1‑143 (MQ=255)
                                                                 tgatgtCGACTTGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTc                                                                 <  1:1640/138‑1 (MQ=255)
                                                                  tataTCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAGGATCGAATCCTTGAGTATCTTGCGGGTCa                                                                >  2:403857/1‑143 (MQ=255)
                                                                    tatCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAATGACCTGCGTCAGGAGCAGGAATTTCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGACTCCTTGAGTATCTTGCGGTTCaaa                                                              <  1:362217/143‑1 (MQ=255)
                                                                     atCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGGCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTGTGGGCTGGAGCGCGTGATAGATCGGAGCCTTGGGTATTTTGCGGTTCGAAg                                                             >  1:496363/1‑143 (MQ=255)
                                                                      tCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCCAGGTCAAAAAAGACCTGCGTCAGGCGCAGGCAATCCTTGATACCGACCATTATGGACTGGAGCGCGTGAAAGATCAAATCCTTGAGTCACTTGCGGGTCAAAGc                                                            >  2:1152721/1‑143 (MQ=255)
                                                                      tCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGc                                                            >  1:1469902/1‑143 (MQ=255)
                                                                         aCTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCgt                                                         >  1:235731/1‑143 (MQ=255)
                                                                           gggAGGGTACAGGTGCCGTGGAATGCGCGTAGCGAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCgtgt                                                       <  2:516090/142‑1 (MQ=255)
                                                                              aTGGTACAGGTGCGGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTcaa                                                    <  2:1443204/143‑1 (MQ=255)
                                                                              aTGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACTTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTcaa                                                    >  2:102015/1‑143 (MQ=255)
                                                                              aTGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGACTCCTTGAGTATCTTGCGGTTCAAAGCCGTGTcaa                                                    >  2:783636/1‑143 (MQ=255)
                                                                              aTGGTAACGGTGGCGTGGAATGCGAGCAGCAAGGTCGAAAAAGCCCTGCGTCAGGACGAGGAAATCCTAGAGACCGACGATGACGGTCTGGGGCGCGTGAAAGCTCGAATCCGTGAGTATCTTGCGGTTCAAAGCCGTGTcaa                                                    <  2:511419/143‑1 (MQ=255)
                                                                                   aCAGGTGCCGTGGAATGCGCGTAGCAAGGACAAAACAGGCCTGCGTGAGGAGCAGGAAATCCTTGCTATGAGGTATTATGGTCTGGAGCGCGGGAAGGATCGAGTGCTTGAGTATCTTGCGGTTCAAGGCGGTGTCAACaaaa                                               >  2:1042569/1‑143 (MQ=255)
                                                                                    cAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAAt                                              <  1:1125960/143‑1 (MQ=255)
                                                                                     aGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGGTCAAAGCCGTGTCAACAAAATc                                             >  1:1026331/1‑143 (MQ=255)
                                                                                       gTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTCTGTTGCGGTTCAAAGCCGTGTCAACAAAATcga                                           >  1:392660/1‑141 (MQ=255)
                                                                                       gTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCaa                                           <  1:1040932/143‑1 (MQ=255)
                                                                                       gTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGACCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCaa                                           <  1:9800/143‑1 (MQ=255)
                                                                                        tGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGCCAACAAAATCaac                                          >  2:327853/1‑142 (MQ=255)
                                                                                            gTGGAATGCGCGTAGCAAGGTCAAAAAGGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTc                                                                 <  2:1130185/116‑1 (MQ=255)
                                                                                            gTGGAATGCGCGTAGCAAGGTCAAAAAGGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTc                                                                 >  1:1130185/1‑116 (MQ=255)
                                                                                                           caggctCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTCGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTgg                       <  2:1267882/138‑1 (MQ=255)
                                                                                                           cAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTgg                       >  1:1048248/1‑143 (MQ=255)
                                                                                                           cAAGGCCAAAAACGACCTGCGTCAGGAGCAGAAAATCCTTGATACCGACCATTATGGGCTGGGGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTGCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTgg                       <  2:1109901/143‑1 (MQ=255)
                                                                                                             agctCAAACAAGCCCCGCGGCAGGAGCAGGAAATCCTTGCTACCGACCATTATGGTCTGGACCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTa                     <  1:1364808/140‑1 (MQ=255)
                                                                                                             aGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGTACCGATCCTCTGTCTGGtc                     >  1:517286/1‑142 (MQ=255)
                                                                                                                ccAAAAAAGCTCTGCGTCAGGATCAGGAAATCGTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAggg                  <  2:1125938/142‑1 (MQ=255)
                                                                                                                      aaGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATAATGGGGTGGAGCGCGTGAGAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCGAGGGACCGATCGTCTGCCTGGTAGGGccgccg            >  1:775764/1‑143 (MQ=255)
                                                                                                                       aGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTATGGAGCGCGTGAAAGATCGAATCCTAGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCgg           >  1:1029091/1‑143 (MQ=255)
                                                                                                                               gTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtaggt   <  1:90340/143‑1 (MQ=255)
                                                                                                                                tCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGTCGCGCTTGAAAGATCGAATCCTCGAGTCTCTTGCGGTTCAAAGCCGTGTCAACAAAATCCAGGGCCCGATCCTCTGCCTGGAAGGGCCGCCGGGggtaggtg  >  1:1007039/1‑142 (MQ=255)
                                                                                                                                tCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtaggta  >  1:925139/1‑143 (MQ=255)
                                                                                                                                      |                                                                                                                                         
GAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGGTA  >  NC_000913/459701‑459971

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