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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913459,8340CA99.3% 201.2 / ‑6.5 70A316E (GCG→GAG) lonLon protease
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (36/33);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGG  >  NC_000913/459695‑459969
                                                                                                                                           |                                                                                                                                       
aaaGCGGAAGCAGAGGTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCa                                                                                                                                      <  1:674811/143‑1 (MQ=255)
    cGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGGCAGAAAAGACCGGCGTCAGGAGCAAGac                                                                                                                                  >  1:692998/1‑142 (MQ=255)
    cGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGGTATAACGACGGGATGGTAGAGGTGACGAGGAATGCGCGTAGCAACGACAAAAAAGACACGCGTCAGGAGCAGGac                                                                                                                                  >  1:407054/1‑142 (MQ=255)
    cGGAAGCAAAGTTGCAGAAGCTGAAAAGGAGGGCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGaa                                                                                                                                  <  2:229779/143‑1 (MQ=255)
       aaGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCCGGAGCAGGACATc                                                                                                                               >  1:1055175/1‑143 (MQ=255)
        aGCAGAGTTGCATAAGCTGAAAATGGTGGCTCCGGTGGCGGCAGAAGCGGCCGGAGTGGGTGGTTATATCGACTTGATTGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCTCCTGACGCAggtctttt                                                                                                                              <  2:887230/143‑7 (MQ=255)
          gagagTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGATACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCtt                                                                                                                            <  1:976912/142‑1 (MQ=255)
               ttGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATAc                                                                                                                       >  2:477131/1‑143 (MQ=255)
                   aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCt                                                                                                                             >  2:677642/1‑133 (MQ=255)
                   aGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCt                                                                                                                             <  1:677642/133‑1 (MQ=255)
                          gAAGATGATGTCCTCGATGTCGGCAGAAGCGCCCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATg                                                                                                            <  1:839916/143‑1 (MQ=255)
                          gAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATg                                                                                                            <  1:578810/143‑1 (MQ=255)
                           aaaaTGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATgg                                                                                                           >  2:416431/1‑143 (MQ=255)
                               tgatgTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGGTACCGACCATTATGGTCTg                                                                                                       >  1:701316/1‑143 (MQ=255)
                               tgatgTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTg                                                                                                       <  1:595792/143‑1 (MQ=255)
                                  tgTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAg                                                                                                    <  2:395572/143‑1 (MQ=255)
                                      tcCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGcgcg                                                                                                >  1:281221/1‑143 (MQ=255)
                                             tCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAGGa                                                                                         >  1:281622/1‑143 (MQ=255)
                                             tCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACAGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGa                                                                                         >  1:728036/1‑143 (MQ=255)
                                             tCGGCAGAAGCGACCGTAGTGCGGGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGCTACCGCCCATTGTGGTCCGGAGCGCGTGAAAGa                                                                                         >  1:233729/1‑143 (MQ=255)
                                               ggCCGAAGCGACCTGAGTGCGTTGTTGTATCGGCTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                       <  2:937539/143‑1 (MQ=255)
                                               ggCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATc                                                                                       <  1:1045516/143‑1 (MQ=255)
                                                    aaGCGACCGTAGTGCGTGGTTATAGCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATc                                                                                  <  1:843664/143‑1 (MQ=255)
                                                     aGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATcc                                                                                 <  2:896732/143‑1 (MQ=255)
                                                     aGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATcc                                                                                 >  1:388292/1‑143 (MQ=255)
                                                         gccGTAGTGCGGTGTTATAGCGAATGGACGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGa                                                                             <  2:868140/142‑1 (MQ=255)
                                                              attGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATc                                                                        <  1:1133533/141‑1 (MQ=255)
                                                                tGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCtt                                                                      >  1:540204/1‑143 (MQ=255)
                                                                 gCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTACGGTCTGGAGCGCGCGAACGAGCGCATCCTTGAgcctactg                                                                     >  2:371851/1‑136 (MQ=255)
                                                                  cGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGc                                                                    >  1:604005/1‑143 (MQ=255)
                                                                  cGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGGGAAAGCTCGCATCCTCGAGTATCTTGc                                                                    >  1:831853/1‑143 (MQ=255)
                                                                   gTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCg                                                                   <  1:532653/143‑1 (MQ=255)
                                                                    tGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCgg                                                                  <  1:477131/143‑1 (MQ=255)
                                                                       ttATATGGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCCTTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTc                                                               <  1:390868/143‑1 (MQ=255)
                                                                          gatCGACTGTAGTGTACAGGGGCCGTGAAATGGGCTTACCAAGGTCAAAAAAGACCTGCGTTAGGAGCAGGAATGCCTTGATACCGACTATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCaaa                                                            <  2:434144/142‑1 (MQ=255)
                                                                            tCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGCGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGATAAAAGc                                                          >  1:804119/1‑143 (MQ=255)
                                                                               aCTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCgt                                                       >  1:1113149/1‑143 (MQ=255)
                                                                                  ggATGGTACGGGTGGGGGGGAAGGCGCGGAGCGAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAGATCCTGGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAAGCCTTGAGTATCTTGCGGTTCAAAGCCGTGTc                                                    <  2:260464/143‑1 (MQ=255)
                                                                                   gATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTca                                                   <  1:999190/143‑1 (MQ=255)
                                                                                       gTACAGGTGCCGGGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTcaacaa                                               <  2:426917/143‑1 (MQ=255)
                                                                                        tACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACaaa                                              >  1:716747/1‑143 (MQ=255)
                                                                                        tACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACaaa                                              <  1:416431/143‑1 (MQ=255)
                                                                                          cAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGCTACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCGTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAAt                                            >  1:61945/1‑143 (MQ=255)
                                                                                              ggCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAg                                        <  2:944424/142‑1 (MQ=255)
                                                                                                  gTGGAATGCGCGTAGCAAGGTCAAAAAAGACCAGCGTCAGGAGCAGGAAATCCTTGATACCGGCCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTCTCTTGCGGTTCCCAGCCGTGTCAACCAAATCAAAGGAc                                    >  1:889589/1‑143 (MQ=255)
                                                                                                       aTGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGAGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAACGCCGTGTCAACAAAATCAAGGGACCGATc                               >  2:410020/1‑143 (MQ=255)
                                                                                                          cgcgTAGCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGAGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCctc                            >  2:727409/1‑143 (MQ=255)
                                                                                                             ttaccAAGGTTAAAATATACCTGCGTTAGGAACAGTAAATCCTTTTTACCTATTATTATTGTCTAGAGCGCGGGAAAGAACGAAACCGTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGc                         <  1:657808/139‑1 (MQ=255)
                                                                                                               aGCAAGGTGAAGAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCt                       <  1:921402/143‑1 (MQ=255)
                                                                                                                gCAAGGTCAAAAAAGGCCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTg                      <  1:1081043/143‑1 (MQ=255)
                                                                                                                gCAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGGTACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTg                      <  2:701316/143‑1 (MQ=255)
                                                                                                                 cAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACCAAATCAAGGGACCGACCCCATGCCCgg                     >  2:910747/1‑143 (MQ=255)
                                                                                                                 cAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTgg                     >  1:92148/1‑143 (MQ=255)
                                                                                                                 cAAGGTCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTgg                     >  1:241793/1‑143 (MQ=255)
                                                                                                                     atCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAgg                 <  1:870215/142‑1 (MQ=255)
                                                                                                                      tCAAAAAAGACCTGCGTCAGGAGCAGGAAATCCGTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAggg                <  2:585411/143‑1 (MQ=255)
                                                                                                                       cAAAAAAGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGCATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGgc               >  1:378808/1‑143 (MQ=255)
                                                                                                                        aaaaaaGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTCTGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAGAGCCGGGTCAAGAAAATCACGGAACCGATACTCTGCATGGGTGGgcg              >  2:1014113/1‑142 (MQ=255)
                                                                                                                         aaaaaGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGgccg             >  2:814534/1‑143 (MQ=255)
                                                                                                                         aaaaaGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAACGGACCGATCCTCTGCCTGGTAGGgccg             >  2:545029/1‑143 (MQ=255)
                                                                                                                           aaaGACCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGgccgcc           >  2:440415/1‑143 (MQ=255)
                                                                                                                           aaaGACCTGCGTCAGGAGCAGGAAATCCGTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAACATCAAGGGACAGATCCGCTGCCTGGTACGgccgcc           >  2:688341/1‑143 (MQ=255)
                                                                                                                             aGACCTGCGCCAGGAGCAGGAACTCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCgg         <  1:630024/143‑1 (MQ=255)
                                                                                                                               gccTGCGACAGGAGCAGGTGTTGGTTGATCGCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTAGCGTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCgggg       <  1:414299/142‑1 (MQ=255)
                                                                                                                               aCCTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAAGCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCgggg       >  1:767108/1‑143 (MQ=255)
                                                                                                                                ccTGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCggggg      >  1:930050/1‑143 (MQ=255)
                                                                                                                                  tGCGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggta    >  1:992080/1‑143 (MQ=255)
                                                                                                                                  tGCGTCAGCAGCAGGCAATCCTGGCTACGGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggta    <  1:724190/143‑1 (MQ=255)
                                                                                                                                   gCGTCAGGAGCTGGAAATCGTTGATGCCGACCATGTTGGTCTGGAGGGCGTGAAAGATCGAGTCCTTGAGTAGCTTGGGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCAATCCTCTGCCTGGTAGGGCCGCCGGGggtag   <  2:971631/143‑1 (MQ=255)
                                                                                                                                    cGTCAGGAGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGggtagg  >  1:17059/1‑143 (MQ=255)
                                                                                                                                           |                                                                                                                                       
AAAGCGGAAGCAGAGTTGCAGAAGCTGAAAATGATGTCTCCGATGTCGGCAGAAGCGACCGTAGTGCGTGGTTATATCGACTGGATGGTACAGGTGCCGTGGAATGCGCGTAGCAAGGTCAAAAAAGACCTGCGTCAGGCGCAGGAAATCCTTGATACCGACCATTATGGTCTGGAGCGCGTGAAAGATCGAATCCTTGAGTATCTTGCGGTTCAAAGCCGTGTCAACAAAATCAAGGGACCGATCCTCTGCCTGGTAGGGCCGCCGGGGGTAGG  >  NC_000913/459695‑459969

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: