Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,606,094 (G)8→7 noncoding (72/87 nt) leuV ← tRNA‑Leu

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,606,0870G.96.8% 274.9 / 0.2 63noncoding (79/87 nt)leuVtRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base G (0/1);  major base . (33/27);  minor base C (0/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 2.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTT  >  NC_000913/4605940‑4606219
                                                                                                                                                    |                                                                                                                                    
tgCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggg                                                                                                                                   >  1:202878/1‑151 (MQ=255)
    gTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggg                                                                                                                                 <  2:510844/149‑1 (MQ=255)
      aTCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggact                                                                                                                             >  1:768273/1‑148 (MQ=255)
        cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTg                                                                                                                          <  2:769220/151‑1 (MQ=255)
        cacaAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTCATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAgggggggactt                                                                                                                            >  1:338860/1‑146 (MQ=255)
             aCAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAAccc                                                                                                                     >  2:742138/1‑151 (MQ=255)
              cAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAAcccc                                                                                                                    >  1:824264/1‑151 (MQ=255)
              cAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAAcccc                                                                                                                    >  1:255120/1‑151 (MQ=255)
                gTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCa                                                                                                                  <  2:137994/151‑1 (MQ=255)
                  tGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCAc                                                                                                                 <  2:315857/150‑1 (MQ=255)
                  tGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACg                                                                                                                >  1:281160/1‑151 (MQ=255)
                    cGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTc                                                                                                              >  1:374586/1‑151 (MQ=255)
                    cGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTc                                                                                                              <  2:36875/151‑1 (MQ=255)
                     ggCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTcc                                                                                                             <  2:253169/151‑1 (MQ=255)
                     ggCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTcc                                                                                                             >  2:332922/1‑151 (MQ=255)
                           aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAgg                                                                                                       >  1:834952/1‑151 (MQ=255)
                           aCATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAg                                                                                                        <  1:252512/150‑1 (MQ=255)
                             aTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                     <  1:389345/151‑1 (MQ=255)
                             aTTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGac                                                                                                     <  1:738071/151‑1 (MQ=255)
                                   cTACAATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAac                                                                                               >  1:660001/1‑151 (MQ=255)
                                       cATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCt                                                                                           >  1:649232/1‑151 (MQ=255)
                                                 ggCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAg                                                                                   <  1:629244/149‑1 (MQ=255)
                                                  gcagaaaTTTCTGGCATCATGCTAATGATAATGAATTTGAAGCAGGGACCAAAAACAAGAGAAGAAACAACAAAAGAAGTAAAGACAATTAGTTCGACGGAGGGGACTTGAACCCCCACGTCCGTAAGGACAATAACGCCTGAAGCTa                                                                                    <  1:779636/142‑1 (MQ=19)
                                                     cTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTc                                                                             >  1:89782/1‑151 (MQ=255)
                                                        aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                          <  2:438225/151‑1 (MQ=255)
                                                        aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                          >  2:268832/1‑151 (MQ=255)
                                                        aTTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTAc                                                                          <  1:44910/151‑1 (MQ=255)
                                                              ggCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAAtt                                                                     <  2:733305/150‑1 (MQ=255)
                                                                    atGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCAc                                                              >  2:672033/1‑151 (MQ=255)
                                                                         gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTc                                                          <  1:596837/150‑1 (MQ=255)
                                                                         gatgatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCATCTACCAATTCCGCCACCtt                                                           <  1:412839/149‑1 (MQ=255)
                                                                            gatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                      >  2:803890/1‑151 (MQ=255)
                                                                            gatAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGca                                                       <  1:460757/150‑1 (MQ=255)
                                                                            gatAATGAACTTGAATCGGGGAGGAAGAACAGAGGACTTACCAATAAAGGAATCAAAGAGAATGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGAACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCtt                                                           <  1:107699/146‑1 (MQ=11)
                                                                             atAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGcac                                                      >  2:637329/1‑150 (MQ=255)
                                                                               aaTGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAg                                                    <  1:742095/150‑1 (MQ=255)
                                                                                    aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                              <  2:649274/151‑1 (MQ=255)
                                                                                    aTTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCt                                              <  2:429932/151‑1 (MQ=255)
                                                                                       tGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTa                                            <  1:122161/150‑1 (MQ=255)
                                                                                        gAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTGATCtt                                             >  2:127891/1‑148 (MQ=255)
                                                                                         aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtt                                         >  2:151726/1‑151 (MQ=255)
                                                                                         aaTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACtt                                         >  1:20887/1‑151 (MQ=255)
                                                                                          aTCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttt                                        >  2:481102/1‑151 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt                                      <  2:116812/151‑1 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt                                      >  2:346958/1‑151 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttttt                                      <  1:505680/151‑1 (MQ=255)
                                                                                            cAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACttt                                        >  1:637491/1‑149 (MQ=255)
                                                                                                 aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATATTACTTTTTTTGa                                  >  2:114534/1‑150 (MQ=255)
                                                                                                 aTGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGCCATCTTACTTTTATTGa                                  >  1:465766/1‑150 (MQ=255)
                                                                                                   gAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGatat                               <  1:729064/151‑1 (MQ=255)
                                                                                                    aaaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATc                              <  2:703979/151‑1 (MQ=255)
                                                                                                     aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGACGGTGGGTACTTGAATCCCCACGTCCGTAAGCACACTAACACCTGACGCTATCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGata                                <  1:13054/149‑1 (MQ=25)
                                                                                                     aaaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCg                             >  1:825608/1‑151 (MQ=255)
                                                                                                      aaaTCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGatat                               >  1:363987/1‑148 (MQ=255)
                                                                                                          cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCTCAATACtttttaagatttcctatc                         >  1:168820/1‑137 (MQ=18)
                                                                                                          cAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTc                         <  1:757846/150‑1 (MQ=255)
                                                                                                                 agaaACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGt                  <  2:617780/151‑1 (MQ=255)
                                                                                                                     aCAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCg               >  2:266228/1‑150 (MQ=255)
                                                                                                                       aaGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAg             >  2:473115/1‑150 (MQ=255)
                                                                                                                         gAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggg          >  2:731581/1‑151 (MQ=255)
                                                                                                                          aaaGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAggggg         <  2:530878/151‑1 (MQ=255)
                                                                                                                           aaGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAgggggg        <  1:305473/151‑1 (MQ=255)
                                                                                                                            aGGAAGTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTATTGATATCGCCTCGCTTGGTGCGAggcggg        >  1:559579/1‑150 (MQ=255)
                                                                                                                                 gTAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACtt  >  2:33578/1‑151 (MQ=255)
                                                                                                                                  tAAAGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACtt  >  2:266080/1‑150 (MQ=255)
                                                                                                                                    aaGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTTCTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACtt  >  2:589428/1‑148 (MQ=255)
                                                                                                                                    aaGATAATTGGTGCGAGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTATCCTTACTTTAAATGAAATCCTATCGCTGGGGGCGAGTCGGGGGACtt  >  1:555518/1‑148 (MQ=14)
                                                                                                                                                    |                                                                                                                                    
TGCTGTATCACAAACAGTTGCGGCATTACATTTTGCTACCATCGTGCCTGGCACTTATTGCTGGCATCATGCTGATGATAATGAATTTGAATCAGGGATGAAAAATCAAGGAAGAAACAAGAAAGGAAGTAAAGATAATTGGTGCGAGGGGGGGGACTTGAACCCCCACGTCCGTAAGGACACTAACACCTGAAGCTAGCGCGTCTACCAATTCCGCCACCTTCGCACAGTCATCTTACTTTTTTTGATATCGCCTCGTTTGGTGCGAGGGGGGGGACTT  >  NC_000913/4605940‑4606219

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