Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,093,303 +G coding (187/705 nt) yggS → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,093,2991.G84.3% 632.6 / 112.2 217coding (183/705 nt)yggShypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (15/19);  new base G (85/98);  total (100/117)
Fisher's exact test for biased strand distribution p-value = 8.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 249 total aligned reads displayed.

CTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAGCCGCC  >  NC_000913/3093207‑3093390
                                                                                             |                                                                                           
ctgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGaa                                                                                              <  2:795544/93‑1 (MQ=255)
ctgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGaa                                                                                              <  2:4058872/93‑1 (MQ=255)
ctgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGaa                                                                                              <  1:6594912/93‑1 (MQ=255)
 tgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAg                                                                                            >  2:2491971/1‑93 (MQ=255)
 tgctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAg                                                                                            <  2:3819780/93‑1 (MQ=255)
  gctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgg                                                                                           >  2:5547261/1‑93 (MQ=255)
   ctTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAggg                                                                                          <  1:532063/93‑1 (MQ=255)
    ttGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggg                                                                                         <  1:5833369/93‑1 (MQ=255)
    ttGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggg                                                                                         <  1:6956734/93‑1 (MQ=255)
    ttGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggg                                                                                         >  1:2256291/1‑93 (MQ=255)
      gCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggggt                                                                                       <  1:3188658/93‑3 (MQ=255)
      gCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggggt                                                                                       <  1:3261788/93‑3 (MQ=255)
      gCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggggt                                                                                       >  1:4924395/1‑91 (MQ=255)
      gCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggggt                                                                                       <  1:4295555/93‑3 (MQ=255)
        agtcagtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAgggggtag                                                                                     >  1:3954150/1‑89 (MQ=255)
          tcagtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGAt                                                                                    >  1:1836110/1‑93 (MQ=255)
             gtAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATaaaa                                                                                >  1:2728353/1‑93 (MQ=255)
              tAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATaaaa                                                                                >  2:1727283/1‑93 (MQ=255)
               aaaacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAAt                                                                               >  1:1452597/1‑93 (MQ=255)
                aaacaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAAtt                                                                              <  2:3399495/93‑1 (MQ=255)
                  acaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCg                                                                            >  1:4154020/1‑93 (MQ=255)
                  acaaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGCAGATAAAATTCg                                                                            >  2:621552/1‑93 (MQ=255)
                   caaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGcc                                                                          >  1:6449236/1‑93 (MQ=255)
                    aaaacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCa                                                                         <  2:4442466/93‑1 (MQ=255)
                      aacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCAc                                                                        >  1:4240278/1‑93 (MQ=255)
                      aacCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCAc                                                                        <  1:4738409/93‑1 (MQ=255)
                       acCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACt                                                                       <  1:4860954/93‑1 (MQ=255)
                        ccTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACtt                                                                      >  1:5804609/1‑93 (MQ=255)
                        ccTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACtt                                                                      <  1:11983/93‑1 (MQ=255)
                          tGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTc                                                                    <  2:2302118/93‑1 (MQ=255)
                           gCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCa                                                                   <  2:1324691/93‑1 (MQ=255)
                           gCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCa                                                                   >  2:5604206/1‑93 (MQ=255)
                              aGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGa                                                                >  1:2858287/1‑93 (MQ=255)
                              aGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGa                                                                >  2:142275/1‑93 (MQ=255)
                                 gcCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACt                                                             <  2:5275760/93‑1 (MQ=255)
                                 gcCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTg                                                            <  1:2528207/93‑1 (MQ=255)
                                 gcCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTg                                                            >  1:1314861/1‑93 (MQ=255)
                                   cATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTgg                                                           <  1:452979/93‑1 (MQ=255)
                                   cATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTgg                                                           <  2:2890903/93‑1 (MQ=255)
                                    aTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTggg                                                          >  1:1282634/1‑93 (MQ=255)
                                    aTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTggg                                                          >  2:1781849/1‑93 (MQ=255)
                                    aTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTggg                                                          <  2:4826631/93‑1 (MQ=255)
                                    aTCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTggg                                                          >  2:1232986/1‑93 (MQ=255)
                                      cGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCg                                                        <  2:7229482/93‑1 (MQ=255)
                                       gCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGt                                                       >  1:5667299/1‑93 (MQ=255)
                                          gAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAAc                                                    >  1:4084005/1‑93 (MQ=255)
                                           aaGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACa                                                   >  2:2707193/1‑93 (MQ=255)
                                           aaGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACa                                                   <  1:5683821/93‑1 (MQ=255)
                                           aaGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACa                                                   <  2:3666149/93‑1 (MQ=255)
                                            aGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAg                                                  <  1:5405612/93‑1 (MQ=255)
                                             gCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAgg                                                 >  2:6843580/1‑93 (MQ=255)
                                               cATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGAt                                               <  2:4238892/93‑1 (MQ=255)
                                               cATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGAt                                               <  2:6753521/93‑1 (MQ=255)
                                               cATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGAtt                                              <  1:597280/93‑1 (MQ=255)
                                                 ttGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTa                                             <  1:1538354/93‑1 (MQ=255)
                                                    aTGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGaa                                          <  2:5718299/93‑1 (MQ=255)
                                                     tGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATg                                        <  1:1995933/93‑1 (MQ=255)
                                                       ccGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATgg                                       >  2:5467421/1‑93 (MQ=255)
                                                        cGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGc                                      <  2:3329556/93‑1 (MQ=255)
                                                         gggCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCa                                     >  1:5570732/1‑93 (MQ=255)
                                                         gggCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCa                                     >  2:3094270/1‑93 (MQ=255)
                                                          ggCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCAt                                    >  2:4814319/1‑93 (MQ=255)
                                                           gCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCAtt                                   <  1:2555944/93‑1 (MQ=255)
                                                             aGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCAtttt                                 <  1:1563692/93‑1 (MQ=255)
                                                              gCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTa                                >  2:3226430/1‑93 (MQ=255)
                                                               cGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTAtt                              >  1:4485735/1‑93 (MQ=255)
                                                                gTCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTAtt                              >  2:7202543/1‑93 (MQ=255)
                                                                 tCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTg                             >  2:6061956/1‑93 (MQ=255)
                                                                 tCAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTg                             >  2:7515942/1‑93 (MQ=255)
                                                                  cAATTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTgg                            <  1:1623770/93‑1 (MQ=255)
                                                                  cAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGc                           >  1:4641801/1‑93 (MQ=255)
                                                                   aaTTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGc                           >  2:194650/1‑93 (MQ=255)
                                                                   aaTTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGc                           >  1:3407920/1‑93 (MQ=255)
                                                                   aaTTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGc                           <  2:4198428/93‑1 (MQ=255)
                                                                    aTTTGGTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGcc                          <  2:2574465/93‑1 (MQ=255)
                                                                         gTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTg                     <  1:5176222/93‑1 (MQ=255)
                                                                         gTGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTg                     >  2:4135817/1‑93 (MQ=255)
                                                                         gTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGc                    >  2:1573562/1‑93 (MQ=255)
                                                                          tGAAAACTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGc                    >  2:2771760/1‑93 (MQ=255)
                                                                            aaaaCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAg                  <  1:3667797/93‑1 (MQ=255)
                                                                             aaaCTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTc                <  1:2919028/93‑1 (MQ=255)
                                                                              aaCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTc                <  2:87198/93‑1 (MQ=255)
                                                                              aaCTACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTc                >  1:5101482/1‑93 (MQ=255)
                                                                               aCTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCta              >  1:2207561/1‑93 (MQ=255)
                                                                                 tACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaa             >  2:2249748/1‑93 (MQ=255)
                                                                                 tACGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaa             <  1:746317/93‑1 (MQ=255)
                                                                                  aCGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGTATTTTATTGGCCCGTTGCAGTCtaat            >  1:6048424/1‑93 (MQ=255)
                                                                                   cGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaata           <  2:2911306/93‑1 (MQ=255)
                                                                                   cGTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaata           <  2:7391156/93‑1 (MQ=255)
                                                                                    gTTCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCtaataa          >  2:1019007/1‑93 (MQ=255)
                                                                                     ttCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTGATaaa         >  2:4973228/1‑93 (MQ=255)
                                                                                     ttCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATaaa         <  1:7163966/93‑1 (MQ=255)
                                                                                     ttCAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATaaa         >  2:1190277/1‑93 (MQ=255)
                                                                                     ttCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATaaaa        >  2:4629161/1‑93 (MQ=255)
                                                                                      tCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAg       >  1:2704563/1‑93 (MQ=255)
                                                                                       cAGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAg       <  1:6685942/93‑1 (MQ=255)
                                                                                       cAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAgc      <  1:3480149/93‑1 (MQ=255)
                                                                                        aGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAgc      >  1:6349871/1‑93 (MQ=255)
                                                                                        aGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAgc      <  1:5524900/93‑1 (MQ=255)
                                                                                        aGGAAGGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAgc      <  2:2020272/93‑1 (MQ=255)
                                                                                             |                                                                                           
CTGCTTGCAGTCAGTAAAACAAAACCTGCGAGCGCCATCGCAGAAGCCATTGATGCCGGGCAGCGTCAATTTGGTGAAAACTACGTTCAGGAAGGGGTAGATAAAATTCGCCACTTTCAGGAACTGGGCGTAACAGGATTAGAATGGCATTTTATTGGCCCGTTGCAGTCTAATAAAAGCCGCC  >  NC_000913/3093207‑3093390

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