Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 85,338 A→C 7.0% intergenic (+26/‑292) leuO → / → ilvI global transcription factor/acetolactate synthase 3 large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091385,3380AC 7.0% 299.2 / 10.6 86intergenic (+26/‑292)leuO/ilvIglobal transcription factor/acetolactate synthase 3 large subunit
Reads supporting (aligned to +/- strand):  ref base A (35/45);  new base C (3/3);  total (38/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GTATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGA  >  NC_000913/85190‑85459
                                                                                                                                                    |                                                                                                                         
gTATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGa                                                                                                                           >  3:569528/1‑149 (MQ=255)
 tATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                          >  3:629487/1‑149 (MQ=255)
  aTTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAAt                                                                                                                         <  3:1070900/149‑1 (MQ=255)
       cGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTg                                                                                                                    <  2:738072/149‑1 (MQ=255)
        gCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                                   <  2:226409/149‑1 (MQ=255)
        gCTGCCGTTAAAACAAAAAAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                                   <  3:1117018/149‑1 (MQ=255)
            ccGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTa                                                                                                               <  4:878083/149‑1 (MQ=255)
                  aaacaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGc                                                                                                         >  2:168003/1‑149 (MQ=255)
                   aacaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                                   >  3:869504/1‑138 (MQ=255)
                   aacaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                                   <  4:869504/138‑1 (MQ=255)
                   aacaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCa                                                                                                        >  3:1466670/1‑149 (MQ=255)
                    acaaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAt                                                                                                       >  3:609204/1‑149 (MQ=255)
                     caaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                                   <  1:1425853/136‑1 (MQ=255)
                     caaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTgg                                                                                                                   >  2:1425853/1‑136 (MQ=255)
                     caaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAtt                                                                                                      >  3:222871/1‑149 (MQ=255)
                      aaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAttt                                                                                                     <  4:1211958/149‑1 (MQ=255)
                      aaaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAttt                                                                                                     <  4:932914/149‑1 (MQ=255)
                       aaacaGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAtttt                                                                                                    <  4:1098018/149‑1 (MQ=255)
                          cagcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTa                                                                                                 <  3:1089892/149‑1 (MQ=255)
                          cagcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTa                                                                                                 <  4:222871/149‑1 (MQ=255)
                           agcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTAt                                                                                                >  4:1186945/1‑149 (MQ=255)
                            gcagAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTAtt                                                                                               <  4:569528/149‑1 (MQ=255)
                               gAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGaa                                                                                            <  4:629487/149‑1 (MQ=255)
                                  ccTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAAt                                                                                                                         <  2:1063179/117‑1 (MQ=255)
                                  ccTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAAt                                                                                                                         >  1:1063179/1‑117 (MQ=255)
                                   cTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCa                                                                                                        <  4:586545/133‑1 (MQ=255)
                                   cTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCa                                                                                                        >  3:586545/1‑133 (MQ=255)
                                   cTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTa                                                                                        <  4:609204/149‑1 (MQ=255)
                                    tGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTAtt                                                                                               <  3:262702/141‑1 (MQ=255)
                                    tGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTAtt                                                                                               >  4:262702/1‑141 (MQ=255)
                                      ttATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGaa                                                                                     <  4:1466670/149‑1 (MQ=255)
                                        aTCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                          >  2:570789/1‑110 (MQ=255)
                                        aTCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGaa                                                                                                                          <  1:570789/110‑1 (MQ=255)
                                        aTCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt                                                                                                             <  1:151396/123‑1 (MQ=255)
                                        aTCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt                                                                                                             >  2:151396/1‑123 (MQ=255)
                                             tcCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTAtt                                                                              <  1:168003/149‑1 (MQ=255)
                                                    aTGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAttt                                                                                                     >  4:686612/1‑119 (MQ=255)
                                                    aTGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCAttt                                                                                                     <  3:686612/119‑1 (MQ=255)
                                                    aTGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTa                                                                                >  2:433858/1‑140 (MQ=255)
                                                    aTGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTa                                                                                <  1:433858/140‑1 (MQ=255)
                                                      gAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGa                                                                          <  3:977157/144‑1 (MQ=255)
                                                      gAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGa                                                                          >  4:977157/1‑144 (MQ=255)
                                                                   gcgATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTggg                                                                 >  4:406195/1‑140 (MQ=255)
                                                                   gcgATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTggg                                                                 <  3:406195/140‑1 (MQ=255)
                                                                   gcgATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTAtt                                                        <  3:1186945/149‑1 (MQ=255)
                                                                     gATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt                                                                                                             >  2:1208716/1‑94 (MQ=255)
                                                                     gATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTAtt                                                                                                             <  1:1208716/94‑1 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGaa                                                                         >  2:1025522/1‑129 (MQ=255)
                                                                      aTAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGaa                                                                         <  1:1025522/129‑1 (MQ=255)
                                                                                   gTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTa                                                                                        >  3:121882/1‑101 (MQ=255)
                                                                                   gTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTa                                                                                        <  4:121882/101‑1 (MQ=255)
                                                                                      gatggaAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGattaa                                              <  2:107289/140‑1 (MQ=255)
                                                                                      gatggaAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGattaa                                              >  1:107289/1‑140 (MQ=255)
                                                                                             gagCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCa                                                                                                        >  4:997854/1‑75 (MQ=255)
                                                                                             gagCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCa                                                                                                        <  3:997854/75‑1 (MQ=255)
                                                                                                  aTTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGtt                                                                                                                 >  4:1440707/1‑61 (MQ=255)
                                                                                                  aTTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGtt                                                                                                                 <  3:1440707/61‑1 (MQ=255)
                                                                                                       tctcAATTTGCTAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTaaaaaaa                                     <  2:299855/132‑1 (MQ=255)
                                                                                                       tctcAATTTGCCAACGCTATCTGATTGCAGAATAGGTCAGACATGAATGTCTGATTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTAt                                                         >  2:1056461/1‑112 (MQ=255)
                                                                                                       tctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTAt                                                         <  1:1056461/112‑1 (MQ=255)
                                                                                                       tctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTaaaaaaa                                     >  1:299855/1‑132 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGatt                                                <  3:1341206/120‑1 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGatt                                                >  4:1341206/1‑120 (MQ=255)
                                                                                                        ctcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAg                   >  3:1256883/1‑149 (MQ=255)
                                                                                                         tcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTg                                                                                                                    >  1:1406359/1‑51 (MQ=255)
                                                                                                         tcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTg                                                                                                                    <  2:1406359/51‑1 (MQ=255)
                                                                                                         tcAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGt                  >  3:978709/1‑149 (MQ=255)
                                                                                                           aaTTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTc                                                                                                                      >  1:1110019/1‑47 (MQ=255)
                                                                                                           aaTTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTc                                                                                                                      <  2:1110019/47‑1 (MQ=255)
                                                                                                                         aaCTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGa  <  4:978709/149‑1 (MQ=255)
                                                                                                                           cTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGt                      >  4:605231/1‑127 (MQ=255)
                                                                                                                           cTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGt                      <  3:605231/127‑1 (MQ=255)
                                                                                                                            tGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGa                                                  <  2:65263/98‑1 (MQ=255)
                                                                                                                            tGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGa                                                  >  1:65263/1‑98 (MQ=255)
                                                                                                                                     aaTAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtgt                                                                    >  1:314156/1‑71 (MQ=255)
                                                                                                                                     aaTAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtgt                                                                    <  2:314156/71‑1 (MQ=255)
                                                                                                                                           ccAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTa                                                     <  3:730939/79‑1 (MQ=255)
                                                                                                                                           ccAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTa                                                     >  4:730939/2‑80 (MQ=255)
                                                                                                                                           ccAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAgg                                                   <  2:282092/81‑1 (MQ=255)
                                                                                                                                           ccAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATGTTATTCTGAATGTGTGGGCTCTCTATTTTAgg                                                   >  1:282092/2‑82 (MQ=255)
                                                                                                                                            cAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTaaaaaaa                                     >  4:799845/1‑95 (MQ=255)
                                                                                                                                            cAGACATTCATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTaaaaaaa                                     <  3:799845/95‑1 (MQ=255)
                                                                                                                                               aCATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGattaa                                              >  4:436641/1‑83 (MQ=255)
                                                                                                                                               aCATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGattaa                                              <  3:436641/83‑1 (MQ=255)
                                                                                                                                               aCATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGatt                                                <  1:380071/81‑1 (MQ=255)
                                                                                                                                               aCATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGatt                                                >  2:380071/1‑81 (MQ=255)
                                                                                                                                               aCATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTaaaaaaa                                     <  4:530009/92‑1 (MQ=255)
                                                                                                                                               aCATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTaaaaaaa                                     >  3:530009/1‑92 (MQ=255)
                                                                                                                                                 aTGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtg                                                                     <  4:747630/58‑1 (MQ=255)
                                                                                                                                                 aTGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAAtgtg                                                                     >  3:747630/1‑58 (MQ=255)
                                                                                                                                                    aaTGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATa                                   <  4:1410706/89‑1 (MQ=255)
                                                                                                                                                    aaTGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATa                                   >  3:1410706/1‑89 (MQ=255)
                                                                                                                                                    |                                                                                                                         
GTATTACCGCTGCCGTTAAAACAAAACAGCAGAACCTGTTATCTCTCCTGGCATGAAGCTGCCGGGCGCGATAAAGGCCATCAGTGGATGGAAGAGCAATTAGTCTCAATTTGCAAACGCTAACTGATTGCAGAATAGGTCAGACATGAATGTCTGGTTTATTCTGCATTTTTTATTGAATGTAGAATTTTATTCTGAATGTGTGGGCTCTCTATTTTAGGATTAATTAAAAAAATAGAGAAATTGCTGTAAGTTGTGGGATTCAGCCGA  >  NC_000913/85190‑85459

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