Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,358,926 T→C T3A (ACC→GCC)  b1297 ← putative glutamine synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,358,9260TC99.5% 648.9 / ‑5.1 207T3A (ACC→GCC) b1297putative glutamine synthetase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (120/85);  total (122/85)
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.09e-01

Only 100 of 214 total aligned reads displayed.

GCTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGTTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTAA  >  NC_000913/1358836‑1359015
                                                                                          |                                                                                         
 cTCCGGTTAGTGTTTCACTCCCTGCAAACGGGCCCTCCCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTc                                                                                        <  2:4312374/93‑1 (MQ=255)
 cTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTc                                                                                        >  1:5579421/1‑93 (MQ=255)
      ggTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGa                                                                                   >  1:560751/1‑93 (MQ=255)
      ggTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGa                                                                                   >  1:2926104/1‑93 (MQ=255)
       gTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGAt                                                                                  >  2:1309587/1‑93 (MQ=255)
       gTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGAt                                                                                  >  1:416539/1‑93 (MQ=255)
             tttCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCg                                                                            >  2:3448100/1‑93 (MQ=255)
                 aCTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCtt                                                                        >  1:2791763/1‑93 (MQ=255)
                   ttCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTg                                                                      >  1:2062413/1‑93 (MQ=255)
                        gCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttg                                                                 >  1:3854269/1‑93 (MQ=255)
                          gTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttgtt                                                               >  1:6104335/1‑93 (MQ=255)
                          gTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttgtt                                                               >  1:2784193/1‑93 (MQ=255)
                          gTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttgtt                                                               >  2:4927893/1‑93 (MQ=255)
                           tAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttgttt                                                              >  2:7228893/1‑93 (MQ=255)
                           tAAAGGCGATACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGgtttgttt                                                              >  2:6303095/1‑93 (MQ=255)
                            aaaGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGtttt                                                             >  1:4897438/1‑93 (MQ=255)
                                  gcTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTc                                                       >  2:1618053/1‑93 (MQ=255)
                                    tACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGt                                                     >  2:3148790/1‑93 (MQ=255)
                                     cccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGtt                                                    <  2:5124515/92‑1 (MQ=255)
                                      ccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTa                                                   <  1:5838886/93‑1 (MQ=255)
                                      ccccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTa                                                   <  1:2598191/93‑1 (MQ=255)
                                       cccGCCTCTCGTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTAt                                                  <  2:2405539/93‑1 (MQ=255)
                                        ccGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATc                                                 >  1:4747551/1‑93 (MQ=255)
                                         cGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCa                                                >  1:7671122/1‑93 (MQ=255)
                                         cGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCa                                                <  2:7556218/93‑1 (MQ=255)
                                         cGCCTCTCTTCTGACTGCTGAACAAAGTCCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCa                                                <  1:3155218/93‑1 (MQ=255)
                                          gCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTTGCTTCCATGATTCTTCGGCTTTGGTTTGTTTTCCGCTCTTTATcca                                               >  2:5957853/1‑91 (MQ=255)
                                          gCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCaa                                               <  1:5716012/93‑1 (MQ=255)
                                          gCCTCTCTTCTGACTGCTGAACAAAGCTCTCTACTTCCACGATATTGGCTTCCCTGATTGTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCaa                                               >  1:697041/1‑93 (MQ=255)
                                            ctctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAg                                             <  1:7882519/93‑1 (MQ=255)
                                            ctctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAg                                             <  2:7625914/93‑1 (MQ=255)
                                            ctctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAg                                             >  2:5797089/1‑93 (MQ=255)
                                             tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGc                                            >  2:2194873/1‑93 (MQ=255)
                                             tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGc                                            >  2:3486859/1‑93 (MQ=255)
                                             tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGc                                            >  1:3714356/1‑93 (MQ=255)
                                             tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGc                                            <  2:3639581/93‑1 (MQ=255)
                                             tctctTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATGCTTCGCCTTTGGGTTGTTTTCCGCTCGTTATTACAGc                                            >  1:7806289/1‑93 (MQ=255)
                                               tcttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGt                                          >  1:6740863/1‑93 (MQ=255)
                                                cttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTa                                         >  1:3454240/1‑93 (MQ=255)
                                                cttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTa                                         <  1:2695614/93‑1 (MQ=255)
                                                 ttctGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTaa                                        >  2:2594237/1‑93 (MQ=255)
                                                    tGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAat                                     <  2:4205372/93‑1 (MQ=255)
                                                    tGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAat                                     >  2:859021/1‑93 (MQ=255)
                                                     gACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAata                                    >  2:3368802/1‑93 (MQ=255)
                                                       ctgctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAAtata                                  >  1:7523923/1‑93 (MQ=255)
                                                        tgctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATaa                                 >  2:2647316/1‑93 (MQ=255)
                                                        tgctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATaa                                 >  1:1039742/1‑93 (MQ=255)
                                                        tgctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATaa                                 >  1:7227644/1‑93 (MQ=255)
                                                        tgctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATaa                                 <  1:6191602/93‑1 (MQ=255)
                                                        tgctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATaa                                 >  1:1625706/1‑93 (MQ=255)
                                                         gctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAAt                                <  1:5084258/93‑1 (MQ=255)
                                                          ctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCGAAGCGTAAAATATAAtt                               >  2:6596077/1‑92 (MQ=255)
                                                          ctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATg                               >  1:2424216/1‑93 (MQ=255)
                                                          ctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATg                               <  2:203650/93‑1 (MQ=255)
                                                          ctgAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATg                               <  1:4139227/93‑1 (MQ=255)
                                                            gAACAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGac                             <  1:4788166/93‑1 (MQ=255)
                                                              aCAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGacca                           <  2:4790906/93‑1 (MQ=255)
                                                              aCAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGacca                           <  1:3525669/93‑1 (MQ=255)
                                                               cAAAGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGaccac                          >  1:5193146/1‑93 (MQ=255)
                                                                 aaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGAccacca                        <  1:3008597/93‑1 (MQ=255)
                                                                 aaGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGATTGTTTTCCGCTCGTTATAAAAGCGTAAAATATAATGAtcacca                        >  2:703153/1‑93 (MQ=255)
                                                                  acttCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCAt                       <  1:4042997/91‑1 (MQ=255)
                                                                  aGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCAt                       <  2:5754717/93‑1 (MQ=255)
                                                                  aGTTCTCTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCAt                       <  1:6651181/93‑1 (MQ=255)
                                                                    ttCTCTACTTCAACGATATTGGTTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTc                     >  2:6319062/1‑93 (MQ=255)
                                                                     tctcTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCg                    >  2:5301264/1‑93 (MQ=255)
                                                                     tctcTACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCg                    >  2:7804109/1‑93 (MQ=255)
                                                                       tctACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGaa                  >  2:6743056/1‑93 (MQ=255)
                                                                       tctACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGaa                  <  1:1697903/93‑1 (MQ=255)
                                                                        ctACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAAt                 <  1:3631394/93‑1 (MQ=255)
                                                                        ctACTTCAACGATATTGGCTTCCATAATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAAt                 <  2:4545248/93‑1 (MQ=255)
                                                                         tACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATc                <  1:4710699/93‑1 (MQ=255)
                                                                         tACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATc                <  2:3286898/93‑1 (MQ=255)
                                                                         tACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATc                >  2:6432808/1‑93 (MQ=255)
                                                                         tACTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATc                <  1:981775/93‑1 (MQ=255)
                                                                         tACTGCAGCGATATTTGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATc                <  1:214898/93‑1 (MQ=255)
                                                                          tcTTCAACGATATTGGCTTACATGATTCATCGCCTTTGGTTTGTTTTCCGCTCGTTATCACAGCGTAAAATATCATGACCACCATTCGAATCt               <  1:7024628/92‑1 (MQ=255)
                                                                          aCTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCt               <  2:2235749/93‑1 (MQ=255)
                                                                           ctgcAACGATATTGGCGTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTg              <  1:1271409/90‑1 (MQ=255)
                                                                           cTTCAACGATATTTGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTg              <  2:1872084/93‑1 (MQ=255)
                                                                           cTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTg              >  2:7373390/1‑93 (MQ=255)
                                                                           cTTCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTg              <  2:1532988/93‑1 (MQ=255)
                                                                            ttCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGt             >  1:308195/1‑93 (MQ=255)
                                                                             tCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGCATCTGTa            >  2:1731721/1‑93 (MQ=255)
                                                                             tCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTa            <  2:5039659/93‑1 (MQ=255)
                                                                             tCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTa            <  1:5935678/93‑1 (MQ=255)
                                                                             tCAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTa            <  2:5840109/93‑1 (MQ=255)
                                                                              cAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTAt           <  2:2813345/93‑1 (MQ=255)
                                                                              cAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTAt           <  1:5811427/93‑1 (MQ=255)
                                                                              cAACGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTAt           <  1:7062042/93‑1 (MQ=255)
                                                                               aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg          >  2:5500899/1‑93 (MQ=255)
                                                                               aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg          >  1:6030588/1‑93 (MQ=255)
                                                                               aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg          >  2:4321453/1‑93 (MQ=255)
                                                                               aaCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATg          >  2:816741/1‑93 (MQ=255)
                                                                                aCGATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGc         >  1:3247286/1‑93 (MQ=255)
                                                                                  gATATTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCaa       >  2:697442/1‑93 (MQ=255)
                                                                                    tatTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAAc     >  2:3758250/1‑93 (MQ=255)
                                                                                     atTGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACt    >  1:7645865/1‑93 (MQ=255)
                                                                                      ttGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTa   >  1:1819909/1‑93 (MQ=255)
                                                                                       tGGCTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTaa  >  1:5893693/1‑93 (MQ=255)
                                                                                          |                                                                                         
GCTCCAGGTAGCGTTTCACTTCCTGCGTAAAGGCGCTACCCCGCCTCTCTTCTGACTGCTGAACAAAGTTCTCTACTTCAACGATATTGGTTTCCATGATTCTTCGCCTTTGGTTTGTTTTCCGCTCGTTATCAAAGCGTAAAATATAATGACCACCATTCGAATCTGTATGCAAACTAA  >  NC_000913/1358836‑1359015

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑