Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,162,811 A→G 9.6% I224T (ATA→ACA)  plsC ← 1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,162,8110AG 9.6% 179.4 / 7.1 62I224T (ATA→ACA) plsC1‑acyl‑sn‑glycerol‑3‑phosphate acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base G (3/3);  total (27/35)
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.86e-01

GACATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGG  >  NC_000913/3162665‑3162956
                                                                                                                                                  |                                                                                                                                                 
gACATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATc                                                                                                                                                 <  2:334696/149‑1 (MQ=255)
gACATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATc                                                                                                                                                 <  1:574666/149‑1 (MQ=255)
 aCATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCg                                                                                                                                                >  4:404413/1‑149 (MQ=255)
 aCATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCg                                                                                                                                                >  2:394416/1‑149 (MQ=255)
   atatTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGaa                                                                                                                                              >  2:554222/1‑149 (MQ=255)
   atatTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGaa                                                                                                                                              <  1:11000/149‑1 (MQ=255)
      tttGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGa                                                                                                                                           >  2:565343/1‑149 (MQ=255)
          cTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAAt                                                                                                                                       <  2:55266/149‑1 (MQ=255)
                 tAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGcagc                                                                                                                                <  2:273833/149‑1 (MQ=255)
                 tAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGcagc                                                                                                                                <  2:567331/149‑1 (MQ=255)
                        ggCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCa                                                                                                                         <  3:79997/149‑1 (MQ=255)
                                   tACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATc                                                                                                              <  1:76979/149‑1 (MQ=255)
                                      cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGCTTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                           <  1:510387/149‑1 (MQ=255)
                                      cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                           <  4:264403/149‑1 (MQ=255)
                                      cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                           <  1:554222/149‑1 (MQ=255)
                                      cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCttt                                                                                                           <  1:308347/149‑1 (MQ=255)
                                      cGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTCCCATTATCAAACGACAATGGGCAGCCAGCTCACGAACCTGATCgtt                                                                                                           <  4:340391/149‑1 (MQ=255)
                                          cAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCa                                                                                                       <  2:253711/149‑1 (MQ=255)
                                           aGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCAt                                                                                                      >  2:50725/1‑149 (MQ=255)
                                                 gtggTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTAATCGAGCTCGGCGATTGTTTGTTCCATTATCGAACGACAATGGGCAGCGCGCTCACGAACCTGATCTTTGCCATactgac                                                                                                >  4:231331/1‑149 (MQ=255)
                                                    ttCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATActgactga                                                                                             <  3:170234/149‑1 (MQ=255)
                                                      ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                  >  4:285702/1‑107 (MQ=255)
                                                      ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                  <  3:285702/110‑4 (MQ=255)
                                                      ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                  <  3:561599/110‑4 (MQ=255)
                                                      ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                  >  4:561599/1‑107 (MQ=255)
                                                      ccccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                  <  1:241063/110‑4 (MQ=255)
                                                      acccGCAAAAAGACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTGTCGTTCGATAATggaac                                                                                                                                  >  2:241063/2‑107 (MQ=255)
                                                       cccGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGt                                                                                          >  4:53956/1‑149 (MQ=255)
                                                          gCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCaa                                                                                       >  1:241745/1‑149 (MQ=255)
                                                          gCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCaa                                                                                       >  1:578336/1‑149 (MQ=255)
                                                            aaaaaTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAAtt                                                                                     <  3:404413/149‑1 (MQ=255)
                                                               aaTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggc                                                                                  <  2:346648/149‑1 (MQ=255)
                                                                aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                 <  3:272351/149‑1 (MQ=255)
                                                                aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                 >  3:38670/1‑149 (MQ=255)
                                                                aTACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTggcg                                                                                 >  1:395015/1‑149 (MQ=255)
                                                                   cAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGgca                                                                              <  2:179317/149‑1 (MQ=255)
                                                                       ttGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGACCGACAATGGGCAGCCAGCTCACGAGCCTGATCTTTGCCATACTGACTGCCGCCAATTGACGtcagtca                                                                          <  4:132561/149‑4 (MQ=255)
                                                                           cGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCg                                                                      >  2:574690/1‑149 (MQ=255)
                                                                           cGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCg                                                                      >  3:138690/1‑149 (MQ=255)
                                                                              cTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACa                                                                   <  4:169222/149‑1 (MQ=255)
                                                                                    cTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCAcc                                                             >  1:75573/1‑149 (MQ=255)
                                                                                        tCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGAc                                                         <  2:226338/149‑1 (MQ=255)
                                                                                         ccGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGAcc                                                        >  2:417920/1‑149 (MQ=255)
                                                                                         ccGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGAcc                                                        >  2:409456/1‑149 (MQ=255)
                                                                                              ggcTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTtgt                                                   >  4:231421/1‑149 (MQ=255)
                                                                                                    gcgTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCTATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTc                                             <  1:477603/149‑1 (MQ=255)
                                                                                                    gcgTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTc                                             >  1:105510/1‑149 (MQ=255)
                                                                                                      gTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCga                                           <  4:209867/149‑1 (MQ=255)
                                                                                                       ttCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgat                                          <  3:53956/149‑1 (MQ=255)
                                                                                                       ttCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgat                                          <  1:565343/149‑1 (MQ=255)
                                                                                                               cttcttTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaagatta                                  >  2:334806/1‑149 (MQ=255)
                                                                                                                 tcttTATCGAGCTCGGCGATTTTTTGATCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAAt                                <  3:299062/149‑1 (MQ=255)
                                                                                                                   tttATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaaga                                     >  3:325528/1‑142 (MQ=255)
                                                                                                                   tttATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCgattaaga                                     <  4:325528/142‑1 (MQ=255)
                                                                                                                     tATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGgca                                                                              <  1:157406/99‑1 (MQ=255)
                                                                                                                     tATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGgca                                                                              >  2:157406/1‑99 (MQ=255)
                                                                                                                      aTCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTAt                           <  3:388262/149‑1 (MQ=255)
                                                                                                                       tCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTAtt                          <  4:38670/149‑1 (MQ=255)
                                                                                                                             tCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGt                    >  1:372492/1‑149 (MQ=255)
                                                                                                                                              cATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACgg   >  3:478083/1‑149 (MQ=255)
                                                                                                                                              cATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCAGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACgg   <  1:105186/149‑1 (MQ=255)
                                                                                                                                               aTTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACggg  >  4:36112/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GACATATTTGCTCCATGTAAAACTGGCTAAAGTAATACCGTGCAGGAAAGTGTTCCCCGCAAAAATACAGATTGCCGACTTAAACTTTTCCGGCGGCTTCGCGTTCTGCGACTTCTTTATCGAGCTCGGCGATTTTTTGTTCCATTATCGAACGACAATGGGCAGCCAGCTCACGAACCTGATCTTTGCCATACTGACTGACGTCAATTGGCGGCAGCATTTCGACAATCACCAGACCGTTGTGCAGTCGATTAAGATTAATCTTATTCGAAGTTGTAGAGACGCACACGGG  >  NC_000913/3162665‑3162956

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