Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,171,379 C→T A240V (GCC→GTC)  birA → biotin‑[acetyl‑CoA carboxylase] holoenzyme synthetase and biotin operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,171,3790CT95.5% 1598.6 / 25.3 495A240V (GCC→GTC) birAbiotin‑[acetyl‑CoA carboxylase] holoenzyme synthetase and biotin operon repressor
Reads supporting (aligned to +/- strand):  ref base C (15/7);  new base T (275/197);  total (290/205)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-01

Only 100 of 522 total aligned reads displayed.

CAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAAT  >  NC_000913/4171289‑4171470
                                                                                          |                                                                                           
cAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCa                                                                                           <  1:1274141/93‑1 (MQ=255)
cAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCa                                                                                           <  2:2675554/93‑1 (MQ=255)
cAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCa                                                                                           <  1:3718563/93‑1 (MQ=255)
 aaTGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCAt                                                                                          <  1:1705577/93‑1 (MQ=255)
 aaTGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCAt                                                                                          >  1:2177972/1‑93 (MQ=255)
  aTGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATg                                                                                         >  2:5300035/1‑93 (MQ=255)
  aTGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATAAATCTCGAGCGTAATACGTTGGCGGTCATg                                                                                         >  1:3189288/1‑93 (MQ=255)
   tGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGc                                                                                        >  1:7083455/1‑93 (MQ=255)
   tGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGc                                                                                        <  1:4577171/93‑1 (MQ=255)
   tGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGc                                                                                        >  1:1108904/1‑93 (MQ=255)
       ccGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAAt                                                                                    <  1:4245391/93‑1 (MQ=255)
        cGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATa                                                                                   <  1:3360770/93‑1 (MQ=255)
         gtgtTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATAc                                                                                  <  2:1469227/93‑1 (MQ=255)
         gtgtTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATAc                                                                                  >  1:918634/1‑93 (MQ=255)
          tgtTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGGCATGCTAATACg                                                                                 >  2:6307692/1‑93 (MQ=255)
            ttGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTg                                                                               >  1:6362098/1‑93 (MQ=255)
              gAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGaa                                                                             >  1:5032560/1‑93 (MQ=255)
               aaGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAAt                                                                            >  2:5673805/1‑93 (MQ=255)
               aaGAGAGTGTCGTTAATCAGGGGGGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAAt                                                                            >  2:10045354/1‑93 (MQ=255)
                   gagTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACg                                                                        >  2:5470092/1‑93 (MQ=255)
                   gagTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACg                                                                        <  1:4733382/93‑1 (MQ=255)
                   gagTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACg                                                                        <  1:4790006/93‑1 (MQ=255)
                     gtgtCGTTAATCAGGGGTGGATCACGCTGCAGGAGGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATTCTAATACGTGAATTACGtg                                                                      <  1:196085/93‑1 (MQ=255)
                     gtgtCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGtg                                                                      <  2:3725604/93‑1 (MQ=255)
                       gtCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGtgct                                                                    >  1:6679825/1‑93 (MQ=255)
                        tCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGtgctg                                                                   <  1:9752758/93‑1 (MQ=255)
                           ttAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGt                                                                <  2:1210564/93‑1 (MQ=255)
                           ttAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTAAATTACGTGCTGCGt                                                                <  2:5065386/93‑1 (MQ=255)
                            tAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGtt                                                               >  2:957080/1‑93 (MQ=255)
                              aTCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTgg                                                             <  2:2469658/93‑1 (MQ=255)
                               tCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCTTGCTAATACGTGAATTAAGTGCTGCGGTGGa                                                            >  1:7421717/1‑93 (MQ=255)
                               tCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGa                                                            >  2:4243663/1‑93 (MQ=255)
                               tCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGa                                                            >  2:7252430/1‑93 (MQ=255)
                               tCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGa                                                            >  1:7160169/1‑93 (MQ=255)
                               tCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACCTGAATTACGTGCTGCGTTGGa                                                            >  1:8360121/1‑93 (MQ=255)
                                  ggggTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAAct                                                         >  1:1578272/1‑93 (MQ=255)
                                  ggggTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAAct                                                         >  1:2582652/1‑93 (MQ=255)
                                  ggggTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAAct                                                         >  1:7861656/1‑93 (MQ=255)
                                   gggTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAActc                                                        <  1:4104756/93‑1 (MQ=255)
                                    ggTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAActct                                                       >  2:9367845/1‑93 (MQ=255)
                                    ggTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAActct                                                       >  1:9371004/1‑93 (MQ=255)
                                      tGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTc                                                     <  1:8629637/93‑1 (MQ=255)
                                         aTCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGaa                                                  >  1:7288274/1‑93 (MQ=255)
                                          tCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAAc                                                 >  1:3520752/1‑93 (MQ=255)
                                           cACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACa                                                <  1:7745756/93‑1 (MQ=255)
                                           cACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACa                                                <  1:5404663/93‑1 (MQ=255)
                                            aCGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACaa                                               >  2:8371935/1‑93 (MQ=255)
                                              gCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACaaga                                             <  2:2920099/93‑1 (MQ=255)
                                              gCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACaaga                                             <  2:8148414/93‑1 (MQ=255)
                                                tGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACaagaag                                           <  1:4377518/93‑1 (MQ=255)
                                                 gCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGCACTCTTCGAACAAGAAgg                                          >  1:5644972/1‑93 (MQ=255)
                                                 gCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAgg                                          >  1:2361779/1‑93 (MQ=255)
                                                 gCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAgg                                          >  1:5761569/1‑93 (MQ=255)
                                                   aGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGAt                                        <  2:3416722/93‑1 (MQ=255)
                                                    ggAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGCAGGAtt                                       >  2:9746209/1‑93 (MQ=255)
                                                    ggAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGAtt                                       <  2:6773979/93‑1 (MQ=255)
                                                    ggAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGAtt                                       >  2:4024526/1‑93 (MQ=255)
                                                    ggAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGAtt                                       >  2:5939136/1‑93 (MQ=255)
                                                     gAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTg                                      >  2:1666313/1‑93 (MQ=255)
                                                     gAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTg                                      >  2:4755774/1‑93 (MQ=255)
                                                      aaGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTgg                                     <  1:6491775/93‑1 (MQ=255)
                                                       aGCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTTGAACTCTTCGAACAAGAAGGATTgtc                                    >  2:9018942/1‑91 (MQ=255)
                                                        gCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTGGAACAAGAAGGATTGGCa                                   <  2:5728075/93‑1 (MQ=255)
                                                        gCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCa                                   >  2:86386/1‑93 (MQ=255)
                                                        gCGGGGATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCCTCGAACCAGAAGGAttatca                                   >  1:6208111/1‑89 (MQ=255)
                                                           gggATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCt                                <  2:2291594/93‑1 (MQ=255)
                                                           gggATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCt                                >  2:2861707/1‑93 (MQ=255)
                                                           gggATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCt                                <  2:4330283/93‑1 (MQ=255)
                                                             gATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTa                              <  1:5955439/93‑1 (MQ=255)
                                                             gATCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTa                              >  1:8178615/1‑93 (MQ=255)
                                                              aTCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTAt                             >  2:641126/1‑93 (MQ=255)
                                                              aTCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTAt                             >  2:6888177/1‑93 (MQ=255)
                                                              aTCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTAt                             >  1:8031678/1‑93 (MQ=255)
                                                              aTCAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTAt                             >  2:4627681/1‑93 (MQ=255)
                                                                cAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACCAGAAGGATTGGCACCTTATCt                           >  1:6008281/1‑93 (MQ=255)
                                                                cAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCt                           <  2:6125216/93‑1 (MQ=255)
                                                                cAATCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCt                           >  1:3136046/1‑93 (MQ=255)
                                                                 aaTCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGACACCTTATCTg                          >  2:3183678/1‑93 (MQ=255)
                                                                  aTCTCGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGt                         <  2:7528693/93‑1 (MQ=255)
                                                                   tctcGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTc                        >  1:9521254/1‑93 (MQ=255)
                                                                   tctcGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTc                        <  2:384337/93‑1 (MQ=255)
                                                                   tctcGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTc                        >  1:2851959/1‑93 (MQ=255)
                                                                   tctcGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTc                        >  2:4973912/1‑93 (MQ=255)
                                                                      cGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTcgcg                     <  2:3773003/93‑1 (MQ=255)
                                                                      cGATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACATTATCTGTcgcg                     >  2:7600201/1‑93 (MQ=255)
                                                                       gATCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCgcgc                    >  2:8736842/1‑93 (MQ=255)
                                                                       gATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCgcgc                    >  2:2027074/1‑93 (MQ=255)
                                                                         tCGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTg                  >  1:6783987/1‑93 (MQ=255)
                                                                         tCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGATCTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTg                  <  1:693421/93‑1 (MQ=255)
                                                                          cGTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTgg                 >  1:5141113/1‑93 (MQ=255)
                                                                           gTAATACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTggg                >  2:3018039/1‑93 (MQ=255)
                                                                              aTACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaa             >  1:8899864/1‑93 (MQ=255)
                                                                              aTACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaa             >  2:9236581/1‑93 (MQ=255)
                                                                               tACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaaa            <  1:8953650/93‑1 (MQ=255)
                                                                               tACGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGaaaa            >  2:702669/1‑93 (MQ=255)
                                                                                 cGTTGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCCCCTTATCTGTCGCGCTGGGAAAAGc          <  1:2079797/93‑1 (MQ=255)
                                                                                    tGGCGGTCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTgg       <  2:1302655/93‑1 (MQ=255)
                                                                                     ggcggtCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGa      >  1:6995270/1‑93 (MQ=255)
                                                                                      gcggtCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGAt     >  2:7360321/1‑93 (MQ=255)
                                                                                       cggtCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATa    >  2:4419395/1‑93 (MQ=255)
                                                                                          |                                                                                           
CAATGCGCCGTGTTGAAGAGAGTGTCGTTAATCAGGGGTGGATCACGCTGCAGGAAGCGGGGATCAATCTCGATCGTAATACGTTGGCGGCCATGCTAATACGTGAATTACGTGCTGCGTTGGAACTCTTCGAACAAGAAGGATTGGCACCTTATCTGTCGCGCTGGGAAAAGCTGGATAAT  >  NC_000913/4171289‑4171470

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