Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,074,945 T→G 6.0% N1313T (AAC→ACC)  putA ← fused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,074,9450TG 6.0% 300.2 / 9.7 83N1313T (AAC→ACC) putAfused DNA‑binding transcriptional regulator/proline dehydrogenase/pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (41/37);  new base G (3/2);  total (44/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCT  >  NC_000913/1074799‑1075092
                                                                                                                                                  |                                                                                                                                                   
gACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTa                                                                                                                                                   <  4:360705/149‑1 (MQ=255)
           aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                                                                                                     >  1:908379/1‑136 (MQ=255)
           aTGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGt                                                                                                                                                     <  2:908379/136‑1 (MQ=255)
                                   ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                                                                                                 <  4:60783/116‑1 (MQ=255)
                                   ttGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAcc                                                                                                                                                 >  3:60783/1‑116 (MQ=255)
                                     gAAAACGTGCGGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCt                                                                                                              >  3:89396/1‑149 (MQ=255)
                                                                gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                                                         <  2:240087/143‑1 (MQ=255)
                                                                gggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                                                         >  1:240087/1‑143 (MQ=255)
                                                                 ggATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt                                                                                  <  4:89396/149‑1 (MQ=255)
                                                                       gACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGc                                                                                                                                               <  3:582335/82‑1 (MQ=255)
                                                                       gACCACGTTAAAGACGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCAAACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTc                                                                            >  3:100581/1‑149 (MQ=255)
                                                                       gACCACGTTAAACATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGc                                                                                                                                               >  4:582335/1‑82 (MQ=255)
                                                                         ccACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGc                                                                          >  2:104733/1‑149 (MQ=255)
                                                                           aCGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCa                                                                        >  1:708150/1‑149 (MQ=255)
                                                                                                gTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaa                                                   >  2:435367/1‑149 (MQ=255)
                                                                                                 ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                                                                  <  2:1060133/53‑1 (MQ=255)
                                                                                                 ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTAc                                                                                                                                                  >  1:1060133/1‑53 (MQ=255)
                                                                                                 ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                                                                    <  2:1102167/115‑1 (MQ=255)
                                                                                                 ttGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGata                                                                                    >  1:1102167/1‑115 (MQ=255)
                                                                                                    tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                                                                     <  1:789841/127‑1 (MQ=255)
                                                                                                    tAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                                                                     >  2:789841/1‑127 (MQ=255)
                                                                                                      aCATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAAtt                                             <  1:435367/149‑1 (MQ=255)
                                                                                                        aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt                                           <  2:708150/149‑1 (MQ=255)
                                                                                                        aTCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt                                           >  2:351685/1‑149 (MQ=255)
                                                                                                         tcctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaaa                                                               >  3:498978/1‑128 (MQ=255)
                                                                                                         tcctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCaaaa                                                               <  4:498978/128‑1 (MQ=255)
                                                                                                           ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTAtt                                                                                                                               >  1:1464671/1‑62 (MQ=255)
                                                                                                           ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTAtt                                                                                                                               <  2:1464671/62‑1 (MQ=255)
                                                                                                           ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTca                                                                                                                             >  3:1341418/1‑64 (MQ=255)
                                                                                                           ctccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTca                                                                                                                             <  4:1341418/64‑1 (MQ=255)
                                                                                                            tccGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCa                                       >  3:531858/1‑149 (MQ=255)
                                                                                                              cGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATc                                     >  3:301880/1‑149 (MQ=255)
                                                                                                               ggCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGa                                                                                                           <  1:296331/78‑1 (MQ=255)
                                                                                                               ggCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGa                                                                                                           >  2:296331/1‑78 (MQ=255)
                                                                                                                  tACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaaa                                                  <  2:318285/132‑1 (MQ=255)
                                                                                                                  tACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGaaa                                                  >  1:318285/1‑132 (MQ=255)
                                                                                                                  tACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCcgc                                 >  4:1218937/1‑149 (MQ=255)
                                                                                                                      tGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCgg                             <  1:351685/149‑1 (MQ=255)
                                                                                                                      tGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCgg                             >  3:1206453/1‑149 (MQ=255)
                                                                                                                       gTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                                                                     >  3:309192/1‑108 (MQ=255)
                                                                                                                       gTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACgg                                                                     <  4:309192/108‑1 (MQ=255)
                                                                                                                        tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAAc                                                                                                                 >  1:192253/1‑63 (MQ=255)
                                                                                                                        tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAAc                                                                                                                 <  2:192253/63‑1 (MQ=255)
                                                                                                                        tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                             <  3:1215411/83‑1 (MQ=255)
                                                                                                                        tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                             >  4:1215411/1‑83 (MQ=255)
                                                                                                                        tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCAc                                                       >  1:66082/1‑121 (MQ=255)
                                                                                                                        tttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCAc                                                       <  2:66082/121‑1 (MQ=255)
                                                                                                                         ttAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          <  4:301880/149‑1 (MQ=255)
                                                                                                                           aaCCTATAGTTATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          >  1:227210/1‑147 (MQ=255)
                                                                                                                           aaCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          <  2:227210/147‑1 (MQ=255)
                                                                                                                           aaCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          <  4:514242/147‑1 (MQ=255)
                                                                                                                           aaCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCaa                          >  3:514242/1‑147 (MQ=255)
                                                                                                                            aCCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTg                       >  1:1335621/1‑149 (MQ=255)
                                                                                                                             ccTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGc                      <  3:1218937/149‑1 (MQ=255)
                                                                                                                              cTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCt                     >  2:438212/1‑149 (MQ=255)
                                                                                                                               tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                             >  3:7490/1‑76 (MQ=255)
                                                                                                                               tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTcc                                                                                             <  4:7490/76‑1 (MQ=255)
                                                                                                                               tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCtt                    <  2:1335621/149‑1 (MQ=255)
                                                                                                                               tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCtt                    <  1:438212/149‑1 (MQ=255)
                                                                                                                               tataGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCAAACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGTATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCtt                    <  4:100581/149‑1 (MQ=255)
                                                                                                                                  aGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTcac                 >  2:992324/1‑149 (MQ=255)
                                                                                                                                   gTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTcaca                <  1:104733/149‑1 (MQ=255)
                                                                                                                                   gTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTcaca                >  3:1459766/1‑149 (MQ=255)
                                                                                                                                    tCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                                                         <  4:1204786/75‑1 (MQ=255)
                                                                                                                                    tCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGaa                                                                                         >  3:1204786/1‑75 (MQ=255)
                                                                                                                                      aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGc                      >  3:1271557/1‑140 (MQ=255)
                                                                                                                                      aTTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGc                      <  4:1271557/140‑1 (MQ=255)
                                                                                                                                       ttAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAg                                                                                        <  3:1302316/73‑1 (MQ=255)
                                                                                                                                       ttAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAg                                                                                        >  4:1302316/1‑73 (MQ=255)
                                                                                                                                        tAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATACTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGcc                                        >  4:1041312/1‑120 (MQ=255)
                                                                                                                                        tAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATACTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGcc                                        <  3:1041312/120‑1 (MQ=255)
                                                                                                                                         aaGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt                                                                                  >  4:1383638/1‑77 (MQ=255)
                                                                                                                                         aaGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt                                                                                  <  1:801102/77‑1 (MQ=255)
                                                                                                                                         aaGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt                                                                                  <  3:1383638/77‑1 (MQ=255)
                                                                                                                                         aaGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATAtt                                                                                  >  2:801102/1‑77 (MQ=255)
                                                                                                                                          cgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTc                   >  4:1314519/2‑139 (MQ=255)
                                                                                                                                          cgCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTc                   <  3:1314519/138‑1 (MQ=255)
                                                                                                                                          aGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATgcg         >  3:110605/1‑149 (MQ=255)
                                                                                                                                           gCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt                                           <  1:167093/114‑1 (MQ=255)
                                                                                                                                           gCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATtgt                                           >  2:167093/1‑114 (MQ=255)
                                                                                                                                           gCTGGCGGTAACGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATgcgc        >  1:480871/1‑149 (MQ=255)
                                                                                                                                              ggCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGCACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGaa     <  4:110605/149‑1 (MQ=255)
                                                                                                                                                cGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAAtt                                             <  4:1138113/107‑1 (MQ=255)
                                                                                                                                                cGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAAtt                                             >  3:1138113/1‑107 (MQ=255)
                                                                                                                                                 gTTACCGCCAGCGGCAGCGGTATTCACACTCAGTGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCt  >  2:1028975/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GACGTTGCGCGATGAGGTATTTTTCAATCAAACGGTTGAAAACGTGCAGCGGATTATTATTGAGGGGATTCGACCACGTTAAAGATGCCGGAGGAGGTTGTAACATCCTCCGGCTACCTGTTTAACCTATAGTCATTAAGCTGGCGTTACCGCCAGCGGCAGCGGTATTCACACTCAGCGAACGCTCGATATACAGCCGTTCCAGAAGGATATTGCTTTCGCCACGGGCAAAACCCTGCACCGAAACAATTGTGCCATCCCGCGCGGCAACTGCTTCACACAATGCGCGAAGCT  >  NC_000913/1074799‑1075092

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