Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,685,542 (C)6→5 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,685,5370C.98.4% 258.5 / ‑2.7 62coding (1052/1404 nt)fumCfumarase C
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (29/32);  minor base A (1/0);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685399‑1685668
                                                                                                                                          |                                                                                                                                   
cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                    <  2:408411/140‑1 (MQ=255)
cGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGcc                                                                                                                                    <  1:164163/140‑1 (MQ=255)
 gCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGTGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGccc                                                                                                                                   <  1:821773/140‑1 (MQ=255)
            ggTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATc                                                                                                                       >  1:1337443/1‑140 (MQ=255)
                   tACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCAcgt                                                                                                                >  1:1319689/1‑140 (MQ=255)
                        aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTc                                                                                                           >  2:764807/1‑140 (MQ=255)
                        aCTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTc                                                                                                           >  2:1049947/1‑140 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         <  2:440684/140‑1 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         >  1:915565/1‑140 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         <  1:54321/140‑1 (MQ=255)
                          tGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTccc                                                                                                         >  1:440684/1‑140 (MQ=255)
                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                        >  1:1071641/1‑140 (MQ=255)
                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                        >  2:388561/1‑140 (MQ=255)
                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                        >  2:318892/1‑140 (MQ=255)
                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                        >  2:116407/1‑140 (MQ=255)
                           gcgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTcccc                                                                                                        >  1:820158/1‑140 (MQ=255)
                            cgcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCa                                                                                                       >  1:1160089/1‑140 (MQ=255)
                             gcAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                      >  2:1215087/1‑140 (MQ=255)
                               aGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCa                                                                                                    <  2:195692/140‑1 (MQ=255)
                                  gtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCTCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                 >  2:537822/1‑140 (MQ=255)
                                  gtttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCt                                                                                                 <  2:1160089/140‑1 (MQ=255)
                                   tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                      <  2:1316282/134‑1 (MQ=255)
                                   tttgttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCAt                                                                                                      >  1:1316282/1‑134 (MQ=255)
                                      gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAca                                                                                             >  1:255208/1‑140 (MQ=255)
                                      gttAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAca                                                                                             >  1:124632/1‑140 (MQ=255)
                                         aaaaCTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGAcagca                                                                                          <  1:722180/140‑1 (MQ=255)
                                             cTTTCCATGCCATCTGCCAGCAAGCTCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                      <  1:537822/140‑1 (MQ=255)
                                             cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                      <  1:636534/140‑1 (MQ=255)
                                             cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                      <  1:764807/140‑1 (MQ=255)
                                             cTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGc                                                                                      >  2:529817/1‑140 (MQ=255)
                                                 ccATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATgg                                                                                  >  2:203712/1‑140 (MQ=255)
                                                  cATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGt                                                                                 <  2:1071641/140‑1 (MQ=255)
                                                    tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                               <  1:1170518/140‑1 (MQ=255)
                                                    tGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTa                                                                               <  1:435112/140‑1 (MQ=255)
                                                         tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                          <  1:116407/140‑1 (MQ=255)
                                                         tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                          <  1:1264714/140‑1 (MQ=255)
                                                         tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                          <  2:1319689/140‑1 (MQ=255)
                                                         tCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGcc                                                                          <  2:955286/140‑1 (MQ=255)
                                                          cTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCt                                                                         >  2:1077040/1‑140 (MQ=255)
                                                             ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                      <  1:1049947/140‑1 (MQ=255)
                                                             ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                      <  1:318892/140‑1 (MQ=255)
                                                             ccAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcac                                                                      <  1:406935/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:36611/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  2:1337443/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  2:820158/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:1400314/140‑1 (MQ=255)
                                                              cAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:529817/140‑1 (MQ=255)
                                                                gCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACt                                                                   >  2:834670/1‑140 (MQ=255)
                                                                 cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     >  2:618470/1‑137 (MQ=255)
                                                                 cAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTcaca                                                                     <  1:618470/137‑1 (MQ=255)
                                                                         ccGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGt                                                          <  2:124632/140‑1 (MQ=255)
                                                                              tGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACt                                                     <  1:1215087/140‑1 (MQ=255)
                                                                                 aGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTc                                                  <  1:1077040/140‑1 (MQ=255)
                                                                                           ggATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        <  1:388561/140‑1 (MQ=255)
                                                                                            gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        <  1:567029/139‑1 (MQ=255)
                                                                                            gATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGa                                        >  2:567029/1‑139 (MQ=255)
                                                                                                 ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGc                                  >  2:191779/1‑140 (MQ=255)
                                                                                                 ccATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGGGTTGGGTTCACTTTCCCCGGCATGATTGAGc                                  >  2:1100873/1‑140 (MQ=255)
                                                                                                            aaGACGTTCAGTTCAAAGTTACCGGAAGCACCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCa                       >  2:750913/1‑140 (MQ=255)
                                                                                                                cGTTCAGTTCAAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCGGCATGATTGAACTGCCCGGCTCATTTTc                  >  2:418428/1‑140 (MQ=255)
                                                                                                                   tcagttcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCg                <  2:709603/140‑1 (MQ=255)
                                                                                                                        tcaAAGTTACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGAtt           >  1:1286160/1‑140 (MQ=255)
                                                                                                                              ttACCGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGAtt     <  1:203712/140‑1 (MQ=255)
                                                                                                                                 ccGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  >  2:1365081/1‑140 (MQ=255)
                                                                                                                                 ccGGAAGCGCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCa  >  2:96242/1‑140 (MQ=255)
                                                                                                                                          |                                                                                                                                   
CGCTCACGATTCGGTTCAATACCCACTGCGCAGTGTTTGTTAAAACTTTCCATGCCATCTGCCAGCAAGCGCACCGATTGCAGGAAATTGTGGATCACCATTGGACGGAAGACGTTCAGTTCAAAGTTACCGGAAGCGCCCCCCATGTTGATCGCCACGTCGTTCCCCATCACCTGACAGCAGAGCATGGTTAATGCCTCACACTGTGTTGGGTTCACTTTCCCCGGCATGATTGAGCTGCCCGGCTCATTTTCCGGGATTGAGATTTCA  >  NC_000913/1685399‑1685668

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