Predicted mutation
evidence seq id position mutation annotation gene description
RA CP000730 1,963,768 C→G V365L (GTA→CTA)  hemG ← protoporphyrinogen oxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0007301,963,7680CG97.1% 217.1 / NA 70V365L (GTA→CTA) hemGprotoporphyrinogen oxidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (32/35);  minor base A (1/1);  total (34/36)
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.99e-01

TGGCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTTCATC  >  CP000730/1963625‑1963912
                                                                                                                                               |                                                                                                                                                
tgGCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATAc                                                                                                                                              >  2:159458/1‑148 (MQ=255)
 ggCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACg                                                                                                                                             <  1:107482/148‑1 (MQ=255)
 ggCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACg                                                                                                                                             >  1:160370/1‑148 (MQ=255)
   cATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTGAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCa                                                                                                                                           >  1:379971/1‑148 (MQ=255)
         tttCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAATTACGCTCGAATc                                                                                                                                     >  2:491916/1‑148 (MQ=255)
            cGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATCATCTGTATGATCATCAACCACAGTATAACCTGGTTTACCTAGATACTCACGAGTCaaa                                                                                                                                  >  2:596429/1‑148 (MQ=255)
               cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacc                                                                                                                               >  1:478978/1‑148 (MQ=255)
               cAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAAGGTGATCATTTGACTTAAATCTCTACGGACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAacc                                                                                                                               <  1:243938/148‑1 (MQ=255)
                       tGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTc                                                                                                                       >  2:588201/1‑148 (MQ=255)
                        gACCATTGTAAATTCAGGATCACCTTTAAATGTCATCATGTGACTTAAATCTCTACGTCCAATCGATACTAATTCATCATCTGTACGATCATCAACCACAGTATCACCTGGTTTACCCAGATACGCACGAATCAAAACCTTACCTTcc                                                                                                                      <  2:162906/148‑1 (MQ=255)
                               ggaaaTTCAGGATCACCTTTAAATTTCATCATTTGACTTAAATCTCTACGTACACTCGATACTAATTCATTATCTGTATGATCATCAACCCCAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                               <  2:268232/146‑1 (MQ=255)
                               gaaaaTGCAGGTTCACCTTGAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACCAATTCATTATCTGTATGATCATCACCCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                               <  2:285505/146‑1 (MQ=255)
                               gTAAATTCAGGATCACCTTTAAATGTCATGATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                               >  1:325058/1‑148 (MQ=255)
                               gTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTgtag                                                                                                               <  1:202873/148‑1 (MQ=255)
                                 aaaTTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGGTACTAATTCATTATATGTATGATCATCAACCACAGTATCACCTGGTTTACTTAGATACCCACGAATCAAAACCTTACCTACCGgtgatgta                                                                                                             >  1:17402/1‑143 (MQ=255)
                                                tttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                              <  1:312442/148‑1 (MQ=255)
                                                tttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                              >  2:251204/1‑148 (MQ=255)
                                                tttAAATGTCATAATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTc                                                                                              <  2:514611/148‑1 (MQ=255)
                                                 ttgagtggcatgatTTGACTTAAGTCTCTGCGTACAATCGATACTAATTCATTATCTGGATGATCAGCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAGCCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                             <  2:31229/140‑1 (MQ=255)
                                                 ttAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCg                                                                                             <  2:31216/148‑1 (MQ=255)
                                                     aTGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGt                                                                                         >  1:5443/1‑148 (MQ=255)
                                                      tGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTc                                                                                        <  1:499120/148‑1 (MQ=255)
                                                         catcctTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTCTTCGATGTCCaa                                                                                     <  2:546870/148‑1 (MQ=255)
                                                         catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTTCGATGTCCaa                                                                                     <  1:347979/148‑1 (MQ=255)
                                                         catcatTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                                     <  2:278177/148‑1 (MQ=255)
                                                         catcatTTGACTTAAACCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCaa                                                                                     <  1:440851/148‑1 (MQ=255)
                                                               ttGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACAt                                                                               >  1:521213/1‑148 (MQ=255)
                                                                   ctttattctcTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                           <  1:422625/142‑1 (MQ=255)
                                                                   cTTAACTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                           <  2:177367/148‑1 (MQ=255)
                                                                   cTTAAATCTCTACGTACAATCGGTACTAATTCATTTTCTTTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                           <  2:503721/148‑1 (MQ=255)
                                                                   cTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTCTCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGTCCATTTTTTCGATGTCCAAGTACATGCgg                                                                           <  1:325165/148‑1 (MQ=255)
                                                                   cTTAAATCGCTACTTACAAGCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCgg                                                                           <  1:480159/148‑1 (MQ=255)
                                                                      aaaTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTaa                                                                        <  1:461602/148‑1 (MQ=255)
                                                                        aTCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCAGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTACAAGTACATGCGGTAATg                                                                      >  2:296350/1‑148 (MQ=255)
                                                                              aCGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                                <  2:160370/148‑1 (MQ=255)
                                                                              aCGTACAATCGAGACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTa                                                                >  1:428849/1‑148 (MQ=255)
                                                                                 tACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCACGTACATGCGGTAATGTCTGTATCa                                                             >  2:266551/1‑148 (MQ=255)
                                                                                  aCAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCAc                                                            >  1:363210/1‑148 (MQ=255)
                                                                                   cactCGATACTAATTCATTATCGTTATGATACTCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACt                                                           <  2:353942/145‑1 (MQ=255)
                                                                                          tACTAATTCATTATCTGGATGATCATCAACCACAGTATCACCTGGTTTACCTAGGTACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTtctc                                                    <  2:68466/148‑1 (MQ=255)
                                                                                              aaTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTACTGTCTGTATCACTCGTTCTCGCaa                                                >  2:341999/1‑148 (MQ=255)
                                                                                              aaTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCaa                                                >  1:14777/1‑148 (MQ=255)
                                                                                              aaTTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCAAAGTACATGCGGTAATGTCTGTATCATTCGTTCTCGCaa                                                >  2:76224/1‑148 (MQ=255)
                                                                                                   aTTCTCTGTATGATCATCAACCACGGTAGAACCTGGTGGACCTAGATACGCACGAATCAAAACCTTACCTTGCGGTCTAGTAAATGGCCATTTTTTCCATGTCCAACTACATGCGGTAATGTCTGTATCACTCGTTGTTGCCATTAca                                           >  2:163892/1‑147 (MQ=255)
                                                                                                    ttATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGa                                          >  1:107945/1‑148 (MQ=255)
                                                                                                      accTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACCTGCATAAAGGTCTGTATCACTCGTTCTCGCAATTACGAAg                                        <  1:625055/146‑1 (MQ=255)
                                                                                                       tCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGc                                       >  1:456181/1‑148 (MQ=255)
                                                                                                        cTGTATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAACTTGCAATTTAACCCAAATCTAATTACATGCGGGATTGTCTGTAACACTCGTTCTCGCAATTACGAAGcc                                      >  1:118650/1‑148 (MQ=255)
                                                                                                         tGTATGATCATCAACCACAGTATCACCCGGCGAACGTATATACGCACGAAACAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCa                                     >  2:289397/1‑148 (MQ=255)
                                                                                                           tATGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTCTCACTCGTTCTCGCAATTACGAAGCCAGt                                   >  1:262819/1‑148 (MQ=255)
                                                                                                             tGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTat                                 >  2:146499/1‑147 (MQ=255)
                                                                                                             tGATCATCAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTGTTCTTGTTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTAc                                 >  1:145989/1‑148 (MQ=255)
                                                                                                               atcatcAACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCGATTACGAAGCCAGTACca                               >  1:46704/1‑148 (MQ=255)
                                                                                                                tcatcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcat                              <  1:492857/148‑1 (MQ=255)
                                                                                                                   tcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatctt                           >  1:631800/1‑146 (MQ=255)
                                                                                                                   tcaACCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACcatcat                           >  2:108671/1‑148 (MQ=255)
                                                                                                                      aCCACAGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATggg                        <  1:137836/148‑1 (MQ=255)
                                                                                                                         acaGTATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCCTGGGTAt                     >  2:58752/1‑148 (MQ=255)
                                                                                                                            gTATCTCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAtttt                  <  2:99588/148‑1 (MQ=255)
                                                                                                                            gTATCACCTGGTTTACCGAGATCCGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTAtttt                  <  2:412926/148‑1 (MQ=255)
                                                                                                                             tATCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGAGTAGTAAATGGCCATTTTTTCGATGTCCAAGTAAATGCGAAAATGTCTGTATCTCTCGTTTTGGAAATTACGGTGTCTGTACCTTCTTGGGTTTTTTc                 <  2:17402/148‑1 (MQ=255)
                                                                                                                               tccccTGGTTTACCTAAATACGCACGAATCAAACTCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCaa               <  1:115509/145‑1 (MQ=255)
                                                                                                                               tCACCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGTTTTTTTTTTACATGCGAGAATGTCTGTATCACTCGTTCTCGCAATTACGAGGCCACTACCATCATTGGTATTTTCaa               >  2:94557/1‑148 (MQ=255)
                                                                                                                                 aCCTGGTTTACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTTTTTTTCTTTGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCACTTACGAAGCCAGTACCATCATGGGTATGTCCaact             >  2:276213/1‑146 (MQ=255)
                                                                                                                                        ttACCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCttttt      <  2:274310/148‑1 (MQ=255)
                                                                                                                                          aCCTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTtca    <  1:67964/148‑1 (MQ=255)
                                                                                                                                           ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTg                                                                  <  2:7290/85‑1 (MQ=255)
                                                                                                                                           ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTg                                                                  >  1:7290/1‑85 (MQ=255)
                                                                                                                                           ccTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTtcat   <  2:521213/148‑1 (MQ=255)
                                                                                                                                            cTAGATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTtcatc  <  2:267752/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                
TGGCATACTTTTCGGCAAACGATTGACAATTGTAAATTCAGGATCACCTTTAAATGTCATCATTTGACTTAAATCTCTACGTACAATCGATACTAATTCATTATCTGTATGATCATCAACCACAGTATCACCTGGTTTACCTACATACGCACGAATCAAAACCTTACCTTCCGGTGTAGTAAATGGCCATTTTTTCGATGTCCAAGTACATGCGGTAATGTCTGTATCACTCGTTCTCGCAATTACGAAGCCAGTACCATCATGGGTATTTTCAATGTCTTTTTCATC  >  CP000730/1963625‑1963912

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