Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,987,935 G→A R156C (CGC→TGC)  hemC ← porphobilinogen deaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,987,9350GA99.6% 825.0 / ‑5.7 230R156C (CGC→TGC) hemCporphobilinogen deaminase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (116/113);  total (117/113)
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.64e-01

Only 100 of 247 total aligned reads displayed.

GACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGT  >  NC_000913/3987845‑3988023
                                                                                          |                                                                                        
gACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGcaca                                                                                        >  1:1317882/1‑93 (MQ=255)
gACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGcaca                                                                                        >  2:3053119/1‑93 (MQ=255)
gACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGcaca                                                                                        >  1:94322/1‑93 (MQ=255)
 aCCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAg                                                                                       <  1:994715/93‑1 (MQ=255)
  ccTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAgg                                                                                      >  2:3295561/1‑93 (MQ=255)
  ccTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAgg                                                                                      >  1:1490565/1‑93 (MQ=255)
   cTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAggg                                                                                     >  2:3171635/1‑93 (MQ=255)
   cTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAggg                                                                                     <  1:4111444/93‑1 (MQ=255)
      aaCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGc                                                                                  >  1:1574036/1‑93 (MQ=255)
      aaCGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGc                                                                                  <  1:397714/93‑1 (MQ=255)
         gTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGa                                                                               >  2:206132/1‑93 (MQ=255)
          ttttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAt                                                                              >  1:2644512/1‑93 (MQ=255)
          ttttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAt                                                                              >  1:649428/1‑93 (MQ=255)
          ttttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGAt                                                                              >  1:1592189/1‑93 (MQ=255)
            ttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATaa                                                                            <  2:2191909/93‑1 (MQ=255)
            ttAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATaa                                                                            >  1:2821332/1‑93 (MQ=255)
            ttAGTCCGGCTACGGAAAGAATGATGGCTCCGTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTTCCGCACAGGGAGCGGATaa                                                                            <  2:1265361/93‑1 (MQ=255)
               gTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCa                                                                         <  2:1315385/93‑1 (MQ=255)
               gTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCa                                                                         >  2:3124675/1‑93 (MQ=255)
                 ccGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGa                                                                       <  1:1846716/93‑1 (MQ=255)
                  cGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGAt                                                                      >  2:600709/1‑93 (MQ=255)
                    gCTACAGCACGAATGAAGGCATCGTATGAGGCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGACGCCGCACAGGGAGCGGATAATCAGAttt                                                                    >  2:1520101/1‑91 (MQ=255)
                        cgtcAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGAc                                                                <  1:3907350/90‑1 (MQ=255)
                           cAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCAGGACGGc                                                             >  2:4144307/1‑93 (MQ=255)
                            aaGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGTCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCg                                                            >  1:781066/1‑93 (MQ=255)
                             aGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGt                                                           <  1:1896487/93‑1 (MQ=255)
                               aaTGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTc                                                         <  2:2534848/93‑1 (MQ=255)
                                atgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCa                                                        >  2:48744/1‑93 (MQ=255)
                                atgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCa                                                        >  1:2788660/1‑93 (MQ=255)
                                atgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGCAGCGGATAATCAGATCCGGACGGCGTTCa                                                        >  1:2737795/1‑93 (MQ=255)
                                 tgatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAg                                                       >  1:753396/1‑93 (MQ=255)
                                  gatgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGc                                                      <  2:2286245/93‑1 (MQ=255)
                                   atgGCATCTTATTCTCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGcc                                                     <  1:4441831/93‑1 (MQ=255)
                                    tgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCa                                                    <  2:701974/93‑1 (MQ=255)
                                    tgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCa                                                    <  1:102369/93‑1 (MQ=255)
                                    tgGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCa                                                    <  2:1877212/93‑1 (MQ=255)
                                     ggCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAg                                                   <  1:520199/93‑1 (MQ=255)
                                     ggCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAg                                                   >  1:2325921/1‑93 (MQ=255)
                                      gCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGt                                                  <  1:2278795/93‑1 (MQ=255)
                                      gCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGt                                                  <  1:3657324/93‑1 (MQ=255)
                                       cATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGtt                                                 >  1:827259/1‑93 (MQ=255)
                                        aTCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTg                                                <  2:2647546/93‑1 (MQ=255)
                                        aTCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTg                                                >  1:3042205/1‑93 (MQ=255)
                                        aTCGTATTAGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTg                                                >  2:4554773/1‑93 (MQ=255)
                                           gTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCa                                             >  2:1164176/1‑93 (MQ=255)
                                            tATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCAc                                            >  1:4147365/1‑93 (MQ=255)
                                               tCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTgg                                         <  1:401209/93‑1 (MQ=255)
                                                cGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGc                                        <  2:2590404/93‑1 (MQ=255)
                                                 gCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCg                                       >  2:3801074/1‑93 (MQ=255)
                                                 gCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCg                                       >  2:2452365/1‑93 (MQ=255)
                                                  ccGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGa                                      <  1:905858/93‑1 (MQ=255)
                                                  ccGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGa                                      >  1:1467920/1‑93 (MQ=255)
                                                  ccGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGa                                      <  2:1587865/93‑1 (MQ=255)
                                                     ttATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGt                                   >  2:428090/1‑93 (MQ=255)
                                                     ttATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGt                                   <  2:229037/93‑1 (MQ=255)
                                                      tATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTa                                  >  2:4420663/1‑93 (MQ=255)
                                                       aTCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaa                                 >  1:2182773/1‑93 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaaa                                >  1:199122/1‑93 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaaa                                >  1:4210807/1‑93 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTaaa                                <  2:3421832/93‑1 (MQ=255)
                                                        tCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGGTCAGCCAGTTGGTACGGGCGACGTgaa                                >  2:1366482/1‑93 (MQ=255)
                                                         ccAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAc                               <  1:583295/93‑1 (MQ=255)
                                                         ccAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAc                               <  2:79164/93‑1 (MQ=255)
                                                         ccAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAAc                               <  1:3418672/93‑1 (MQ=255)
                                                          cAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACt                              <  1:2599289/93‑1 (MQ=255)
                                                          cAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACt                              <  2:535816/93‑1 (MQ=255)
                                                          cAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACt                              <  1:391932/93‑1 (MQ=255)
                                                            gTTTGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTgg                            <  2:2334092/93‑1 (MQ=255)
                                                             tttGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGa                           <  2:2925281/93‑1 (MQ=255)
                                                             tttGCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGa                           <  2:1665371/93‑1 (MQ=255)
                                                                gCTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGt                        >  2:3702850/1‑93 (MQ=255)
                                                                 cTCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTc                       <  2:4369936/93‑1 (MQ=255)
                                                                  tCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTcc                      >  1:1664174/1‑93 (MQ=255)
                                                                  tCAGGCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTcc                      <  1:433653/93‑1 (MQ=255)
                                                                   cAGCCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTccc                     <  2:868050/93‑1 (MQ=255)
                                                                     ggCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCcga                   >  1:3178690/1‑93 (MQ=255)
                                                                     ggCGAGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTCCAGCCAGTTGGCACTGGCGACGTAAACTGGCGCTCCcga                   >  1:1188833/1‑93 (MQ=255)
                                                                         aGTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCCGCCAGTTGGCACTGGCGCAGTAAACTGGACGTCCCGAAGAt               >  2:1226187/1‑93 (MQ=255)
                                                                          gTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATa              <  1:4603985/93‑1 (MQ=255)
                                                                          gTGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATa              <  2:2079195/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  2:584302/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  2:3566438/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  1:4494088/93‑1 (MQ=255)
                                                                           tGCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAc             <  2:1846310/93‑1 (MQ=255)
                                                                            gCCGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATAct            >  1:3725963/1‑93 (MQ=255)
                                                                             ccGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActg           >  1:1105186/1‑93 (MQ=255)
                                                                              cGACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgc          <  2:3522714/93‑1 (MQ=255)
                                                                               gACGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcc         <  2:3615735/93‑1 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        >  2:3824764/1‑93 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        <  2:2172306/93‑1 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        >  2:3764094/1‑93 (MQ=255)
                                                                                aCGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcct        <  1:1367958/93‑1 (MQ=255)
                                                                                 cGTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctg       >  2:3610686/1‑93 (MQ=255)
                                                                                  gTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctgc      >  1:3292637/1‑93 (MQ=255)
                                                                                  gTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctgc      >  1:1825711/1‑93 (MQ=255)
                                                                                  gTTGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATActgcctgc      <  2:742307/93‑1 (MQ=255)
                                                                                   ttGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACtgcctgcc     <  2:1967635/93‑1 (MQ=255)
                                                                                   ttGCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACtgcctgcc     >  2:1030644/1‑93 (MQ=255)
                                                                                     gCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCgg   >  2:2849335/1‑93 (MQ=255)
                                                                                     gCCGCACAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCgg   >  2:2898493/1‑93 (MQ=255)
                                                                                          |                                                                                        
GACCTAAACGTTTTAGTCCGGCTACGGCAAGAATGATGGCATCGTATTCGCCGTTATCCAGTTTGCTCAGGCGAGTGCCGACGTTGCCGCGCAGGGAGCGGATAATCAGATCCGGACGGCGTTCAGCCAGTTGGCACTGGCGACGTAAACTGGACGTCCCGACGATACTGCCTGCCGGT  >  NC_000913/3987845‑3988023

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