Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,858,802 T→C P185P (CCA→CCG b1776 ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,858,8020TC84.7% 583.0 / 90.7 255P185P (CCA→CCGb1776putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (18/21);  new base C (99/117);  total (117/138)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 256 total aligned reads displayed.

GGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCAAGTAC  >  NC_000913/1858711‑1858894
                                                                                           |                                                                                            
 gAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCgg                                                                                            >  2:1123442/1‑93 (MQ=255)
 gAAGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCgg                                                                                            >  2:4887900/1‑93 (MQ=255)
   aGCGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGc                                                                                          <  2:2899051/93‑1 (MQ=255)
     cGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCa                                                                                        <  1:1361453/93‑1 (MQ=255)
     cGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCa                                                                                        <  2:3710287/93‑1 (MQ=255)
     cGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCa                                                                                        <  2:4685318/93‑1 (MQ=255)
     cGGACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATGCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCa                                                                                        <  1:5432965/93‑1 (MQ=255)
      ggACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAg                                                                                       <  2:1831805/93‑1 (MQ=255)
      ggACCGTCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAg                                                                                       <  2:2737879/93‑1 (MQ=255)
          cgtcacaTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGcc                                                                                   <  1:4299569/90‑1 (MQ=255)
            tcacaTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGa                                                                                 >  2:5908288/2‑93 (MQ=255)
            tcacaTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGa                                                                                 <  2:1430984/92‑1 (MQ=255)
              acaTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATg                                                                               >  1:1893141/1‑93 (MQ=255)
               caTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATga                                                                              <  1:1788579/93‑1 (MQ=255)
               caTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATga                                                                              >  1:7105043/1‑93 (MQ=255)
               caTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATga                                                                              <  1:1820116/93‑1 (MQ=255)
               caTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATga                                                                              >  2:1645119/1‑93 (MQ=255)
               caTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATga                                                                              >  2:6084921/1‑93 (MQ=255)
                aTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATgac                                                                             <  2:6982601/93‑1 (MQ=255)
                aTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATgac                                                                             >  1:1067228/1‑93 (MQ=255)
                aTCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCAATTTGTATGGAGAACAGCCCGAGCGGGCCAGTGCCGATgac                                                                             <  2:5162669/93‑1 (MQ=255)
                 tcttctTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATgacg                                                                            >  2:2430751/1‑93 (MQ=255)
                  cttcttGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATgacga                                                                           <  2:3320316/93‑1 (MQ=255)
                   ttcttGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATgacgac                                                                          <  2:2629021/93‑1 (MQ=255)
                   ttcttGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATgacgac                                                                          >  2:1003458/1‑93 (MQ=255)
                      ttGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACcac                                                                       >  1:1516566/1‑93 (MQ=255)
                      ttGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTCTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACcac                                                                       <  1:5973748/93‑1 (MQ=255)
                       tGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACcaca                                                                      <  2:1256429/93‑1 (MQ=255)
                       tGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACcaca                                                                      <  1:5757416/93‑1 (MQ=255)
                        gCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACAt                                                                     <  1:4088461/93‑1 (MQ=255)
                          aGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATcc                                                                   >  2:4948407/1‑93 (MQ=255)
                          aGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATcc                                                                   <  1:707834/93‑1 (MQ=255)
                             ccAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGa                                                                >  2:5865941/1‑93 (MQ=255)
                              cAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGAc                                                               >  1:1918370/1‑93 (MQ=255)
                              cAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGAc                                                               <  1:6241966/93‑1 (MQ=255)
                               aaCGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGAcc                                                              <  2:2521446/93‑1 (MQ=255)
                                acgacgacgATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCGGTGCCGACGACGACCACATCCTGACCa                                                             <  2:6416788/93‑1 (MQ=255)
                                acgacgacgATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCa                                                             <  1:3309427/93‑1 (MQ=255)
                                acgacgacgATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCa                                                             <  2:777060/93‑1 (MQ=255)
                                acgacgacgATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCa                                                             >  2:1694319/1‑93 (MQ=255)
                                 cgacgacgaTATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCag                                                            <  1:548982/93‑1 (MQ=255)
                                   acgacgaTATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCagga                                                          >  1:6549226/1‑93 (MQ=255)
                                   acgacgaTATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCagga                                                          <  2:2366037/93‑1 (MQ=255)
                                    cgacgaTATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCaggaa                                                         <  2:3440450/93‑1 (MQ=255)
                                      acgaTATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCaggaagg                                                       <  1:4626928/93‑1 (MQ=255)
                                      acgaTATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCaggaagg                                                       <  2:5426358/93‑1 (MQ=255)
                                         atatTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAAt                                                    >  2:7326111/1‑93 (MQ=255)
                                          tatTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAAtt                                                   >  2:7262551/1‑93 (MQ=255)
                                            tttACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTc                                                 <  1:2511235/93‑1 (MQ=255)
                                            tttACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTc                                                 <  1:1011906/93‑1 (MQ=255)
                                              tACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGCGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGa                                               <  1:7301876/93‑1 (MQ=255)
                                              tACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGa                                               <  1:3901286/93‑1 (MQ=255)
                                                ccGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGAtt                                             <  1:3117771/93‑1 (MQ=255)
                                                ccGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGAtt                                             <  1:494769/93‑1 (MQ=255)
                                                    ccccATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATtgctg                                         >  1:7397938/1‑93 (MQ=255)
                                                    ccccATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATtgctg                                         <  1:6596805/93‑1 (MQ=255)
                                                    ccccATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATtgctg                                         <  2:4384015/93‑1 (MQ=255)
                                                     cccATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTACGATtgctgc                                        <  2:6310545/93‑1 (MQ=255)
                                                     cccATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATtgctgc                                        >  2:6879888/1‑93 (MQ=255)
                                                       cATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTgcgc                                      >  2:6528066/1‑93 (MQ=255)
                                                         tAATTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGa                                    >  2:5924342/1‑93 (MQ=255)
                                                          aaTTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGAt                                   >  1:5556282/1‑93 (MQ=255)
                                                          aaTTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGAt                                   >  2:2348085/1‑93 (MQ=255)
                                                           aTTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATg                                  >  1:784509/1‑93 (MQ=255)
                                                           aTTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATg                                  >  1:5169643/1‑93 (MQ=255)
                                                           aTTCGCGCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATg                                  <  2:5197414/93‑1 (MQ=255)
                                                              cgcgCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGAt                               <  2:6746199/93‑1 (MQ=255)
                                                              cgcgCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGAt                               <  2:5590775/93‑1 (MQ=255)
                                                              cgcgCCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGAt                               <  1:4177592/93‑1 (MQ=255)
                                                               gcgcCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGAtt                              <  2:4274183/93‑1 (MQ=255)
                                                               gcgcCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGAtt                              <  2:2035417/93‑1 (MQ=255)
                                                                cgcCATTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGAttt                             >  1:4370924/1‑93 (MQ=255)
                                                                    aTTTGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAg                         >  1:1216809/1‑93 (MQ=255)
                                                                     tttGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAgg                        >  1:6906355/1‑93 (MQ=255)
                                                                     tttGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAgg                        >  1:7494651/1‑93 (MQ=255)
                                                                       tGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCa                      <  1:5312613/93‑1 (MQ=255)
                                                                       tGTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCa                      >  1:639301/1‑93 (MQ=255)
                                                                        gTACGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCAt                     >  1:2970250/1‑93 (MQ=255)
                                                                           cGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTAc                  <  2:1531364/93‑1 (MQ=255)
                                                                           cGGAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTAc                  >  2:1164814/1‑93 (MQ=255)
                                                                            ggAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACa                 >  2:3894754/1‑93 (MQ=255)
                                                                            ggAGAACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACa                 >  2:7405956/1‑93 (MQ=255)
                                                                              agaACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGa               >  1:2718342/1‑93 (MQ=255)
                                                                              agaACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGa               <  1:4994943/93‑1 (MQ=255)
                                                                               gaACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGAt              >  2:736159/1‑93 (MQ=255)
                                                                               gaACAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGAt              >  2:5744317/1‑93 (MQ=255)
                                                                                aaCAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATa             <  1:5090591/93‑1 (MQ=255)
                                                                                aaCAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATa             <  1:4767398/93‑1 (MQ=255)
                                                                                aaCAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATa             <  2:223832/93‑1 (MQ=255)
                                                                                aaCAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATa             >  2:3749364/1‑93 (MQ=255)
                                                                                  cAGCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAgg           <  1:4033850/93‑1 (MQ=255)
                                                                                    gCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGt         <  1:6705878/93‑1 (MQ=255)
                                                                                    gCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGt         >  1:1457244/1‑93 (MQ=255)
                                                                                    gCCCGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGt         <  2:5448086/93‑1 (MQ=255)
                                                                                      ccGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCa       >  1:413206/1‑93 (MQ=255)
                                                                                       cGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCaa      <  1:3398255/93‑1 (MQ=255)
                                                                                       cGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCaa      <  2:5009037/93‑1 (MQ=255)
                                                                                       cGAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCaa      <  2:206860/93‑1 (MQ=255)
                                                                                        gAGCGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCAAg     >  1:4377673/1‑93 (MQ=255)
                                                                                           tGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCAAGTAc  <  1:197737/93‑1 (MQ=255)
                                                                                           |                                                                                            
GGAAGCGGACCGCCACATCTTCTTGCAGACCAACGACGACGATATTTACCGCCCCCATAATTCGCGCCATTTGTACGGAGAACAGCCCGAGTGGGCCAGTGCCGATGACGACCACATCCTGACCAGGAAGGAATTTCGATTGCTGCGCGATGGATTTGTAGGCATTACAGATAGGGTCAAGTAC  >  NC_000913/1858711‑1858894

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