Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,200,901 A→G *466W (TGA→TGG) 
T2A (ACG→GCG) 		 hydH →
hydG →		 HydH
HydG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,200,9010AG99.7% 987.5 / NA 282*466W (TGA→TGG) 
T2A (ACG→GCG) 		hydH
hydG		HydH
HydG
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (141/140);  total (141/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

TACAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAAC  >  NC_000913/4200811‑4200990
                                                                                          |                                                                                         
 aCAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGgcgc                                                                                        <  2:3775052/93‑1 (MQ=255)
 aCAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGgcgc                                                                                        >  2:1700123/1‑93 (MQ=255)
 aCAATTCAGGGCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGgcgc                                                                                        <  1:5369079/93‑1 (MQ=255)
  cAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCa                                                                                       >  1:4569685/1‑93 (MQ=255)
   aaTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCAc                                                                                      <  1:3930561/93‑1 (MQ=255)
    aTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACg                                                                                     >  2:1833052/1‑93 (MQ=255)
    aTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACg                                                                                     >  2:1362164/1‑93 (MQ=255)
    aTTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACg                                                                                     >  2:4885657/1‑93 (MQ=255)
     ttCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGa                                                                                    >  2:5725762/1‑93 (MQ=255)
     ttCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGa                                                                                    >  2:1147146/1‑93 (MQ=255)
        aGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataa                                                                                 >  2:2395135/1‑93 (MQ=255)
        aGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataa                                                                                 >  2:5589609/1‑93 (MQ=255)
         ggTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACTataat                                                                                <  1:4286210/93‑1 (MQ=255)
         ggTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGacaat                                                                                >  2:4750234/1‑93 (MQ=255)
          gTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataata                                                                               <  2:2247318/93‑1 (MQ=255)
          gTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGataata                                                                               >  1:4618540/1‑93 (MQ=255)
            cGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATc                                                                             <  2:4026394/93‑1 (MQ=255)
                aGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGata                                                                         <  1:1288821/93‑1 (MQ=255)
                  ccAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATAtt                                                                       <  1:2348581/93‑1 (MQ=255)
                   cAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTc                                                                      <  1:3333937/93‑1 (MQ=255)
                   cAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTc                                                                      >  2:5553452/1‑93 (MQ=255)
                   cAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTc                                                                      >  1:5160049/1‑93 (MQ=255)
                    aggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCt                                                                     <  2:2923008/93‑1 (MQ=255)
                    aggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCt                                                                     <  1:4477876/93‑1 (MQ=255)
                     ggaggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtg                                                                    >  2:4647293/1‑93 (MQ=255)
                       aggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggt                                                                  <  1:3766942/93‑1 (MQ=255)
                       aggGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggt                                                                  >  2:4346753/1‑93 (MQ=255)
                        gggAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtg                                                                 <  2:2769980/93‑1 (MQ=255)
                        gggAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtg                                                                 >  2:4264250/1‑93 (MQ=255)
                        gggAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtg                                                                 <  1:1968668/93‑1 (MQ=255)
                          gAAAAGGCTCAGCGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggt                                                               <  1:607538/93‑1 (MQ=255)
                          gAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggt                                                               >  2:2224994/1‑93 (MQ=255)
                           aaaaGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggtg                                                              >  2:4058471/1‑93 (MQ=255)
                            aaaGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCtggtggtgg                                                             >  1:2971757/1‑93 (MQ=255)
                             agggCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGga                                                            <  1:15410/91‑1 (MQ=255)
                             aaGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGga                                                            >  2:5695014/1‑93 (MQ=255)
                              aGGCTCAACGTTCGCCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgat                                                           <  2:5998178/93‑1 (MQ=255)
                              aGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgat                                                           <  2:6347041/93‑1 (MQ=255)
                              aGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgat                                                           <  2:2669122/93‑1 (MQ=255)
                              aGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgag                                                           >  1:797159/1‑92 (MQ=255)
                               ggCTCAACGTTCACCCTCTGGCTCCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgatg                                                          <  1:166388/93‑1 (MQ=255)
                               ggCTCAACGTGCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgatg                                                          <  2:2971760/93‑1 (MQ=255)
                                 cTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGgatgat                                                        <  1:271876/93‑1 (MQ=255)
                                    aaCGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGAc                                                     <  1:4366168/93‑1 (MQ=255)
                                     aCGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACa                                                    <  2:5300207/93‑1 (MQ=255)
                                      cGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACAt                                                   <  2:840550/93‑1 (MQ=255)
                                      cGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACAt                                                   <  2:186979/93‑1 (MQ=255)
                                       gTTCACCCTCTGGCTTCCGGTCACTATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACAtt                                                  >  2:5089301/1‑93 (MQ=255)
                                        ttCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTa                                                 <  1:5802129/93‑1 (MQ=255)
                                        ttCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTa                                                 <  2:5874116/93‑1 (MQ=255)
                                         tCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAg                                                <  2:614527/93‑1 (MQ=255)
                                           tcccTCTGGCTTCCGGTCGATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGcc                                              <  1:4972671/92‑1 (MQ=255)
                                             ccTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCAc                                            <  1:3847127/93‑1 (MQ=255)
                                             ccTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCAc                                            >  1:3290769/1‑93 (MQ=255)
                                              ctctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACt                                           >  2:1754141/1‑93 (MQ=255)
                                              ctctGGCTCCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACt                                           <  1:4622759/93‑1 (MQ=255)
                                                ctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGc                                         >  1:4427875/1‑93 (MQ=255)
                                                ctGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGc                                         >  2:5634889/1‑93 (MQ=255)
                                                   gCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACt                                      >  1:4066310/1‑93 (MQ=255)
                                                     ttCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAt                                    >  1:6338075/1‑93 (MQ=255)
                                                      tCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAtt                                   <  2:4687637/93‑1 (MQ=255)
                                                      tCCGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAtt                                   >  2:1317217/1‑93 (MQ=255)
                                                        cGGTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTAtttt                                 <  2:5489560/93‑1 (MQ=255)
                                                         ggTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTg                                >  2:3390010/1‑93 (MQ=255)
                                                          gTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGc                               >  2:1090484/1‑93 (MQ=255)
                                                          gTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGc                               <  1:778562/93‑1 (MQ=255)
                                                          gTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGCGGATGATGACATTCGCCACTGCACTATTTTGc                               >  1:2088126/1‑93 (MQ=255)
                                                          gTCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATACTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGc                               >  1:2512983/1‑93 (MQ=255)
                                                           tCAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCa                              <  2:5379900/93‑1 (MQ=255)
                                                            cAATATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAg                             >  1:523406/1‑93 (MQ=255)
                                                             aaTATTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAgg                            >  1:2579244/1‑93 (MQ=255)
                                                              atatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGc                           >  1:6165417/1‑93 (MQ=255)
                                                              atatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGc                           >  1:4991169/1‑93 (MQ=255)
                                                              atatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGc                           <  2:6218506/93‑1 (MQ=255)
                                                               tatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATGTCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCt                          <  1:1462448/93‑1 (MQ=255)
                                                               tatTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCt                          <  2:6094221/93‑1 (MQ=255)
                                                                atTACGCGTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCtt                         >  1:4699257/1‑93 (MQ=255)
                                                                  taagcgTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTa                       <  1:3227792/90‑1 (MQ=255)
                                                                       gTAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcg                  <  2:4467696/93‑1 (MQ=255)
                                                                        tAAGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTgcgc                 >  1:1346114/1‑93 (MQ=255)
                                                                          aGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCgg               >  2:3300226/1‑93 (MQ=255)
                                                                          aGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCgg               >  2:519702/1‑93 (MQ=255)
                                                                          aGGACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCgg               >  1:6437693/1‑93 (MQ=255)
                                                                           ggACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGGTGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGc              >  2:3675467/1‑93 (MQ=255)
                                                                           ggACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGc              >  1:3679789/1‑93 (MQ=255)
                                                                           ggACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGc              >  2:4561673/1‑93 (MQ=255)
                                                                           ggACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGc              <  2:749379/93‑1 (MQ=255)
                                                                           ggACCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGc              >  1:5874787/1‑93 (MQ=255)
                                                                             aCCCACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTg            <  1:1619226/93‑1 (MQ=255)
                                                                               ccACAAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTggg          >  1:3590188/1‑93 (MQ=255)
                                                                                cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         <  1:3076939/93‑1 (MQ=255)
                                                                                cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         >  2:3112083/1‑93 (MQ=255)
                                                                                cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         >  2:3537556/1‑93 (MQ=255)
                                                                                cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         >  1:5602773/1‑93 (MQ=255)
                                                                                cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         >  2:4506097/1‑93 (MQ=255)
                                                                                cacaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTgggg         >  1:4210578/1‑93 (MQ=255)
                                                                                 acaAGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGc        >  1:4791022/1‑93 (MQ=255)
                                                                                  aaaGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCt       <  2:2031272/92‑1 (MQ=255)
                                                                                    aGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCtat     >  1:1835514/1‑93 (MQ=255)
                                                                                    aGGATGGCGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCtat     >  2:4962841/1‑93 (MQ=255)
                                                                                          |                                                                                         
TACAATTCAGGTCGCAAGCCAGGAGGGAAAAGGCTCAACGTTCACCCTCTGGCTTCCGGTCAATATTACGCGTAAGGACCCACAAGGATGACGCACGATAATATCGATATTCTGGTGGTGGATGATGACATTAGCCACTGCACTATTTTGCAGGCTTTACTGCGCGGCTGGGGCTATAAC  >  NC_000913/4200811‑4200990

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