Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,683,566 (C)6→7 coding (1047/1404 nt) fumC ← fumarase C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,683,5601.C91.6% 698.4 / 55.3 227coding (1053/1404 nt)fumCfumarase C
Reads supporting (aligned to +/- strand):  ref base . (10/9);  new base C (104/104);  total (114/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 = 2.49e-01

Only 100 of 255 total aligned reads displayed.

TTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCA  >  NC_000913/1683469‑1683651
                                                                                            |                                                                                              
  tCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑ccc                                                                                           >  2:8371227/1‑93 (MQ=255)
   ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑cccc                                                                                         >  1:2173403/1‑93 (MQ=255)
   ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑cccc                                                                                         <  1:977361/93‑1 (MQ=255)
    cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGCACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑ccccc                                                                                        >  1:6328874/1‑93 (MQ=255)
      tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑ccccccc                                                                                      <  2:2057362/93‑2 (MQ=255)
       gCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑ccccccca                                                                                     >  1:210275/1‑91 (MQ=255)
       gCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCCCAt                                                                                     <  2:4662775/93‑1 (MQ=255)
        ccATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑cccccccat                                                                                    >  2:3754929/1‑90 (MQ=255)
         cATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑cccccccatg                                                                                   >  2:8594563/1‑89 (MQ=255)
           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGtt                                                                                  >  1:2521173/1‑93 (MQ=255)
           tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGtt                                                                                  >  1:3040446/1‑93 (MQ=255)
              gCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGAt                                                                               <  2:6230800/93‑1 (MQ=255)
              gCCAGCAAGCGCACCGATTGCAGGAAATTGGTGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGAt                                                                               <  2:2112707/93‑1 (MQ=255)
               ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATc                                                                              >  2:5137565/1‑93 (MQ=255)
                cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCg                                                                             <  2:4836636/93‑1 (MQ=255)
                 aGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGc                                                                            >  1:1223202/1‑93 (MQ=255)
                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGcc                                                                           >  2:1921935/1‑93 (MQ=255)
                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGcc                                                                           <  2:5932539/93‑1 (MQ=255)
                  gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGcc                                                                           <  2:6416193/93‑1 (MQ=255)
                   cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCa                                                                          <  1:3279186/93‑1 (MQ=255)
                   cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCa                                                                          >  1:7422785/1‑93 (MQ=255)
                   cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCCCATGTTGATCGCCAc                                                                         <  1:402661/93‑1 (MQ=255)
                     aGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcg                                                                        >  2:3312110/1‑93 (MQ=255)
                      gcgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgt                                                                       >  2:7607495/1‑93 (MQ=255)
                      gcgcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgt                                                                       <  1:7661502/93‑1 (MQ=255)
                        gcACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgtcg                                                                     <  1:8103218/93‑1 (MQ=255)
                         cACCGATTGCGGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgtcgt                                                                    <  1:4716725/93‑1 (MQ=255)
                         cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgtcgt                                                                    >  1:5354348/1‑93 (MQ=255)
                         cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgtcgt                                                                    >  1:4753069/1‑93 (MQ=255)
                         cACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAcgtcgt                                                                    >  1:5528162/1‑93 (MQ=255)
                          aCCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGtt                                                                   >  2:3561777/1‑93 (MQ=255)
                            cGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTcc                                                                 >  2:5526874/1‑93 (MQ=255)
                             gATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTccc                                                                <  1:8578034/93‑1 (MQ=255)
                              aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTcccc                                                               >  2:116419/1‑93 (MQ=255)
                              aTTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTcccc                                                               <  2:5104829/93‑1 (MQ=255)
                                  cAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCCATGTTGATCGCCACGTCGTTCCCCATc                                                            >  2:7926199/1‑93 (MQ=255)
                                    ggAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCAcc                                                         <  2:1691190/93‑1 (MQ=255)
                                    ggAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCAcc                                                         >  2:698799/1‑93 (MQ=255)
                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                        <  2:23724/93‑1 (MQ=255)
                                     gAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                        >  1:3136756/1‑93 (MQ=255)
                                       aaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                      >  2:5805797/1‑93 (MQ=255)
                                       aaTTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGa                                                      >  1:2545664/1‑93 (MQ=255)
                                         ttGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAca                                                    >  1:3688036/1‑93 (MQ=255)
                                          tgtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcag                                                   >  1:8809204/1‑93 (MQ=255)
                                           gtgtATCCCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                  <  2:1602515/93‑1 (MQ=255)
                                           gtgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagc                                                  >  2:4397517/1‑93 (MQ=255)
                                            tgGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                 >  1:5458915/1‑93 (MQ=255)
                                              gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                               <  1:8708872/93‑1 (MQ=255)
                                              gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCaga                                               <  2:4222145/93‑1 (MQ=255)
                                               gtCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCagag                                              <  2:2595907/92‑1 (MQ=255)
                                               aTCACCATTGGACGGAAGACGTTCAGTTCCAAGTTACCGGAAGCGCCCCCCCATGTTTATCGCCACGTCGTTCCCCATCACCTGCCCGCagag                                              >  1:5946945/1‑93 (MQ=255)
                                               aTCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                             >  1:3914531/1‑93 (MQ=255)
                                                 cACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCa                                            >  2:7401026/1‑93 (MQ=255)
                                                  aCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                           >  1:4479715/1‑93 (MQ=255)
                                                  aCCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCAt                                           >  2:676670/1‑93 (MQ=255)
                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                          >  2:3788866/1‑93 (MQ=255)
                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                          <  1:5726981/93‑1 (MQ=255)
                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                          <  1:51023/93‑1 (MQ=255)
                                                   ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATg                                          <  1:3218708/93‑1 (MQ=255)
                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                         >  1:206086/1‑93 (MQ=255)
                                                    cATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                         >  2:2854238/1‑93 (MQ=255)
                                                     aTTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                        >  1:1651201/1‑93 (MQ=255)
                                                     aTTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                        <  2:364569/93‑1 (MQ=255)
                                                     aTTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                        >  2:7748311/1‑93 (MQ=255)
                                                      ttGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                       <  2:5595458/93‑1 (MQ=255)
                                                      ttGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGtt                                       <  1:644980/93‑1 (MQ=255)
                                                            ggAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                 <  2:4590079/93‑1 (MQ=255)
                                                             gAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                >  1:2872117/1‑93 (MQ=255)
                                                             gAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                >  2:4528574/1‑93 (MQ=255)
                                                               acaCGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTca                              >  2:1169952/3‑93 (MQ=255)
                                                               aGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCAGGTCGTTCCCCATCACCTGCCAGCAGAGCATGGTTAATGCCTca                              >  1:1542814/1‑93 (MQ=255)
                                                                  cGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCacac                           <  1:3568449/93‑1 (MQ=255)
                                                                  cGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCacac                           <  1:1022727/93‑1 (MQ=255)
                                                                  cGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCacac                           <  1:6111422/93‑1 (MQ=255)
                                                                  cGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                          >  2:6558657/1‑93 (MQ=255)
                                                                        gttcaAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGtt                     <  2:3110621/93‑1 (MQ=255)
                                                                          tcaAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                   >  2:1486288/1‑93 (MQ=255)
                                                                          tcaAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTgg                   >  1:7538188/1‑93 (MQ=255)
                                                                           caAAGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTggg                  <  2:4722073/93‑1 (MQ=255)
                                                                            aaaGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt                 <  1:1874931/93‑1 (MQ=255)
                                                                            aaaGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt                 >  2:823186/1‑93 (MQ=255)
                                                                            aaaGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt                 >  2:7672310/1‑93 (MQ=255)
                                                                             aaGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                <  2:4951038/93‑1 (MQ=255)
                                                                             aaGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                <  2:5853411/93‑1 (MQ=255)
                                                                             aaGTTACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGtt                <  1:2956624/93‑1 (MQ=255)
                                                                                ttACCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc             <  1:6284254/93‑1 (MQ=255)
                                                                                ttAACGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCAc             <  2:8753720/93‑1 (MQ=255)
                                                                                  aCCGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACtt           >  1:6739741/1‑93 (MQ=255)
                                                                                  aCCGGAAGCG‑‑CCCACCCATGTTTATCGCCACGTCGGTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACtt           <  2:6334444/93‑1 (MQ=255)
                                                                                   ccGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACttt          >  1:6558789/1‑93 (MQ=255)
                                                                                    cGGAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc         <  2:1835397/93‑1 (MQ=255)
                                                                                     ggAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcc        >  1:3272809/1‑93 (MQ=255)
                                                                                     ggAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcc        <  2:5358566/93‑1 (MQ=255)
                                                                                      gAAGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTccc       <  1:8462733/93‑1 (MQ=255)
                                                                                       aaGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc      >  1:467557/1‑93 (MQ=255)
                                                                                       aaGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc      <  1:4139780/93‑1 (MQ=255)
                                                                                       aaGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTcccc      <  2:8042072/93‑1 (MQ=255)
                                                                                       aaGCGCCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTAAATGCCTCACACTGTGTTGGGTTCACTTTcccc      >  1:7031949/1‑93 (MQ=255)
                                                                                        aGCG‑‑CCACCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTCGGTTCACTTTCCCCg     >  2:8554100/1‑93 (MQ=255)
                                                                                          cg‑‑cCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCa  <  2:41160/93‑1 (MQ=255)
                                                                                            |                                                                                              
TTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCG‑‑CCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCA  >  NC_000913/1683469‑1683651

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