Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,472,190 G→C 28.9% A69G (GCC→GGC)  pyrB ← aspartate carbamoyltransferase, catalytic subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,1900GC28.9% 124.8 / 73.9 83A69G (GCC→GGC) pyrBaspartate carbamoyltransferase, catalytic subunit
Reads supporting (aligned to +/- strand):  ref base G (29/30);  new base C (13/11);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

GCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGT  >  NC_000913/4472091‑4472288
                                                                                                   |                                                                                                  
gCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTgg                                                                                                    <  2:2244303/100‑1 (MQ=255)
 cATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGc                                                                                                   >  2:2050308/1‑100 (MQ=255)
     gACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGccccc                                                                                               >  1:733112/1‑100 (MQ=255)
       cGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAg                                                                                             >  1:1917994/1‑100 (MQ=255)
        gTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCACAACACGCTGGCCCCCAg                                                                                             >  2:1797340/1‑99 (MQ=255)
         tAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGc                                                                                           <  2:1143934/100‑1 (MQ=255)
          aaGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCg                                                                                          <  2:1881832/100‑1 (MQ=255)
             tGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTg                                                                                       <  1:1953575/100‑1 (MQ=255)
             tGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTg                                                                                       >  2:2080158/1‑100 (MQ=255)
              gCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGc                                                                                      >  1:1810082/1‑100 (MQ=255)
              gCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGc                                                                                      <  2:2138005/100‑1 (MQ=255)
              gCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGc                                                                                      <  2:1504161/100‑1 (MQ=255)
               cTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCa                                                                                     <  2:1833574/100‑1 (MQ=255)
                tGATAACCGACATGGGATCGGCCAGCGTTTCCCCCTTTTTACCCAGTGATTTATTGGCGCTGTCGTAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCAt                                                                                    <  2:1806121/100‑1 (MQ=255)
                tGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCAt                                                                                    <  1:1360900/100‑1 (MQ=255)
                tGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCAt                                                                                    <  1:1577006/100‑1 (MQ=255)
                tGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCAt                                                                                    <  1:1057085/100‑1 (MQ=255)
                tGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCAt                                                                                    <  1:367924/100‑1 (MQ=255)
                 gATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATa                                                                                   >  1:835248/1‑100 (MQ=255)
                 gATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATa                                                                                   <  1:484522/100‑1 (MQ=255)
                 gATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATATATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATa                                                                                   <  2:514659/100‑1 (MQ=255)
                      ncGAAATGGTATCGGCCAGCGTTTGGCCCTTTTTACCAAGAGATGTATTAGCGCTGTCGGAGAAGCACACAACGCTGGCCCCCAGGCGGTGCATAGATGt                                                                              >  2:539898/2‑100 (MQ=255)
                      ccGAAATGGTATCGGCCAGCGTTTCGCGCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGAt                                                                                <  2:1536628/98‑1 (MQ=255)
                      ccGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGt                                                                              >  1:651959/1‑100 (MQ=255)
                        gAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGttt                                                                            >  2:1681227/1‑100 (MQ=255)
                        gAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTAGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGt                                                                              >  2:300179/1‑98 (MQ=255)
                         aaaTGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTc                                                                           <  2:91862/100‑1 (MQ=255)
                         aaaTGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTc                                                                           <  1:1086369/100‑1 (MQ=255)
                           aTGGTATCGGCCGGCGTTTGGCCCTTTTTCCCCAGTGGTGTATTGGCGCTGTGGGAGAATTCCACCACTCTGGCCTCCATGCGGTTCATAGAGGTTTCGa                                                                         <  2:301885/100‑1 (MQ=255)
                           aTGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGa                                                                         >  1:1677988/1‑100 (MQ=255)
                              gTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAag                                                                      >  1:1557482/1‑100 (MQ=255)
                              gTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAag                                                                      >  1:1523780/1‑100 (MQ=255)
                              gTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAag                                                                      <  1:1499795/100‑1 (MQ=255)
                                aTCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAagag                                                                    >  1:381925/1‑100 (MQ=255)
                                   ggCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAgg                                                                 >  2:963460/1‑100 (MQ=255)
                                    gCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGc                                                                <  1:464706/100‑1 (MQ=255)
                                        ncgatacGCCCTTTATACCCAGTGATGTATTAGCGCTGACGGAGAAGCCCACCACGCTGGCGCCCAGGCGGTGCATAGATGTTCCGAAAGCGAGGc                                                                >  2:872796/7‑96 (MQ=255)
                                          gTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTAGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGcgggtac                                                          >  2:1199229/1‑100 (MQ=255)
                                              cgtccTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGcgggtacgggt                                                      <  1:446246/97‑1 (MQ=255)
                                              cGCCCTTTTTACCAAGTGATGGATAGGCGCTGTCGGAGAAGCCCATCTCGCTGGCCCCCAGGCGGTGCATAGGTGTTTCGAAAGAGAGGcgggtacgggt                                                      <  2:171574/100‑1 (MQ=255)
                                               gCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGcgggtacgggta                                                     >  2:1260942/1‑100 (MQ=255)
                                               gCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGcgggtacggg                                                       <  2:1357717/98‑1 (MQ=255)
                                                 ccTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTaga                                                   <  1:2139738/100‑1 (MQ=255)
                                                  cTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTaaag                                                  >  1:253159/1‑100 (MQ=255)
                                                  cTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTagag                                                  <  1:969659/100‑1 (MQ=255)
                                                  cTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACAACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTagag                                                  <  2:660317/100‑1 (MQ=255)
                                                   tttttACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAgg                                                 <  2:1847060/100‑1 (MQ=255)
                                                       tACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTCCGGGTAGAGGCtt                                              >  1:1382031/1‑99 (MQ=255)
                                                          ccAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCga                                           <  1:950089/99‑1 (MQ=255)
                                                           cAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCgaa                                          >  2:2195725/1‑99 (MQ=255)
                                                              tGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGCGAGGCGGGGTCGGGTAGTGGCtttggatata                                      >  1:741251/1‑92 (MQ=255)
                                                                 tGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAg                                   >  1:1327664/1‑100 (MQ=255)
                                                                    aTTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAAAAACAACTg                                >  1:788118/1‑100 (MQ=255)
                                                                    aTTGGCGCTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGATACGGGTAGAGGCTTCGAAGACACAACTg                                >  2:547674/1‑100 (MQ=255)
                                                                      tGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGc                              <  1:1106058/100‑1 (MQ=255)
                                                                      tGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGc                              <  1:2234072/100‑1 (MQ=255)
                                                                        gcgcTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCaa                            >  1:2115957/1‑100 (MQ=255)
                                                                        gcgcTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCaa                            >  2:748410/1‑100 (MQ=255)
                                                                         cgcaatCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTGCGAAAGAGAGGCGGGGACGGGGAGAGGGTTCGGAGAAACAGCTGGCaa                            >  2:1328415/6‑99 (MQ=255)
                                                                         cgcTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAAt                           <  1:2050322/100‑1 (MQ=255)
                                                                         cgcTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTCCGGGTCGAGGCTTCGAAGAAACAGCTGGCaa                            >  1:361310/1‑99 (MQ=255)
                                                                         cgcTGTCGGAGAAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAAt                           >  1:1429389/1‑100 (MQ=255)
                                                                             gTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGAGGGGAGGGGTAGAGGGATCGAAGAAATAGCTAGCAATGACt                       >  2:377969/1‑100 (MQ=255)
                                                                                 gagaAGCACACGACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTtgt                   >  2:1177946/1‑100 (MQ=255)
                                                                                  agaAGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTtgtg                  <  1:1736282/100‑1 (MQ=255)
                                                                                   gaAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGc                 <  1:1645686/100‑1 (MQ=255)
                                                                                   gaAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGc                 <  2:1057926/100‑1 (MQ=255)
                                                                                    aaGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGc                 <  1:2116404/99‑1 (MQ=255)
                                                                                    aaGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCt                <  1:1897443/100‑1 (MQ=255)
                                                                                    aaGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCt                <  1:1212377/100‑1 (MQ=255)
                                                                                    aaGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCt                <  1:1681238/100‑1 (MQ=255)
                                                                                    aaGCCCACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCt                >  2:373819/1‑100 (MQ=255)
                                                                                     aacccACCACGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCtt               >  1:1512125/3‑100 (MQ=255)
                                                                                     aGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAACTTGCAATGACTTTGTGCtt               >  1:1956060/1‑100 (MQ=255)
                                                                                     aGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGGACGGGTAGAGGCTTCGCAGAAACAGCTGGCCATGACTTTGTGCtt               >  1:801286/1‑100 (MQ=255)
                                                                                      gCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTc              >  1:1987803/1‑100 (MQ=255)
                                                                                        ccaccaCACTGGCCCCCAGGCGGTGCAAAGATGTTTCGAAAGAGAGGCGGGTATGGGGAAAGGCTTTGAAGAAACAGCTGGGCTTGACTTTGTGCTTCaa            >  2:288490/1‑100 (MQ=255)
                                                                                         caccacGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGTTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAAc           <  2:1189969/100‑1 (MQ=255)
                                                                                         caccacGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGGACGGGGAGAGGCTTCGAAGAAAAAGCTGGCAATGACTGTGTGCTTCaa            >  2:1678795/1‑99 (MQ=255)
                                                                                         caccacGCTGCCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAAc           <  1:865612/100‑1 (MQ=255)
                                                                                          accacGCTTGCCCCCAGGCGATGCATAGATGTTTCGAAAGAGATGCGGGTACGGTTAGAGTCTTCGTAGAAACAGCTTGCAATGACTTTGTGCt                <  2:1902401/94‑1 (MQ=255)
                                                                                          accacGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACa          <  1:1599709/100‑1 (MQ=255)
                                                                                          accacGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGGACGGGGAGAGGCTTCGAAGAAACAGCTTGCAATGACTTTGTGCTTCAACa          >  1:1645135/1‑100 (MQ=255)
                                                                                                cTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTg    >  2:1859154/1‑100 (MQ=255)
                                                                                                  ggCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGt  >  1:115616/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GCATCGACGTAAGTGCTGATAACCGAAATGGTATCGGCCAGCGTTTCGCCCTTTTTACCCAGTGATGTATTGGCGCTGTCGGAGAAGCCCACCACGCTGGCCCCCAGGCGGTGCATAGATGTTTCGAAAGAGAGGCGGGTACGGGTAGAGGCTTCGAAGAAACAGCTGGCAATGACTTTGTGCTTCAACAGCTCTGGT  >  NC_000913/4472091‑4472288

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