Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,885,719 T→C D206G (GAT→GGT)  cysH ← 3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,885,7190TC88.0% 535.3 / 52.9 192D206G (GAT→GGT) cysH3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base T (11/12);  new base C (94/75);  total (105/87)
Fisher's exact test for biased strand distribution p-value = 5.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 199 total aligned reads displayed.

TTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATA  >  NC_000913/2885629‑2885809
                                                                                          |                                                                                          
tttAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCAccc                                                                                          <  1:502054/93‑1 (MQ=255)
 ttAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCAcccc                                                                                         <  1:2332093/93‑1 (MQ=255)
  tAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATTCTTCACCCCa                                                                                        >  1:5216080/1‑93 (MQ=255)
   aaGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCAt                                                                                       <  2:5371787/93‑1 (MQ=255)
   aaGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCAt                                                                                       <  2:2421277/93‑1 (MQ=255)
     accAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATaa                                                                                     <  1:5448360/92‑1 (MQ=255)
      ccAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAAt                                                                                    >  1:5197457/1‑93 (MQ=255)
      ccAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAAt                                                                                    >  2:5688973/1‑93 (MQ=255)
       cAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATg                                                                                   >  1:6017739/1‑93 (MQ=255)
       cAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATg                                                                                   >  2:3735223/1‑93 (MQ=255)
        aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATgg                                                                                  <  1:5575146/93‑1 (MQ=255)
        aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATgg                                                                                  >  1:5747601/1‑93 (MQ=255)
        aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATgg                                                                                  >  2:501750/1‑93 (MQ=255)
        aaaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATgg                                                                                  >  1:1225966/1‑93 (MQ=255)
         aaGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATggg                                                                                 >  2:4062549/1‑93 (MQ=255)
           gAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTg                                                                               >  2:662201/1‑93 (MQ=255)
            aaaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGa                                                                              <  2:2294961/93‑1 (MQ=255)
            aaaCGCGTTTCTTCTTCCACCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGa                                                                              <  2:270222/93‑1 (MQ=255)
             aaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGat                                                                             >  2:6564337/1‑93 (MQ=255)
             aaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGat                                                                             >  1:1098766/1‑93 (MQ=255)
             aaCGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGat                                                                             <  1:4121939/93‑1 (MQ=255)
                gcgTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATAtt                                                                          >  1:168918/1‑93 (MQ=255)
                   tttCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCa                                                                       >  1:6336238/1‑93 (MQ=255)
                   tttCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCa                                                                       >  2:2925161/1‑93 (MQ=255)
                   tttCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCa                                                                       <  1:3823359/93‑1 (MQ=255)
                   tttCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGCGATATTTCa                                                                       >  2:642843/1‑93 (MQ=255)
                    ttcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAg                                                                      >  1:4027150/1‑93 (MQ=255)
                     tcttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAgg                                                                     >  1:2091822/1‑93 (MQ=255)
                      cttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGc                                                                    <  2:4110704/93‑1 (MQ=255)
                      cttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGc                                                                    <  2:790134/93‑1 (MQ=255)
                       ttcttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGcc                                                                   >  2:5388623/1‑93 (MQ=255)
                         cttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCAt                                                                 <  1:5719334/93‑1 (MQ=255)
                          ttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATg                                                                <  1:2552193/93‑1 (MQ=255)
                          ttcCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATg                                                                >  2:1734275/1‑93 (MQ=255)
                            ccGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtt                                                              <  1:1118179/93‑1 (MQ=255)
                              gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGtttt                                                            <  1:4744679/93‑1 (MQ=255)
                              gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGtttt                                                            <  1:2368018/93‑1 (MQ=255)
                              gCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGtttt                                                            <  2:6518522/93‑1 (MQ=255)
                                   gCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAg                                                       >  2:319875/1‑93 (MQ=255)
                                    ccGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAgg                                                      >  2:5905560/1‑93 (MQ=255)
                                    ccGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAgg                                                      >  1:6120840/1‑93 (MQ=255)
                                      gggTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTa                                                    >  2:3291318/1‑93 (MQ=255)
                                      gggTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTa                                                    >  2:1239955/1‑93 (MQ=255)
                                      gggTTCCCAATTACGGGTTGTATGGGTATCGTCCACCGATAAAAATCCGTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTa                                                    >  2:6655282/1‑93 (MQ=255)
                                       ggTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTAc                                                   <  2:3066488/93‑1 (MQ=255)
                                       ggTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTAc                                                   <  1:3137354/93‑1 (MQ=255)
                                        gTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACt                                                  <  2:2039904/93‑1 (MQ=255)
                                          tCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGa                                                >  2:6333731/1‑93 (MQ=255)
                                           cccATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGAt                                               >  1:5900720/1‑93 (MQ=255)
                                           cccATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGAt                                               >  2:4197717/1‑93 (MQ=255)
                                            ccATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATTTTTTTGCAGGTACTGATa                                              >  2:5823832/1‑93 (MQ=255)
                                             cATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATaa                                             <  2:3981028/93‑1 (MQ=255)
                                              aTTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATaaa                                            <  2:6299617/93‑1 (MQ=255)
                                                ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATa                                          <  2:1935561/93‑1 (MQ=255)
                                                ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATa                                          <  2:3583391/93‑1 (MQ=255)
                                                ttACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATa                                          <  2:3761082/93‑1 (MQ=255)
                                                 tACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAg                                         <  1:3055469/93‑1 (MQ=255)
                                                  ccGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGt                                        <  1:1361036/92‑1 (MQ=255)
                                                   cGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGtt                                       <  2:4818420/93‑1 (MQ=255)
                                                   cGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGtt                                       <  1:3357774/93‑1 (MQ=255)
                                                     ggTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCg                                     >  1:785283/1‑93 (MQ=255)
                                                      gTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCgg                                    <  1:5263235/93‑1 (MQ=255)
                                                      gTTGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCgg                                    <  1:2249935/93‑1 (MQ=255)
                                                        tGTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGtt                                  >  1:6442934/1‑93 (MQ=255)
                                                         gTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTa                                 <  1:4326475/93‑1 (MQ=255)
                                                         gTATGGGTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTa                                 >  2:4894063/1‑93 (MQ=255)
                                                             gggTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATccc                             >  2:2584449/1‑93 (MQ=255)
                                                              ggTATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCa                            <  1:436053/93‑1 (MQ=255)
                                                               gTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAg                           <  2:4912746/93‑1 (MQ=255)
                                                                tATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGt                          <  1:1574898/93‑1 (MQ=255)
                                                                tATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGt                          <  1:885127/93‑1 (MQ=255)
                                                                tATCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGt                          >  2:2645324/1‑93 (MQ=255)
                                                                 aTCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTc                         <  2:6011309/93‑1 (MQ=255)
                                                                 aTCGCCCACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTc                         >  2:1112313/1‑93 (MQ=255)
                                                                     cccACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                     >  2:4632620/1‑93 (MQ=255)
                                                                     cccACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                     >  2:2276543/1‑93 (MQ=255)
                                                                     cccACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATa                     <  1:6530379/93‑1 (MQ=255)
                                                                      ccACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATaa                    >  1:1874261/1‑93 (MQ=255)
                                                                       cACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAAt                   >  2:5581823/1‑93 (MQ=255)
                                                                       cACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAAt                   >  1:3648383/1‑93 (MQ=255)
                                                                       cACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAAt                   <  1:2502171/93‑1 (MQ=255)
                                                                       cACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAAt                   >  2:3401759/1‑93 (MQ=255)
                                                                       cACCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTAGGTTATCCCAGTCGATAAt                   >  2:5654469/1‑93 (MQ=255)
                                                                        aCCGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATc                  <  1:3945525/93‑1 (MQ=255)
                                                                          cGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCgg                >  1:5787663/1‑93 (MQ=255)
                                                                          cGATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAACAGTTCGGTTATCCCAGTCGATAATCgg                <  2:1044550/93‑1 (MQ=255)
                                                                           gATAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgc               >  2:1394458/1‑93 (MQ=255)
                                                                            aTAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgca              >  2:3907977/1‑93 (MQ=255)
                                                                            aTAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgca              >  2:2699080/1‑93 (MQ=255)
                                                                             tAAATATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcag             >  1:3489236/1‑93 (MQ=255)
                                                                              aaaTATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagc            <  1:3691891/93‑1 (MQ=255)
                                                                               aaTATCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGgcagca           >  2:922388/1‑93 (MQ=255)
                                                                                atatCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCAc          <  2:3676293/93‑1 (MQ=255)
                                                                                acatCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCAc          <  1:6379438/91‑1 (MQ=255)
                                                                                 tatCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACt         >  1:3066272/1‑93 (MQ=255)
                                                                                 tatCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACt         <  2:438562/93‑1 (MQ=255)
                                                                                  atCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACtt        <  1:4497076/93‑1 (MQ=255)
                                                                                   tCCTTCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACttt       >  2:131708/1‑93 (MQ=255)
                                                                                      ttCACCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTaaa    >  1:5885176/1‑93 (MQ=255)
                                                                                       tCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAAt   >  1:427276/1‑93 (MQ=255)
                                                                                          |                                                                                          
TTTAAGCCAAAGAAACGCGTTTCTTCTTCCGCCATGCCGGGTTCCCATTTACGGGTTGTATGGGTATCGCCCACCGATAAATATCCTTCATCCCATAATGGGTGATATTTCAGGCCATGTTTTTGCAGGTACTGATAAATAGTTCGGTTATCCCAGTCGATAATCGGCAGCACTTTAAATA  >  NC_000913/2885629‑2885809

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