Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,903,621 G→A 18.5% G27G (GGC→GGT bglB ← cryptic phospho‑beta‑glucosidase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,903,6210GA18.5% 158.4 / 29.2 65G27G (GGC→GGTbglBcryptic phospho‑beta‑glucosidase B
Reads supporting (aligned to +/- strand):  ref base G (19/34);  new base A (7/5);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCG  >  NC_000913/3903483‑3903753
                                                                                                                                          |                                                                                                                                    
ggCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGc                                                                                                                                     >  1:214943/1‑140 (MQ=255)
ggCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGc                                                                                                                                     <  1:152756/140‑1 (MQ=255)
ggCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGc                                                                                                                                     <  1:152584/140‑1 (MQ=255)
       aaCGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTc                                                                                                                              <  2:119167/140‑1 (MQ=255)
       aaCGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTc                                                                                                                              >  1:1422121/1‑140 (MQ=255)
         cgcgATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTGTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCGGAGGTCGAGATCCCTTTGCCATCTTCCt                                                                                                                            >  2:1864320/1‑140 (MQ=255)
                     cGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCtt                                                                                                                <  1:955219/140‑1 (MQ=255)
                          aaCGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCAcc                                                                                                           <  2:1778181/140‑1 (MQ=255)
                           aCGGTGATAAAAATCGATGGCGACATATTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCt                                                                                                          >  2:385662/1‑140 (MQ=255)
                            cGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTg                                                                                                         >  2:1020720/1‑140 (MQ=255)
                             ggTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGa                                                                                                        <  1:1371577/140‑1 (MQ=255)
                                 aTAAAAATCGATGGCGACATCTTTGATATTCTCTATCCCCAGGATGCGCGGTTCCATGTTTCCCATTACGCCATGAGGCAGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCAGATTgg                                                                                                    >  1:1030096/1‑140 (MQ=255)
                                       aTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                                              >  1:1248577/1‑140 (MQ=255)
                                       aTCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTg                                                                                              >  1:1721119/1‑140 (MQ=255)
                                        tCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGt                                                                                             <  1:533335/140‑1 (MQ=255)
                                                 aCATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAATTCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCAcc                                                                                    <  1:1691397/140‑1 (MQ=255)
                                                     cTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCCTGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCaa                                                                                <  1:1034265/140‑1 (MQ=255)
                                                     cTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCaa                                                                                <  1:1034243/140‑1 (MQ=255)
                                                     cTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCaa                                                                                <  1:1034212/140‑1 (MQ=255)
                                                           aataccctTTCCCAAGGATGCGCGGTTCCAATTGTCCCATTTCGCCATGTGGCTGCAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCCCATTCCACCTGATTGGCAGCTGGTGCGCCACCACAACGAAAt                                                                          <  2:1644658/134‑1 (MQ=255)
                                                             ttCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGt                                                                        <  2:1268928/140‑1 (MQ=255)
                                                               ctctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGttt                                                                      >  1:1586181/1‑140 (MQ=255)
                                                               ctctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGttt                                                                      >  1:1086005/1‑140 (MQ=255)
                                                               ctctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGttt                                                                      >  2:1169750/1‑140 (MQ=255)
                                                                tctTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTc                                                                     >  2:642930/1‑140 (MQ=255)
                                                                  tttCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTg                                                                   >  2:138431/1‑140 (MQ=255)
                                                                    tCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGa                                                                 <  1:1616868/140‑1 (MQ=255)
                                                                             tGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTACATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTc                                                        <  2:889454/140‑1 (MQ=255)
                                                                             tGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTc                                                        <  2:540480/140‑1 (MQ=255)
                                                                             tGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTc                                                        <  2:230872/140‑1 (MQ=255)
                                                                                gcgGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATa                                                     >  1:1247798/1‑140 (MQ=255)
                                                                                      ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                               <  1:780967/140‑1 (MQ=255)
                                                                                      ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                               <  1:64584/140‑1 (MQ=255)
                                                                                      ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                               <  1:1694568/140‑1 (MQ=255)
                                                                                      ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                               <  1:1694394/140‑1 (MQ=255)
                                                                                      ccATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAAc                                               <  2:1743317/140‑1 (MQ=255)
                                                                                                   tACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGtt                                  >  2:484655/1‑140 (MQ=255)
                                                                                                     cGCCATGAGGCTGTAAATCTGACGTCGAGATCCCTTTGCCAACTTCCTGCCAGGCACCTTCCACCTGATTGTCATCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAg                                <  1:1706211/140‑1 (MQ=255)
                                                                                                         aTGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGAAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGaa                            <  1:1020720/140‑1 (MQ=255)
                                                                                                          tGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat                           >  2:851340/1‑140 (MQ=255)
                                                                                                          tGAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAat                           >  2:1399838/1‑140 (MQ=255)
                                                                                                           gAGGCTGTAAATCTGAGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatg                          <  1:1010754/140‑1 (MQ=255)
                                                                                                              gCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAatgatg                       <  2:1016762/140‑1 (MQ=255)
                                                                                                                tGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGa                     >  2:825093/1‑140 (MQ=255)
                                                                                                                tGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGa                     >  1:1320328/1‑140 (MQ=255)
                                                                                                                    aaTCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAAc                 <  2:1056737/140‑1 (MQ=255)
                                                                                                                    aaTCTGAGGTCGAGATCCCTTTGCCATCTTCCTGACAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAATAAATGTTTCTTGAAATGCTTTCATAATTAACACCTTTTATCGTTAGCGAATGATGGATAAc                 <  2:79544/140‑1 (MQ=255)
                                                                                                                          aGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt           >  2:771732/1‑140 (MQ=255)
                                                                                                                          aGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt           >  2:1299260/1‑140 (MQ=255)
                                                                                                                          aGGTCGAGATCCCTTTACCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGt           >  2:1299098/1‑140 (MQ=255)
                                                                                                                              cGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCAc       <  2:1450647/140‑1 (MQ=255)
                                                                                                                              cGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCATTTTATCGTTAGCGAATGATGGATAACAGCAGTTCAc       <  2:1753282/140‑1 (MQ=255)
                                                                                                                                agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCt     >  2:295279/1‑140 (MQ=255)
                                                                                                                                agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCt     >  2:295318/1‑140 (MQ=255)
                                                                                                                                agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCt     >  2:1621978/1‑140 (MQ=255)
                                                                                                                                agaTCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCt     >  2:1475998/1‑140 (MQ=255)
                                                                                                                                 caTCCCATTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAAAGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTg    <  2:1881209/139‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCTTCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTATCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:1568477/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGACATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:752428/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:1191583/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:281772/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:1191602/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:752553/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:1191613/140‑1 (MQ=255)
                                                                                                                                   tCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTgcg  <  2:1679133/140‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                    
GGCAAATAACGCGATATCTTCCGGGTAACGGTGATAAAAATCGATGGCGACATCTTTGATATTCTCTTTCCCCAGGATGCGCGGTTCCATTTTTCCCATTACGCCATGAGGCTGTAAATCTGAGGTCGAGATCCCTTTGCCATCTTCCTGCCAGGCACCTTCCACCTGATTGGCAGCTGTTGCGCCACCCCAAAGAAATGTTTCTGGAAATGCTTTCATAATTAACTCCTTTTATCGTTAGCGAATGATGGATAACAGCGGTTCACCTGCG  >  NC_000913/3903483‑3903753

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